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Introduction 



The study of the behaviour of the hadronic interactions at colhders is a fascinating subject and a 
speciahzed field of research that is essential in order to widen our knowledge on the fundamental 
constituents of matter. 

Elementary particle physics, needless to say, is both a theoretical and an experimental science, 
and requires a continuous interplay between the two approaches. Therefore it is important, from 
this viewpoint, to develope, from the theorist's side, theoretical tools and strategies that can be 
used by the experimental community, and this requires considerable effort. 

This requires the investigation of specific processes which can be tested by experiments, 
developing formalism that can be used for further theoretical elaborations, but this also points 
toward the need to develope software that can help the experimental collaborations to proceed 
with their complex phenomenological analysis. 

This thesis is the result of this philosophy, and collects several applications of Quantum 
Chromodynamics, the accepted theory of the strong interactions, to hadron colliders both at high 
and intermediate energies, and on the study of a specific process, termed by us "deeply virtual 
neutrino scattering", where we generalize the formaUsm of deeply virtual Compton scattering to 
neutral currents. The process is potentially relevant for neutrino detection at neutrino factories. 

This dissertation is composed of two parts. In the first section we focus our attention on 
the study of the initial state scaling violations and the evolution of the unpolarized parton 
distributions through Next-to-Next to Leading Order (NNLO) in a^, the strong coupling, suitable 
for precision studies of the parton model at the LHC. Specifically, we will analize the methods 
available to solve the equations and develope the theory that underlines a new proposed method 
which is shown to be highly accurate through NNLO. The theoretical analysis is accompanied 
by the developement of professional software partly documented in this thesis. 

The motivations of this study are in the need to compare traditional strategies in the solutions 
of the DGLAP equations with an independent theoretical approach which allows to reorganize 
all the scaling violations at a fixed perturbative order in terms of logarithms of the coupling 
constant times some scale invariant functions, introduced via factorization. 

We recall that high energy collisions can be well understood by the use of the factorization 
theorems, which allow to separate the perturbative information contained in the hard scatter- 
ings of a given process from the non-perturbative information which is contained in the parton 
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distribution functions. These tell us what is the probability density for finding a given parton 
with a fraction x of the total initial momentum of the incoming hadron at an energy Q, and 
have a formal definition as light-cone correlation functions. 

The knowledge of the PDFs with higher and higher accuracy is an important task for future 
experiments, for the correct interpretation of the experimental data at the LHC, but they are 
also essential for a better description of the internal structure of the hadrons. 

In the first chapter of the thesis we discuss the solutions of the NNLO DGLAP equations using 
x-space methods and for this purpose we show how our approach allows to obtained accurate 
and exact solutions of the evolution equations for the pdf's using an analytical ansatz. We have 
applied the formal developements contained in this analysis is a numerical program, XSIEVE 
written in collaboration with A. Cafarella and C. Coriano that has been used for precision studies 
of two specific processes at hadron colliders. 

The first application has been in the NLO study of the double spin asymmetries in DY proton- 
antiproton collisions, near the kinematical region of the J/^p resonance. This study, relevant for 
the proposed PAX experiment on the measurement of the transverse spin distribution, is analized 
in detail. We present result for the asymmetries of the process and elaborate on their possible 
measurements in the near future. A second application that we discuss is the study of the total 
cross section for Higgs production at the Tevatron and at the LHC, focusing in particular on the 
renormalization/factorization scale dependence of the predictions. While in other applications 
presented in the previous literature this study has been more limited, we perform a numerical 
analysis of the stability of the predictions on a large range of variability. These studies have been 
performed on a cluster. 

The second section of the thesis will be based on applications of perturbative QCD to exclusive 
processes at intermediate energy in presence of weak interactions. 

In particular we will analyse some theoretical aspects of the Generalized Parton Distribution 
functions (GPDs) which are the natural generalization of the PDFs when perturbative QCD is 
applied to exclusive reactions at intermediate energies. We introduce the notion of electroweak 
GPD's and develope the relevant formalism. Then, we will develope a phenomenological analysis 
of the neutrino-nucleon reaction mediated by a Zq boson exchange and we will calculate the 
leading twist amplitudes for the charged current interaction as well, involving bosons. 

The Relevance of the Perturbative QCD to the Precision Studies of the Physics 
at LHC: Motivations 

Precision studies of some hadronic processes in the perturbative regime are going to be very 
important in order to confirm the validity of the mechanism of mass generation in the Standard 
Model at the new collider, the LHC. This program involves a rather complex analysis of the QCD 
background, with the corresponding radiative corrections taken into account to higher orders. 
Studies of these corrections for specific processes have been performed by various groups, at a level 



of accuracy which has reached the next-to-next-to-leading order (NNLO) in a^, the QCD couphng 
constant. The quantification of the impact of these corrections requires the determination of the 
hard scattering of the partonic cross sections up to order , with the matrix of the anomalous 
dimensions of the DGLAP kernels determined at the same perturbative order. 

The study of the evolution of the parton distributions may include both a NNLO analysis 
and, possibly, a resummation of the large logs which may appear in certain kinematic regions 
of specific processes [1]. The questions that we address in this thesis concern the types of 
approximations which are involved when we try to solve the DGLAP equations to higher orders 
and the differences among the various methods proposed for their solution. 

The clarification of these issues is important, since a chosen method has a direct impact 
on the structure of the evolution codes and on their phenomenological predictions. We address 
these questions by going over a discussions of these methods and, in particular, we compare 
those based in Mellin space and the analogous ones based in x-space. Mellin methods have been 
the most popular and have been implemented up to NLO and, very recently, also at NNLO [2]. 
We remark that j;-space methods based on logarithmic expansions have never been thoroughly 
justified in the previous literature even at NLO, in the case of the QCD parton distributions 

(pdf's) [a HIE]. 

We fill this gap and present exact proofs of the equivalence of these methods - in the case of 
the evolution of the QCD pdf's - extending a proof which had been outlined by Rossi |B] at LO 
and by Da Luz Vieira and Storrow at NLO in their study of the parton distributions of the 
photon. 

In more recent times, these studies on the pdf's of the photon have triggered similar studies 
also for the QCD pdf's. The result of these efforts was the proposal of new expansions for the 
quark and gluon parton distributions [HI HI [5] which had to capture the logarithmic behaviour 
of the solution up to NLO. Evidence of the consistency of the ansatz was, in part, based on a 
comparative study of the generic structures of the logarithms that appear in the solution using 
Mellin moments, since the same logarithms of the coupling could be reobtained by recursion 
relations. 

In this thesis we are going to clarify - using the exact solutions of the corresponding recursion 
relations - the role of these previous expansions and present their generalizations. In particular 
we will show that they can be extended to retain higher logarithmic corrections and how they 
can be made exact. It is shown that by a suitable extension of this analysis, all that has been 
known so far in moment space can be reobtained directly from x-space. In the photon case the 
Da Luz Vieira-Storrow solution [7] can be now understood simply as a first truncated ansatz of 
the general truncated solutions that we analize. 



Applications of the Perturbative Technique to the Neutrino Physics Regime 

There is no doubt that the study of neutrino masses and of flavour mixing in the leptonic 
sector will play a crucial role for uncovering new physics beyond the Standard Model and to 
test Unification. In fact, the recent discovery of neutrino oscillations in atmospheric and solar 
neutrinos (see [8j for an overview) has raised the puzzle of the origin of the mass hierarchy among 
the various neutrino flavours, a mystery which, at the moment, remains unsolved. The study 
of the mixing among the leptons also raises the possibility of detecting possible sources of CP 
violation in this sector as well. It seems then obvious that the study of these aspects of flavour 
physics requires the exploration of a new energy range for neutrino production and detection 
beyond the one which is accessible at this time. 

For this purpose, several proposals have been presented recently for neutrino factories, where 
a beam, primarily made of muon neutrinos produced at an accelerator facility, is directed to a 
large volume located several hundreds kilometers away at a second facility. The goal of these 
experimental efforts is to uncover various possible patterns of mixings among flavours - using the 
large distance between the points at which neutrinos are produced and detected - in order to study 
in a more detailed and "artificial" way the phenomenon of oscillations. Detecting neutrinos at 
this higher energies is an aspect that deserves special attention since several of these experimental 
proposals ttHl E] require a nominal energy of the neutrino beam in the few GeV region. We 
recall that the incoming neutrino beam, scattering off deuteron or other heavier targets at the 
detector facility, has an energy which covers, in the various proposals, both the resonant, the 
quasi-elastic (in the GeV range) and the deep inelastic region (DIS) at higher energy. In the 
past, neutrino scattering on nucleons has been observed over a wide interval of energy, ranging 
from few MeV up to 100 GeV, and these studies have been of significant help for uncovering 
the structure of the fundamental interactions in the Standard Model. Generally, one envisions 
contributions to the scattering cross section either in the low energy region, such as in neutrino- 
nucleon elastic scattering, or in the deep inelastic scattering (DIS) region. Recent developments 
in perturbative QCD have emphasized that exclusive (see [12]) and inclusive processes can be 
unified under a general treatment using a factorization approach in a generalized kinematical 
domain. The study of this domain, termed deeply virtual Compton scattering, or DVCS, is an 
area of investigation of wide theoretical interest, with experiments planned in the next few years 
at JLAB and at DESY. The key constructs of the DVCS domain are the non-forward parton 
distributions, where the term non-forward is there to indicate the asymmetry between the initial 
and final state typical of a true Compton process, in this case appearing not through unitarity, 
such as in DIS, but at amplitude level. 



Chapter 1 



The Logarithmic Expansions and Exact 
Solutions of the DGLAP Equations 
from x-Space 

1.1 New Algorithms for Precision Studies at the LHC 

The chapter is organized as follows. After defining the conventions, we bring in a simple example 
that shows how a non-singlet LO solution of the DGLAP is obtained by an x-space ansatz. Then 
we move to NLO and introduce the notion of truncated logarithmic solutions at this order, 
moving afterwards to define exact recursive solutions from x-space. In all these cases, we show 
that these solutions contain exactly the same information of those obtained in Mellin space, 
to which they turn out to be equivalent. The same analysis is then extended to NNLO. The 
approach lays at the foundation of a numerical method -based on x-space - that solves the NNLO 
DGLAP with great accuracy down to very small-x (10~^). The method, therefore, not only does 
not suffer from the usual well known inaccuracy of x-space based approaches at small-x values 
[13| . but is, a complementary way to look at the evolution of the pdf's in an extremely simple 
fashion. We conclude with some comments concerning the timely issue of defining benchmarks 
for the evolution of the pdf's, obtained by comparing solutions extracted by Mellin methods 
against those derived from our approach, in particular for those solutions which retain accuracy 
of a given order in Ug {0{as) accurate solutions), relevant for precise determination of certain 
NNLO observables at the LHC. 



1.2 Definitions and conventions 

Before we start the analysis it is convenient to define here the notations and conventions that 
we will use in the rest of the thesis. 
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We introduce the 3-loop evolution of the coupling via its /^-function 



and its three-loop expansion is 



where 



/3(«s) = -^«, -^a, -^a,+0(aj, (1.2) 



/3o = yiVc - ^T;, 

(1.3) 



are the coefficients of the beta function. In particular, P2 [HI US] and P3 [16] are in the MS 
scheme. We have set 

— I A 1 

where A'^c is the number of colors, is the number of active flavors, that is fixed by the number 
of quarks with ruq < Q. One can obtain either an exact or an accurate (truncated) solution of 
this equation. An exact solution includes higher order effects in Ug, while a truncated solution 
retains contributions only up to a given (fixed) order in a certain expansion parameter. The 
structure of the NLO exact solution of the RGE for the coupling is well known and relates 
as{^J'l) in terms of Q!s(/i|) via an implicit solution 

' . ' +ft m (4V m 1; + I , (1.5) 

where as{fj?) = as(/i^)/(47r). The truncated solution is obtained by expanding up to a given 
order in a small variable 



(1.6) 



where the n\ dependence is shifted into the factor L = h\{^\/ jjL\), and we have used a /3-function 
expanded up to NLO, involving /3o and Exact solutions of the RGE for the running coupling 
are not available (analytically) beyond NLO, while they can be obtained numerically. Truncated 
solutions instead can be obtained quite easily, for instance expanding in terms of the logarithm 



of a specific scale (A) 
47r 



^ _ Pi log La _^ 1 



La PiLj 



Pi 



(log^LA-logLA-l) +(32 



+ 



Li 



(1.7) 



where 



La = log 



MS 



and where A^^^ is calculated using the known value of as(m^) and imposing the continuity 
of Us at the thresholds identified by the quark masses. 



1.3 General Issues 

For an integro-differential equation of DGLAP type, which is defined in a perturbative fashion, 
the kernel P{x) is known perturbatively up to the first few orders in a^, approximations which 
are commonly known as LO, NLO, NNLO [T7] . 

The equation is of the form 

9/(x,Q2) _ 



51nQ2 

with 



P{x,Q')®f{x,Q'), (1.9) 



aix)®b{x)= [ —a(y)b(x/y), (1.10) 

Jo y 

and the expansion of the kernel at LO, for instance, is given by 

p(x,gY° = (^)p(°)(x). (1.11) 

In the case of QCD one equation is scalar, termed non-singlet, the other equation involves 2-by-2 
matrices, the singlet. In other cases, for instance in supersymmetric QCD, both the singlet and 
the non-singlet equations have a matrix structure [I8|. Except for the LO case, exact analytic 
solutions of the singlet equations are not known. However, various methods are available in order 
to obtain a numerical solution with a good accuracy. These methods are of two types: brute 
force approaches based in x-space and those based on the inversion of the Mellin moments. 

A brute force method involves a numerical solution of the PDE based on finite differences 
schemes. One can easily find a stability scheme in which the differential operator on the left- 
hand-side of the equation gets replaced by its finite difference expression. This method has the 
advantage that it allows to obtain the so called "exact" solution of the equation at a given order 
(LO, NLO, NNLO). The only approximation involved in this numerical solution comes from the 
perturbative expansion of the kernels. Solutions of this type are not accurate to the working 



order of the expansion of the kernels, since they retain higher order terms in Og. A short-come 
of brute force methods is the lack of an ansatz for the solution, which could instead be quite 
useful in order to understand the role of the retained perturbative logarithms. The use of Mellin 
inversion allows to extract, in the non singlet case, the exact solution quite immediately up to 
NNLO. 

The Mellin moments are defined as 

a{N) = [ a(x)x^-Mx (1.12) 
Jo 

and the basic advantage of working in moment space is to reduce the convolution product (S) into 
an ordinary product. For instance, at leading order we obtain the LO DGLAP equation (|1.9p in 
moment space 

df{N,as) _ {^)P^'Hn) 



das 



-fiN,as), (1.13) 



2p(0)(jv) ^ 

a,\ /3o .... ^ J 2PW(A^)^ fas 



which is solved by 



f{N,as)=f{N,ao)^-^J = /(iV, oo) exp | bg j^-i ) [> , (1.14) 

where we have used the notation a = a{Q'^) and ao = ck(Qo)- ^^i^ point, to contruct the 
solution in x-space, we need to perform a numerical inversion of the moments, following a contour 
in the complex plane. This method is widely used in the numerical construction of the solutions 
and various optimization of this technique have been proposed [23]. We will show below how one 
can find solutions of any desired accuracy by using a set of recursion relations without the need 
of using numerical inversion of the Mellin moments. 



1.3.1 The logarithmic ansatz in LO 

To illustrate how the logarithmic expansion works and why it can reproduce the same solutions 
obtained from moment space, it is convenient, for simplicity, to work at LO. We try, in the 
ansatz, to organize the logarithmic behaviour of the solution in terms of and its logarithmic 
powers, times some scale-invariant functions An{x), which depend only on Bjorken x. As we are 
going to see, this re-arrangement of the scale dependent terms is rather general for evolution 
equations in QCD. The number of the scale-invariant functions A„ is actually infinite, and they 
are obtained recursively from a given initial condition. 

The expansion that summarizes the logarithmic behaviour of the solution at LO is chosen of 

^It is common, however, to refer to these solutions as to the "exact" ones, though they have no better status 
than the accurate (truncated) ones. In principle, large cancelations between contributions of higher order in the 
perturbative expansion of the kernels beyond NNLO, which are not available, and the known contributions, could 
take place at higher orders and this possibility remains unaccounted for in these "exact" solution. The term 
"exact", though being a misnomer, is however wide spread in the context of perturbative applications and for this 
reason we will use it throughout our work. 



the form 



E 

n=0 



An{x) 



In 



2\\ n " 



(1.15) 



To determine An{x) for every n we introduce the simplified notation 



L = los 



a.(<5o 



2\ ' 



(1.16) 



and insert our ansatz p.l5|) into the DGLAP equation together with the LO expansion of the 
/3-function to get 



^— ' 77,1 Att ^— ' n< J/TT 



(1.17) 



n=0 n=0 

Equating term by term in powers of L we find the recursion relation 

2 



fl.18) 



At this point we need to show that these recursion relations can be solved in terms of some initial 
condition and that they reproduce the exact LO solution in moment space. This can be done 
by taking Mellin moments of the recursion relations and solving the chain of these relations in 
terms of the initial condition Aq{x). At LO the solution of (jl.lSp in moment space is simply 
given by 

AniN) = (--|pW)%(iV,a,(Q2)), (1.19) 

having imposed the initial condition Aq = g(a;, as((5o))- ^^^^ point we plug in this solution 
into (jl.lSp to obtain 



n 



2\' 



which clearly coincides with (I1.14p . after a simple expansion of the latter 



(1.20) 



oo r 

fiN,as) = fiN,ao)Y,{- 



n=0 





2P(0)(A^)' 







log" ( ^ 



(1.21) 



Notice that this non-singlet solution is an exact one. In the singlet case the same approach will 
succeed at the same order and there is no need to introduce truncated solution at this order. As 
expected, however, things will get more involved at higher orders, especially in the singlet case. 

The strategy that we follow in order to construct solutions of the DGLAP equations is all 
contained in this trivial example, and we can summarize our systematic search of logarithmic 
solutions at any order as follows: we 

1) define the logarithmic ansatz up to a certain perturbative order and we insert it into the 
DGLAP equation, appropriately expanded at that order; 



2) derive recurrency relations for the scale invariant coefficients of the expansion; 

3) take the Mellin moments of the recurrency relations and solve them in terms of the moments 
of the initial conditions; 

4) we show, finally, that the solution of the recursion relations, so obtained, is exactly the 
solution of the original evolution equation firstly given in moment space. This approach is 
sufficient to solve all the equations at any desired order of accuracy in the strong coupling, as 
we are going to show in the following sections. Ultimately, the success of the logarithmic ansatz 
lays on the fact that the solution of the DGLAP equations in QCD resums only logarithms of 
the coupling constant. 



1.4 Truncated solution at NLO. Non-singlet 

The extension of our procedure to NLO (non-singlet) is more involved, but also in this case 
proofs of consistency of the logarithmic ansatz can be formulated. However, before starting 
our technical analysis, we define the notion of "truncated solutions" of the DGLAP equations, 
expanding our preliminary discussion of the previous sections. We start with some definitions. 

A truncated solution retains only contributions up to a certain order in the expansion in the 
coupling. We could define a 1-st truncated solution, a 2-nd truncated solution and so on. The 
sequence of truncated solutions is expected to converge toward the exact solution of the DGLAP 
as the number of truncates increases. This can be done at any order in the expansion of the 
DGLAP kernels (NLO,NNLO,NNNLO,...). For instance, at NLO, we can build an exact solution 
in moment space (this is true only in the non-singlet case) but we can also build the sequence of 
truncated solutions. It is convenient to illustrate the kind of approximations which are involved 
in order to obtain these solutions and for this reason we try to detail the derivations. 

Let's consider the NLO non-singlet DGLAP equation, written directly in moment space 
and search for its exact solution, which is given by 

2p(0)(jV) 2p(°)(iV) 4p(^)(iV) 

HN,c.) = /(iV.«„)(-j ( 4./;o + .„ft j ■ 

Notice that equation p.22|) is the exact NLO equation. In particular we have preserved the 
structure of the right-hand side, that involves both the beta function and the NLO kernels and 
is given as a ratio of two polynomials in a<j 

pNLO( \ 



where 



P{x,Q^) 



"^(Q') ^ p{0)( 



x) + 



2tt 



(1.25) 



is the NLO kernel. The factorization of the LO solution from the NLO equation can be obtained 
expanding the ratio P/ (3 in a^, which allows the factorization of a l/os contribution. Equiva- 
lently, one can redefine the integral of the solution in moment space by subtraction of the LO 
part 



da 



(1.26) 



Denoting by bi = Pi/Po, the truncated differential equation can be written as 



dfiN,as 
da. 



Poas 



p(0)(Ar) + 



2lT 
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(1.27) 



which has the solution 



f{N,as) 



2P(0) 



X exp 



7r/?o 



f{N,ao). 



(1.28) 



Notice that this solution of the truncated equation, exactly as in the exact solution (ll.23p . 
contains as a factor the LO solution and therefore can be rewritten in the form 



f{N, as) = exp 
where f^^{N,as) is given by 
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ao) 



|p(°)-P«)}/^°(iV,as), 



(1.29) 
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(1.30) 



Eq. p.29|) exemplifies a typical mathematical encounter in the search of solutions of PDE's of a 
certain accuracy: if we allow a perturbative expansion of the defining equation arrested at a given 
order, the solution, however, is still affected by higher order terms in the expansion parameter 
(in our case a^). To identify the expansion which converges to (|1.29p proceeds as follows. We 
start from the 1-st truncated solution. 



Expanding (I1.29P to first order around the LO solution we obtain 



/(iV,a.) = f^^{N,as) X {l + i^^^ (|pW - P«) } , (1.31) 

which is the expression of the 1-st truncated solution, accurate at order Og. One can already see 
from (ll.3ip that the ansatz which we are looking for should involve a double expansion in two 
values of the coupling constant: Og and ao- This point will be made more clear below. For this 



reason we are naturally lead to study the logarithmic expansion 



n=0 



nl 



BJx) 



n=0 



In 



(1.32) 



which is the obvious generalization of the analogous LO expansion (ll.lSh . 



Inserting this ansatz in the NLO DGLAP equation, we derive the following recursion relations 
for An and -B„ 

An+l = -^P'^''\x)®An{x), 

Bn+1 = -Bn{x) - -P-An+l{x) - ^P'-^Hx) B„(x) - -^P^^\x) An{x), 



(1.33) 



together with the initial condition 



f{x,Ql) = Ao + asiQl)Bo. 



(1.34) 



At this point we need to prove that the recursion relations p.33|) reproduce in moment space 
p.3ip . To do so we rewrite the recursion relations in Mellin-space 



An+i{N) = --p('\N)AniN), 

PO 

Bn+i{N) = -Bn{N) - Jl-An+i{N) - lp(0){N)Bn{N) - -^p(^) {N)An{N), 

4vr/?o PO T^Po 

(1.35) 

and search for their solution over n. After solving these relations with respect to and Bq, it 
is simple to realize that our ansatz p.32|) exactly reproduces the truncated solution (ll.3ip only 
if the condition i?o = is satisfied. In fact, denoting by 



Rq 

Ri 



Po 



P(0)(A^), 



.p(l) 



(1.36) 



the recursive coefficients, we can rewrite the recursion relations as 



An+l — RoAn, 

Bn+1 = (Rq - l)Bn + RlA, 



(1.37) 



Then, observing that 



Bi = {Ro - l)Bo + RiAo, 

B2 = {Ro - if Bo + RiAo{2Ro - 1), 

Bs = {Ro - if Bo + RiAo [{2Ro - l)(i?o - 1) + ^0] 



(1.38) 



we identify the structure of the n^/i iterate in close form 

Brr = {Ro - If Bo + RiAo [R'i - {Ro - If] • 

Substituting the expressions for A„ and Bn so obtained in terms of ^0 and Bo in the initial 
ansatz, and summing the logarithms (a procedure that we call "exponentiation ") we obtain 

Va.^^V^L" = y^a, {{R^-irBo + RiAo[R^o-{R^-m} 

z — / ri.\ ^ — ^ n 



n\ ^ — ' n 
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(1.39) 



expression that can be rewritten as 



f{N,as)=Ao(—Y\asBo(—Y" \ a^RiAo ( —Y' - asRiAo ( —Y' \ (1.40) 
This expression after a simple rearrangement becomes 

f{N, as) = (—) ° [1 + (a. - ao) Ri] f{N, Q^) - (^) ° (a, - ao) {RicxoBo)- (1-41) 

This solution in moment space exactly coincides with the truncated solution (|1.3ip if we impose 
the condition Bq = 0. It is clear that the solution gets organized in the form of a double expansion 
in the two variables a^ and ao- While Og appears explicitely in the ansatz p.32|) . ao appears 
only after the logarithmic summation and the factorization of the leading order solution. An 
obvious question to ask is how should we modify our ansatz if we want to reproduce the exact 
solution of the truncated DGLAP equation in moment space, given by eq. p.28|l . The answer 
comes from a simple extension of our recursive method. 



1.4.1 Higher order truncated solutions 



We start by expanding the solution of the truncated equation (|1.29l) . whose exponential factor 
is approximated by its double expansion in as and ao to second order, thereby identifying the 
approximate solution 



f{N,as 



aoj 
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1 - Ri{ao - as) + ^ii?(ao - a^)^ + Ri{al - a^)^ 
z ovr 



fiN,ao). 



(1.42) 



To generate this solution with the recursive method it is sufficient to introduce the higher order 
(2nd order) ansatz 



L 



f{x, "s) = ^ [An{x) + asBn{x) + a^C„(x)] , 



(1.43) 



n=0 



where we have included some new coefficients C„(x) that will take care of the higher order terms 
we aim to include. Inserted into the NLO DGLAP equation, this ansatz generates an appropriate 
chain of recursion relations 

An+l{x) = -|-P(°)(x) 0^„(x), 

m 

Bn+l{x) = -Bn{x) - ^P^^\x) ® Bn{x) - -^^„+i(x) - -^P^^\x) ® A„(x), 



/3o 



(47r)- 



Cn+l{x) = -2Cn{x) - ^P^^\x) Cnix) - -^Bn+l{x) - -^Bn{x) - -^P^^\x) ® S„(x), 

PO (4vrj (47rj TTfjQ 



(1.44) 



that we solve by going to Mellin space and obtain 



Bn = Ri [R^ - {Ro - 1)"] ^0, 



Cn. 



{Ro - 2)" - [Ro - ir 



RlAo + 



-RiRqAq + —Ribi {Ro — 2)"'j4o — —RiRobiAo, 
z ovr ovr 



(1.45) 



where the initial conditions are f{N,ao) = Ao and -60 = ^^0 = 0. It is a trivial exercise to 
show that the solution of the recursion relation, inserted into (|1.43|) . coincides with (|1.42|) . after 
exponentiation. 

Capturing more and more logs of the truncated logarithmic equation at this point is as easy 
as never before. We can consider, for instance, a higher order ansatz accurate to ©(af) for the 



NLO non-singlet solution 



L 



f{x, "s) = X] [Anix) + asBnix) + a^C„(x) + a^Dnix)] , 



(1.46) 



n=0 



that generates four independent recursion relations. The relations for ^n+i, Bn+i, Cn+i are, for 
this extension, the same as in the previous case and are listed in (|1.45p . Hence, we are left with 
an additional relation for the Dn+i coefficient which reads 



Dn+l{x) = -3Dn{x) - Dn{x) - -^Cn+lix) - -^Cn{x) - -^P^^^ Cn{x). 



(1.47) 



These are solved in Mellin space with respect to Aq, Bq,Co, Dq (with the condition Bq = Co 
Do = 0). We obtain 



An — ^0^0 5 
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Exponentiating we have 

f{x, as) = {1 + as \ I 
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(1.49) 



which is the solution of the truncated equation computed with an 0{al) accuracy. 



1.5 Non-singlet truncated solutions at NNLO 



The generalization of the method that takes to the truncated solutions at NNLO is more involved, 
but to show the equivalence of these solutions to those in Mellin space one proceeds as for the 
lower orders. As we have already pointed out, one has first to expand the ratio P/ (3 at a certain 
order in a^, then solve the equation in moment space - solution that will bring in automatically 
higher powers of Ug - and then reconstruct this solution via iterates. 

At NNLO the kernels are given by 
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pW(x) + 
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(1.50) 



and the equation in Mellin-space is given by 

df{N,as) pNNLO^^^ 



da. 



NNLO 



f{N,as 



We search for solutions of this equation of a given accuracy in and for this purpose we 
truncate the evolution integral of the ratio P/P to O(a^). This is the first order at which the 
P^"^^ component of the kernels appear. As we are going to see, this will generate the first truncate 
for the NNLO case. Therefore, while the first truncate at NLO is of 0(as), the first truncate at 
NNLO is of O(a^). We obtain 
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f PNNLoix,a) _ PLoja) 



-Riao - ^R2al + RiOs + ^R2a1. (1-52) 



At this retained accuracy of the evolution integral, the exact solution of the corresponding 
(truncated) DGLAP equation can be found, in moment space, as in p.72|) 
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(1.53) 



At this order this solution coincides with the exact NNLO solution of the DGLAP equation, 
obtained from an exact evaluation of the integral p.52|) . followed by a double epansion in the 
couplings. Therefore, similarly to (|1.42p . the solution is organized effectively as a double ex- 
pansion in as and ao- This approach remains valid also in the singlet case, when the equations 



assume a matrix form. As we have already pointed out above, all the known solutions of the 
singlet equations in moment space are obtained after a truncation of the corresponding PDE, 
having retained a given accuracy of the ratio P/ (3. For this reason, and to compare with the 
previous literature, it is convenient to rewrite (|1.53p in a form that parallels the analogous singlet 
result [19]. It is not difficult to perform the match of our result with that previous one, which 
takes the form p]. [2]. [20] 
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which becomes, after some manipulations 
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where the functions Ui{N) are defined as 
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We intend to show rigorously that this solution is generated by a simple logarithmic ansatz 
arrested at a specific order. For this purpose we simplify (|1.55p obtaining 
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and the ansatz that captures its logarithmic behaviour can be easily found and is given by 
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Setting the initial conditions as 
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fix, Ql) = Aq{x) + a^Boix) + alCo{x), 



fl.58) 
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and introducing the 3-loop expansion of the /3-function, we derive the following recursion relations 
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We need to show that the solution of the NNLO recursion relations reproduces (11.55(1 in 
Mellin-space, once we have chosen appropriate initial conditions for Ao{N), Bo{N) and Co{N) 



i 



At NLO we have already seen that Bq{N) has to vanish for any A^, i.e. Bq{x) = 0, and we 
try to impose the same condition on Co{N). In this case we obtain the recursion relation for the 
moments 



C„+i(iV) = -2Cn{N) - -^Bn{N) - J^Bn+i{N) 
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(1.61) 



^It can be shown that the infinite set of recursion relations have internal symmetries and different choices of 
initial conditions can bring to the same solution. The choice that we make in our analysis is the simplest one. 



which combined with the relations for An{N) and Bn{N) give 

C„ = I -Rl {Ro -ir- \r2R^ + \ [(i?o - 2)" Rl + {Ro - 2)" R2 + RIR^] | Ao, 
En = [R'S - {Ro - IT] RiAo, 

An = R'SAo, (1.62) 

where the dependence in the coefficients Ri has been suppressed for simplicity. The solution 
is determined exactly as in (I1.37P and (|1.39ll and can be easily brought to the form 

1 /a \^"(^) 

f{N,Q^) = -i — ) [2 - 2Ri{as - ao) + Rl{as - aof + R2{al - a'i)] Ao, (1.63) 

2 \aoJ 

which is the result quoted in eq. p. 571) with ^0 = f{^, ckq)- 

We have therefore identified the correct logarithmic expansion at NNLO that solves in x-space 
the DGLAP equation with an accuracy of order a^. 



1.6 Generalizations to all orders : exact solutions of truncated 
equations built recursively 

We have seen that the order of approximation in Us of the truncated solutions is in direct 
correspondence with the order of the approximation used in the computation of the integral 
on the right-hand-side of the evolution equation (|1.26p . This issues is particularly important 
in the singlet case, as we are going to investigate next, since any singlet solution involves a 
truncation. We have also seen that all the known solutions obtained in moment space can be 
easily reobtained from a logarithmic ansatz and therefore there is complete equivalence between 
the two approaches. We will also seen that the structure of the ansatz is insensitive to whether 
the equations that we intend to solve are of matrix forms or are scalar equations, since the ansatz 
and the recursion relations are linear in the unknown matrix coefficients (for the singlet) and 
generate in both cases the same recursion relations. 

We pause here and try to describe the patterns that we have investigated in some generality. 
We work at a generic order N'^LO, with m=l denoting the NLO, m = 2 the NNLO and so 
on. We have already seen that one can factorize the LO solution, having defined the evolution 
integral 

lN'^Lo{as,ao) = dal— — ^ , (1-64) 

Jao V PN"^Lo{a) PLo{a)J 

and the exact solution can be formally written as 

f{N, as) = f^°{N, as) x e^Nr^Loi<^s,»o) _ (165) 



A Taylor expansion of the integrand in the (ll.64p around = at order k = {m — 1) gives, in 



moment space, 



V PN"^LO[0!s) PL0{0ls) J 

+ ■■■ + . . . , iV)a:, (1.66) 



which at NNLO becomes 

/ PNNLo{N a ) _ PloM\ ^ R^^pm^piD^N) + R^{P^'\P^'\P^^\N)as. (1.67) 
V PNNLoias) PLo{as) J 

Thus, integrating between ag and oq we, obtain the following expression for I^mx^oiois, 0:0) at 

OK) 

■^^^™LoK, ao) = --Riao - -^^200 -RkCHq + RiOs + ^i?2as ^ --R^as ^ (1-68) 

z /°c .z 

where the {P^^\ P^^\ . . . , P^"^\ N) dependence in the Rk coefficients has been omitted. To 
summarize: in order to solve the K-truncated version of the N'^LO DGLAP equation 

df{N, as) -fVLo(^Vi£) ^ 

^ = —3 7 — ^f{N,as), (1.69) 

das PN'^Lo(as) 

which is obtained by a Taylor expansion - around = - of the ratio Pn'^^lo{N, as)//?N™Lo(as), 
we need to solve the equation 

^-^^i^'"'^ = — [Ro + asRi + alR2 + ■■■ + f{N, as), (1.70) 

OUs Us 

where the coefficients R^ have a dependence on P^^^ and P^^^ in the NLO case, and on P^^\ 
and in the NNLO case. Eq. (I1.70p admits an exact solution of the form 

f{N, as) = (^^^ exp {as - ao) + ^i?2 (a^ - ag) + • • • + ^R^ (a^ - aS)| /(iV, oq), 

(1.71) 

having factorized the LO solution. 

At this stage, the Taylor expansion of the exponential around {as,ao) = (0,0) generates an 
expanded solution of the form 

fN^"Lo{N,as) « fLo{N,ao)e^N^Lo 

= /Lo(iV,ao)(l + /£LO + |(4™Lo)' + •••)' (1-72) 



which can be also written as 

fLo{N, ao) (l + /&io + ^ {InIlo) ' + •••) = 
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(1.73) 



where the coefficients Cij are defined in moment space. The K-th truncated solution of the 
equation above 

/Ar™Lo(iV,a.)|o(^j) = /^°(iV,a,) Y.'^'-\\alci^j | , (1.74) 

is therefore accurate to O(a^), and clearly does not retain all the powers of the coupling constant 
which are, instead, part of p.72|) . However, as far as we are interested in an accurate solution 
of order k, we can reobtain exactly the same expression from x-space using the ansatz 



oo 

/jv™Lo(3;,a.)|o(„j') = {^li^)+<^sAl{x)+alAl{x) + --- + a'; Ai{x)^ 

n=0 

(1.75) 

which can be correctly defined to be a truncated solution of order k' of the («:)-truncated equation. 
Since we have truncated the evolution integral at order k, this is also the maximum order at 
which the truncated logarithmic expansion (|1.75p coincides with the exact solution of the full 
equation. This corresponds to the choice n' = k. Notice, however, that the number of coefficients 
that one introduces in the ansatz is unrelated to k and can be larger than this specific value. 
This implies that we obtain an improved accuracy as we let k' in (|1.75p grow. 

If we choose the accuracy of the evolution integral to be k, while sending the index k' in the 
logarithmic expansion of p. 751) to infinity, then the ansatz that accompanies this choice becomes 

oo / oo 
n=0 \l=0 

and converges to the exact solution of the (order k) truncated equation (|1.72l) . Obviously, this 
exact solution starts to difi'er from the exact solution of the exact DGLAP equation at 0(0^"*"^). 
Also in this before, one should notice that the double expansion in as and oq of the 

exact solution of the (K)-truncated equation can be reobtained after using our exponentiation, 
and not before. We remark, if not obvious, that the recursion relations, in this case, need to be 
solved at the chosen order k', as widely shown in the examples discussed before. 

We remark that, as done for the LO, we could also factorize the NLO solution and determine 



1 /^iM^ 
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the N'^LO solution using the integral 
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(1.77) 



and then restart the previous procedure. Obviously, the two approaches imply a resummation 
of the logarithmic behaviour of the pdf 's in the two cases. 

It is convenient to summarize what we have achieved up to now. A truncation of the evolution 
integral introduces an approximation in the search for solutions, which is controlled by the 
accuracy (k) in the expansion of the same integral. The exact solution of the corresponding 
truncated equation, as we have seen from the previous examples, involves all the powers of 
Us and ao stnd, obviously, a further expansion around the point = ao = is needed in 
order to identify a set of truncated solutions which can be reobtained by a logarithmic ansatz. 
This is possible because of the property of analiticity of the solution. Therefore two types of 
truncations are involved in the approximation of the solutions: 1) truncation of the equation and 
2) truncation {k') of the corresponding solution. In the non-singlet case, which is particularly 
simple, one can therefore identify a wide choice of solutions (by varying k and k') that retain 
higher order effects in quite different fashions. Previous studies of the evolution using an ansatz 
a la Rossi-Storrow [6j[7], borrowed from the pdf's of the photon, were therefore quite limited in 
accuracy. Our generalized procedure is the logical step forward in order to equate the accuracy 
of solutions obtained in moment space to those in x-space, without having to rely on purely 
numerical brute force methods§| 



1.6.1 Recursion relations beyond NNLO and for all k's 

In the actual numerical implementations, if we intend to use a generic truncate of the non- 
singlet equation (the result is actually true also for the singlet), it is convenient to work with 
implementations of the recursion relations that are valid at any order. In fact it is not practical 
to rederive them at each new order k' of approximation. Next we are going to show how to do 
it, identifying generic relations that are easier to implement numerically. 

The expression to all orders of the DGLAP kernels is given by 
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(1.78) 



^In [2] a flag variable called IMODEV allows to switch among the exact solution (IM0DEV=1), the exact 
solution of the 0(0?) truncated equation (IM0DEV=2). A third option (IMODEV=0) involves at NLO and at 
NNLO 0{as) and 0(a?), respectively, truncated ansatz that, in our cases, are reconstructed logarithmically. 



and the equation in Mellin-space in Melin space is given by 
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whose exact solution can be formally written as 

V"o/ 
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(1.80) 



where is obtained from the evolution integral and whose specific form is not relevant at this 
point. 

We will use the notation (P, (3) to indicate the sequence of components of the kernels and 
the coefhcients of the /3-function {P^^\P^^\P^'^\ ...). 

Then, the Taylor expansion around = ao of the solution is formally given by 



(1.81) 



?i=0 



for an appropriate ^„{dT,d'^J^,...,d'^J^). is a function that depends over all the partial 
derivative obtained by the Taylor expansion. Since it is calculated for ag = uq, it has a parametric 
dependence only on ao, and we can perform a further expansion around the value ag = 
obtaining 



ml dan n! 



m=0 



.n=0 



(1.82) 



ao=0 



This expression can always be arranged and simplified as follows 
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n=0 



,n-l 



^ci,„(0, Pj)+i -n^ci,„(o, P, 0) + ^d^n{ao,P, 



ao=0. 



1 (n — l)na 



n-2 



a 



n-l 



2 
1 



n! 



-$„(0, P, P) - n^9$„(ao, P, /?) 
n! 



ao=0 



ay'^n{ao,P,(3) 



"0 + (• • • )ao + • • • 



ao=0 . 



l + a,Cf)(P,/3) + ... + a^^W(p,^) 
+ao (e^ /3) + (P, ^) + a^e^') (P, /3) + • • • + a^^^^) (P, /3) 



(e^^ (P, /3) + a,C^) (P, /3) + a^^f ) (pJ) + ... + a^^M (p, /3) 

oo oo 



71=0 771=0 



(1.83) 



where we are formally absorbing all the dependence on both the kernels P and the /3-function /3, 
coming from the functions d"^^n calculated at the point ao = 0, in the coefficients ^i™^ Finally, 
we can reorganize the solution to all orders as 



^) ° f{N,ao)Y,Y.«^^rrHPJ)- (1-84) 

"° ^ n=0 m=0 

With the help of the general notation 

P^,^ = ^pm^pii)^^ 

Pnnlo = {P^'\p^'\P^'^), 
Pnlo = {Po,Pi), 

PNNLO = iPo,Pl,f32), (1.85) 

we can try to indentify, by a formal reasoning, the exact solution up to a fixed - but generic - 
perturbative order of the expansion of the kernels. 

To obtain the NLO/NNLO exact solution it is sufficient to take as null the components 
(p{2)^p{3)^_ g^^^ (/32,/33,...) for the NLO and (P(^\...) and (/^s,-) for the NNLO case. Since 
eq. (ll.84p contains all the powers of aoUs up to a^a"^ (i.e. it is a polynomial expression of order 



a""*"™), if we aim at an accuracy of order a^, we have to arrange the exact expanded solution as 



f{N,as 



«£_ Y 



f{N, ao) Yl c^K-'d'HP, P) + 0{a^) 

n=0 j=0 



(1.86) 



with 0{ag) indicating all the higher-order terms containing powers of the type ao^s +■ ■ • + aQa'^. 
Hence, the index k represents the order at which we truncate the solution. 

As a natural generalization of the cases discussed in the previous sections, we introduce the 
higher-order ansatz (K-truncated solution) 



fiN,as)=Y 



n=0 



.m=0 



n! 



log^ 



(1.87) 



that reproduces the exact solution (|1.86p expanded at order k. 

In fact, inserting this last ansatz in p.79p we generate a generic chain of recursion relations 
of the form 



00+i(7V)=FO(00(pW,/3o)), 
Ol^.iN) = F\OlO^^+„OlPj), 

o«+i(iV) = F«(oo,...,o;;,oo+i,...,o«;J,p,/3). (1.88) 

A deeper look at the explicit structure of (ll.88p . for the NLO non-singlet case, reveals the 
following structures for the generic iterates 



0°+i(iV) 



/3o 
2 

To 

-4o-(iv)-.o^(^) 



p(0)o- (N)-— P«(Ar)Or^(iV) 



(k - 1 



'Aori(iv), 



(1.89) 



at NLO, while for the NNLO case we obtain 

2 r 



0°+i(iV) 



2 r 

To 



p(°)(Af)0°(A^) 
pW(iV)oi(iV) 



1 

7r/?o 



2 r 



o^+i(iv) = -- p(o)(iv)o;:(iv) 

PO L 



27r2/3o 



1 



P«(iV)O0(iV) 



-«0::(iV) - (k - l)^o:~\N) - - 2) 
4vrpo 



/32 



167r2/3o 



Hence, one is able to determine the structure of the n-th recursion relation when the k = and 
K = 1 cases are known. This property is very useful from the computational point of view. 

Passing to the resolution of the recursion relations in moment space, we get the formal 
expansion of each 0'^{N) in terms of the initial condition Oq{N), which reads 



fiN, ao) = Og(iV) + ao + + ... + aS . 



(1.91) 



Here we have set to zero all the higher order terms, as a natural generalization of Bq = Cq = 0..., 
according to what has been discussed above. 

These relations can be solved as we have shown in previous examples, and the generic struc- 
ture of their solution can be identified. If we define Rq = -^P^^\ the expressions of all the 
0!^{N) in terms of f{N,ao) = fo become 

OliN)=G\RS,iRo-ir,Pj)fo, 

Ol{N) = G\R^, (Ro - If, {Ro - 2)", P,/3)/o, 



0;:(7V) = G^{R^, {Ro - ir, {Ro - 2)", {Ro - k)", P, /?)/o 



(1.92) 



and in particular, by an explicit calculation of 0!^{N), one can work out the structure of these 
special functions C^. For instance, we get for m = k the expression 



G'^(Po", (i?o - 1)", {Ro - 2)", {Ro - k)", P, /3) = ^^(Po - jT^i^HP, P) , (1-93) 

j=0 



*Tb.e relations l|1.89p and l|1.90p hold also in the NLO/NNLO singlet case and can be generalized to any 
perturbative order in the expansion of the kernels. 



for suitable coefficients . Substituting the 0^{N) functions with m = 0, . . . k in the higher- 
order ansatz (|1.87|) and performing our exponentiation, we get an expression of the form 




(1.94) 

which can be written as follows by the use of eq. (I1.93P 




(1.95) 



This is the exact solution expanded up to 0{ag) in accuracy. 



1.7 The Search for the exact non-singlet NLO solution 

We have shown in the previous sections how to construct exact solutions of truncated equations 
using logarithmic expansions. We have also shown the equivalence of these approaches with 
the Mellin method, since the recursion relations for the unknown coefficient functions of the 
expansions can be solved to all orders and so reproduce the solution in Mellin space of the 
truncated equation. The question that we want to address in this section is whether we can 
search for exact solutions of the exact (untruncated) equations as well. These solutions are known 
exactly in the non singlet case up to NLO. It is not difficult also to obtain the exact NNLO solution 
in Mellin space, and we will reconstruct the same solutions using modified recursion relations. 
The expansions that we will be using at NLO are logarithmic and solve the untruncated equation. 
The NNLO case, instead, will be treated in a following section, where, again, we will use recursion 
relations to build the exact solution but with a non-logarithmic ansatz 

The exact NLO non-singlet solution has been given in eq. (I1.23|) . The identification of an 
expansion that allows to reconstruct in moment space eq. (|1.23p follows quite naturally once the 
typical properties of the convolution product (X" are identified. For this purpose we define the 
serie of convolution products 

' nl 

n=0 



^In PEGASUS [2] the NNLO non-singlet solution is built by truncation in Mellin space of the evolution 
equation, while the NLO solution is implemented as an exact solution. 



that acts on a given initial function as 

gF^p^(x)®^(^) = (^5(1 -x)®+FaPa(x)0+^f1Pa«> Pa + <A(a:) 

= <Pix) + Fa (Pa ® </>) + ^^1 (Pa ® Pa®) H^) + • • • • (1-97) 

The functions Fa and Fb are parametrically dependent on any other variable except the 
variable x. The proof of the associativity, distributivity and commutativity of the product is 
easily obtained after mapping these products in Mellin space. For instance, for generic functions 
o(x) b{x) and c{x) for which the (gi product is a regular function one has 

M[{a^b)^c]{N) = M[a^{b^c)]{N) 

= a{N)B{N)C{N), (1.98) 

where Ai denotes the Mellin transform and N is the moment variable. Also one obtains 

e^APAW® e^BPflW® ^(^) = ^{FAPA{x)+FsPB{xm ^(^)^ (1.99) 

and 



M 



(AT) = e^-^'^(^)(/)(7V), (1.100) 



with M x-independent, since both left-hand-side and right-hand-side of (|1.99|) can be mapped 
to the same function in Mellin space. Notice that the role of the identity in (g-space is taken by 
the function 5(1 — x). We will also use the notation 

(f2A'nix)FU <p{x) = iAo{x)<E)+FAAi{x)'S) + ...)<J)ix), (1.101) 

\n=0 / (g, 

where the A'n(x) and the An{x) capture the operatorial and the functional expansion - respec- 
tively - and are trivially related 

<(x) (t>{x) = An{x). (1.102) 
To identify the x-space ansatz we rewrite (I1.23P as 

/(JV,c) = nN,ao)e'm'-e'm" 



E^-" E^M"), ,U03, 



„ n! / \ ^ m! 

Kn=0 / \m=0 



where we have introduced the notations 



r 1 

L = log—, 



M = log 

a{N) = 



AirPo + aoPi ' 
2PW(iV) 
Wo ' 



m) = - ■ (1-104) 

Pi 

Our ansatz for the exact solution in x-space is chosen of the form 



n\ I \ ^ — ' ml 

\n=0 / (5?i \m=0 



y y L-M^^-^^^^^^k^ ® fix, Q, 

n!(.-n)! 



2n 



oo s 



s=0 n=0 ^ ^ 

where in the first step we have turned the product of two series into a single series of a combined 
exponent s = n + m, and in the last step we have introduced the functions 

C'^{x)=AUx)0B',_^{x)0f{x,Ql), {n<s). (1.106) 

Setting Q = Qq in p.lQSp we get the initial condition ^o(^) ~ -^o(^) ~ "^(1 ~ ^) equivalently, 

f{x,Ql) = C^{x). (1.107) 

Inserting the ansatz p.lOSp into the NLO DGLAP equation, with the expressions of the kernel 
and beta function included at the corresponding order, we obtain the identity 

s=On=0 ^ ^ J \ J 



y y I f p(°) + JiiMTi, (1.108 



=0 ra=0 

Equating term by term the coefficients of a and a^, we find from this identity the new exact 
recursion relations 

C'ntl = -^P^'^^Ci, (1.109) 
= -C7^+1-^P«®C^, (1.110) 



or, equivalently, 



c„ 



Po 

^"+1 Pi 



a 



s-l 



(1.111) 

(1.112) 



Notice that although the coefficients are convolution products of two functions, the recursion 
relations do not let these products appear explicitely. These relations just written down allow to 
compute all the coefficients (n < s) up to a chosen s starting from Cg, which is given by the 
initial conditions. In particular eq. (11.1111) allows us to move along the diagonal arrow according 
to the diagram reproduced in Table 11.11 eq. (11.1121) instead allows us to compute a coefficient in 
the table once we know the coefficients at its right and the coefficient above it (horizontal and 
vertical arrows). To compute there is a certain freedom, as illustrated in the diagram. For 
instance, to determine C| we can only use (ll.lllj) . and for the coefficients Cq we can only use 
(|1.112p . For all the other coefficients one can prove that using (ll.lllj) or (I1.112|) brings to the 
same determination of the coefficients, and in our numerical studies we have chosen to implement 
p.llip . being this relation less time consuming since it involves P^^^ instead of P^^\ 

The recursion relations defining the iterated solution can be solved as follows. From the first 
relation (jl.llip . keeping the s-index fixed, we have 



Po 



Po 



.p(0) 



-n-l 



(1.113) 



then, since the second relation p.ll2p also holds for n = 0, we can write (using eq. (I1.113|) ) 

4 „n^ „o_i 
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Finally, inserting the above 
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in (|1 
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Ap{o) _ ±p(i) 



Pi 



(1.114) 



(1.115) 



which is the solution we have been searching for. The last step in the proof consists in taking 
the Mellin transform of this operatorial solution and summing the corresponding series 



s=0 n=0 
oo s 

= EE 

s=0 n=0 

that after summation gives 



nl{s — n)l 
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2.p(0) _ ±p(l) 



Po 



Pi 



CO(iV), (1.116) 



fiN,as) = e ^0 



exp 



PO Pi 



(1.117) 



which is exactly the expression in eq. p.l03p . Hence, it is obvious that the exact solution of the 
DGLAP equation (ll.22p can be written in x-space as 



f{x,asiQ' 



fl.118) 



therefore proving that the ansatz (|1. 105(1 reproduces the exact solution of the NLO DGLAP 
equation from x-space. 

A second version of the same ansatz for the NLO exact solution can be built using a factor- 
ization of the NLO DGLAP equation. This strategy is analogous to the method of factorization 
for ordinary PDE's. For this purpose we define a modified LO DGLAP equation, involving (5^^'^ 



dfLo{x,as) 



da. 



a. 



27r/? 



NLO 



P^^\x)<E)fLoix,as), 



(1.119) 



whose solution is given by 



fLo{x,a) = e 



M 



fLo{x,a) = e V y f{x,ao) 



2P(0) 
-2P(0) 



fLo{x,a) 



(1.120) 



where we have introduced the ordinary LO solution /lo; expressed in terms of a typical initial 
condition /(x,ao); and the NLO recursion relations can be obtained from the expansion 



fNLo{x,a) = (f2^4^M^ ) fLo(.x,a) 



\n=0 



ni 



(1.121) 
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Table 1.1: Schematic representation of the procedure followed to compute each coefficient C^. 
Inserting this relation into (|1.22p we obtain the recursion relations 



Bq{x) = (5(1 -x), (1.122) 



which is solved in moment space by 



Bn{N) = (^-±p(^)y Bo{N) 

Bo{N) = I. (1.123) 



The solution eq. (|1.121|) can be re-expressed in the form 



fNLo{x,a) = V ^0 ; e V "0 ; 



e 



^ fix,ao) (1.124) 



which agrees with p.l03p once a{N) and b{N) have been defined as in (|1.104|) . 

We have therefore proved that the exact NLO solution of the DGLAP equation can be 
described by an exact ansatz. Since the ansatz is built by inspection, it is obvious that one 
needs to know the solution in moment space in order to reconstruct the coefficients. Though 
the recursive scheme used to construct the solution in x-space is more complex compared to the 
recursion relations for the truncated solution, its numerical implementation is still very stable 
and very precise, reaching the same level of accuracy of the traditional methods based on the 
inversion of the Mellin moments. 



1.8 Finding the exact non-singlet NNLO solution 

To identify the NNLO exact solution we proceed similarly to the NLO case and start from the 
DGLAP equation in moment space at the corresponding perturbative order 



a/(iv..,)_(ft)P'°'(ivH(a-pWW + a)^K-'(A^) 

After a separation of variables, all the new logarithmic/non logarithmic and dependences come 
from the integral 

fas(Q^) pNNLOl \ 

/ . 3NNLo\^^ ^ (1-126) 

and the solution of p.l25l) is 

\n=0 ■ / \m=0 ■ / \p=0 ^' J 

where we have defined 

C = log—, (1.128) 

_ Idn^Po + 47ra/?i + a^/?2 . . 

167r2/3o + 47rao/3i + «g/32 ' ^ ^ ^ 

Q = — = arctan — ; — ; — — — — — , (1.130) 

V4/3o/32 - Pt 27r(87r/3o + (a + ao)/3i) + aao/?2 

a{N) = (1.131) 
KN) = (1.132) 

PO P2 

c{N) = ^p(0)(Ar)_8p(i)(Ar) + Mlp(2)(^). (1.133) 

Po P2 

Notice that for nj = 6 the solution has a branch point since 4/3o/?2 — /3i < 0. If we increase as 
we step up in the factorization scale then, for nj = 6, Q is replaced by its analytic continuation 



- , ' ■ „ (1.134) 

^yp'f - 4/3o/32 27r(87r/?o + (a + ao)/3i) +aao/32 

Eq. (|1.127|1 incorporates all the nontrivial dependence on the coupling constant as (now deter- 
mined at 3-loop level) into C,M. and Q. 

As a side remark we emphasize that it is also possible to obtain various NNLO exact recursion 
relations using the formalism of the convolution series introduced above. For this purpose it is 
convenient to define suitable operatorial expressions, for instance 

£2 = e^"o ^'^^ (1.135) 



which are manipulated under the prescription that the integral in a is evaluated before that 
any convolution product acts on the initial conditions. The re-arrangement of these operatorial 
expressions is therefore quite simple and one can use simple identities such as 

,a. P«')(x.o) ,,A„p,0,_A„,0,^i_,,p,0| 



= I/"l2^jF^, = V^« 



as 



JnNLO = 1^^ da ^NNLO^a) ^ 



Q ( 2|ipW ^ _8p(l)^^ _ 2 £p(0) ^ + 1 _^p(0)^^ 



(1.136) 



to build the NNLO exact solution using a suitable recursive algorithm. For instace, using (I1.135|) 
one can build an intermediate solution of the equation 

dfNLoix,as) 

= pNNLo f^Loix, as) (1.137) 

given by 

Inlo = £if{x,ao), (1.138) 

and then constructs with a second recursion the exact solution 

f{x,as) = S2fNLO- (1.139) 

A straightforward approach, however, remains the one described in the previous section, that we 
are going now to extend to NNLO. In this case, in the choice of the recursion relations, one is 
bound to equate 3 independent logarithmic powers of £, Ai and Q that appear in the symmetric 
ansazt 



/(.,g^) = (f:^^') ff;^-^") [f^^^i 

\n=0 • / \m=0 ■ / ^ \p=0 J ^ 

t^oht'o n\{t-n)\{s-ty. 



oo s t 

n\{t - ny.{s - t)\ 



s=Q t=0 n=0 



and where 



DlJx) = <(x) B[_^{x) ® ® fix, Ql). 



il.Ul) 



The ansatz is clearly identified quite simply by inspection, once the structure of the solution in 
moment space p.l27p is known explicitely. In (|1.140|) we have at a first step re-arranged the 
product of the three series into a single series with a given total exponent s = n + m + p, and 
we have introduced an index t = n + m. The triple-indexed function Z)|„(x) can be defined also 
as an ordinary product 



Dl^{x) = An{x) Bt^nix) ^ Cs-t{x), {n < t < s), 
where we have absorbed the operator into the definition of A, B and C, 



A{x)B{x)C{x) = A'{x) ® 



B'{x)(^(c'{x)®f{x,Ql))) 



Setting Q = Qq \n (11.1401) we get the initial condition 

f{x,Ql) = DUx). 



(1.142) 



(1.143) 



(1.144) 



Inserting the ansatz (|1.140p into the 3-loop DGLAP equation together with the beta function 
determined at the same order and equating the coefficients of a, and a^, we find the recursion 
relations satisfied by the unknown coefficients Z)f„(x) 
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In the computation of a given coefficient -D|„, if more than one recursion relation is allowed to 
determine that specific coefficient, we will choose to implement the less time consuming path, 
i.e. in the order (11.148(1 . (|1.150p and (ll.MQh . At a fixed integer s we proceed as follows: we 



1. compute all the coefficients with n / using p.l48p : 

2. compute the coefficient D^q using (I1.149|) : 

3. compute the coefficient D^q with t ^ s using (11.150(1 . in decreasing order in t. 
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Table 1.2: Schematic representation of the procedure followed to compute each coefficient 
for s = 4. The underlined coefficients are computed via eq. p.l48p . 



This computational strategy is exemplified in the diagram in Table 11.21 for s = 4, where the 
various paths are highlighted. 

Following a procedure similar to the one used for the NLO case, we can solve the recursion 
relations for the NNLO ansatz with the initial conditions Do q(x). Solving the relations (|1.1481 
ll.lSOj) . we obtain the chain conditions 
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Then, from the second relation we get the additional ones 
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From the last relation we also obtain the relations 
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(1.153) 



which solve the recursion relations in x-space in terms of the initial condition Dq q. Finally, the 
explicit expression of the Df^ coefficient will be given by 
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The solution of the NNLO DGLAP equation reproduced by (|1.140p in Melhn space will then be 
written in the form 



UN. a.) 



EEE 

s=0 t=0 n=0 

oo s t 

EEE 

s=0 t=0 n=0 



DtjN) 



£nj^t-nQS-t 



n\{t - n)\{s - t)\ 

n\{t - n)\{s ~ t)\ 
2f3i 



.Ap{0) 



.8P(1) (iV) + ^P(O) ^ gPi p(2) 



/5o 

s— 2n 



/32 



t~n 



(1.155) 



which is equivalent to 

f{N,as) = 

exp 



^P(0)|10B 



exp 



p(0)(iV) 4P(2)(]v)' 



/3o 



/?2 



log 



167r^/3o + 47ra^/?i + ag/?2 
167r2/3o + inaoPi + ao/32 



PO /52 



: arctan ■ 



DUN), (1.156) 



^ V4/3o/32 - /5i 27r(87r/3o + (as + ao)/?i) + asao/32^ 

and it reproduces the result in (|1.127p . In x-space the above solution can be simply written as 
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We have shown how to obtain exact NNLO solutions of the non-singlet equations using recursion 
relations. It is clear that the solution shown above conceals all the logarithms of the coupling 
constant into more complicated functions of as and therefore it performs an intrinsic resummation 
of all these contributions, as obvious, being the exact solution of the non-singlet equation at 
NNLO. A numerical implementation of the recursion relations associated to these new functions 
of the coupling constants, in this case, is no different from the previous cases, when only functions 
of the form log(a/ao) have been considered, but with a faster convergence rate. 



1.9 Truncated solutions at LO and NNLO in the singlet case 

The proof of the existence of a valid logarithmic ansatz that reproduces the truncated solution of 
the singlet DGLAP equation at NLO is far more involved compared to the non singlet case. Before 
we proceed with this discussion, it is important to clarify some points regarding some known 
results concerning these equations in moment space. First of all, as we have widely remarked 



before, there are no exact solutions of the singlet equations in moment space beyond those known 
at LO, due to the matrix structure of the equations. Therefore, it is no surprise that there is no 
logarithmic ansatz that can't do better than to reproduce the truncated solution, since only these 
ones are available analytically in moment space. If we knew the structure of the exact solution 
in moment space we could construct an ansatz that would generate by recursion relations all 
the moments of that solution, following the same strategy outlined for the non-singlet equation. 
Therefore, inverting numerically the equations for the moments has no advantage whatsoever 
compared to the numerical implementation of the logarithmic series using the algorithm that 
we have developed here. However, we can arbitrarily improve the logarithmic series in order to 
capture higher order contributions in the truncated solution, a feature that can be very appealing 
for phenomenological purposes. 

The proof that a suitable logarithmic ansatz reproduces the truncated solution of the moments 
of the singlet pdf 's at NLO goes as follows. % 



1.9.1 The exact solution at LO 

We start from the singlet matrix equation 

whose well known LO solution in Mellin space can be easily identified 

/(iV, Q,) = L(a„ ao, N)fiN, oq), (1.159) 
/ \Ro(N) 

and where L{as,ao, N) = ( §J ) is the evolution operator. 

Diagonalizing the Rq operator, in the equation above, we can write the evolution operator 

L{as, ao, N) as 

L{as,ao,N) = e+(—)^ + e.(—) , (1.160) 



where r± are the eigenvalues of the matrix Rq = {—2/ I3q)Pq and e+ and e_ are projectors [21 
[22l defined as 



1 



^0 - r^I 



r± - rzp 

Since the e± are projection operators, the following properties hold 



(1.161) 



e+e+ = e+ e_e_ = e_ e+e_ = e_e+ = e+ + e_ = 1 . (1.162) 



Based on the paper [79| 



Hence it is not difficult to see that 



Rq{N) = e+r+ + e_r_ . 



(1.163) 



It is important to note that one can write a solution of the singlet DGLAP equation in a closed 
exponential form only at LO. 

It is quite straightforward to reproduce this exact matrix solution at LO using a logarithmic 
expansion and the associated recursion relations. These are obtained from the ansatz (here 
written directly in moment space) 



subject to the initial condition 



AniN) 



n=0 



f{N,ao) = Ao{N). 



(1.164) 



(1.165) 



Then the recursion relations become 



2 



An+i{N) = - pW{N)An{N) = Ro{N)An{N), 



and can be solved as 



(1.166) 



Ari 



e+rl + e.r^)fiN,ao), 



(1.167) 



having used eq. (I1.162p . Inserting this expression into eq. (|1.164p we easily obtain the relations 



[Ro{N) 



n=0 



In 



a. 



AoiN) 



ao 



RoW 



fiN,ao) 



(1.168) 



in agreement with eq. (I1.159|) . 



1.9.2 The standard NLO solution from moment space 



Moving to NLO, one can build a truncated solution in moment space of eq. (I1.158|) by a series 
expansion around the lowest order solution. 

We start from the truncated version of the vector equation (jl.lSSp 



da. 



Potts 



pio) . ili / p(i) _ ^ p(o) 
27r V 2 



(1.169) 



that we re-express in the form 

df{N,as) 
das 



1 



1 

a. 



Rq + a. 



\ vr/3o 47r 
f{N,as), 



f{N,as 



(1.170) 



and with the i?i operator defined as 

RiiN) = - 



(1.171) 



We use [211 [22] a truncated vector solution of p.l70p - accurate at 0{as) - of the form 
f{N,as) = U{as,N)L{as,ao,N)U-^(ao,N)f{N,ao) 



l + asUi{N) L{as,ao,N) l + aoC/i(iV) f{N,ao), 



-1 



(1.172) 



where we have expanded in powers of Og the operators U{as,N) and Lf ^{ao,N). Inserting 
(|1.172p in eq. (|1.170p . we obtain the commutation relations involving the operators Ui, Rq and 



RqjUi 



Ui-Ri, 



(1.173) 



which appear in the solution in Mellin space (SUES]- Then, using the properties of the projection 
operators 

Ui = e+Uie+ + e+Uie^ + e_J7ie+ + e_C/ie_ = {7i++ + J7f + + + , (1.174) 
and inserting this relation in the commutator (|1.173|1 we easily derive the relation 

e+Rie- e-Rie+ 



e+i?ie+ + e-RiC- 



r_ — 1 r_ — r-L. — 1 



(1.175) 



Finally, expanding the term 



1 + QqC/i 



in eq. p.l72p we arrive at the solution [2ll [22] 



L + agUiL — a^LUi 



fiN,ao) 



(1.176) 



where the {us, uq, N) dependence has been dropped. Such solution can be put in a more readable 



form as 



"0 



e+ + {as — ao) e+i?ie++ 



ao - as \ — ] 



+ (+---)^/(A^,«o), (1.177) 



which can be called the standard NLO solution, having been introduced in the literature about 
20 years ago [22]. It is obvious that this solution is a (first) truncated solution of the NLO singlet 
DGLAP equation, with the equation truncated at the same order. 



1.9.3 Reobtaining the standard NLO solution using the logarithmic expan- 
sion 

Having worked out the well-known NLO singlet solution in moment space, our aim is to show 
that the same solution can be reconstructed using a logarithmic ansatz. This fills a gap in the 
previous literature on this types of ansatze for the QCD pdf's. To facilitate our duty, we stress 
once more that the type of recursion relations obtained in the non-singlet and singlet cases are 
similar. In fact the matrix structure of the equations doesn't play any role in the derivation due 
to the linearity of the ansatz in the (vector) coefficient functions that appear in it. 

Our NLO singlet ansatz has the form 

. Vao/J ^ n! [ V«o/J 

with An and -B„ now being vectors involving the singlet components. The recursion relations 

are 

An+i{N) = Ro{N)MN), 

Bn+i{N) = -Bn{N) - ^An+i{N) + Ro{N)Bn{N) - -l-p^^){N)An{N), (1.179) 

47r Trpo 

subjected to the initial condition 

/(A^,ao) = MN) + aoMN). (1.180) 

The solution of (|1.179|1 in moment space can be easily found and is given by 

^(iV) = [e+{r+r + e^r^r] M^)- (1.181) 

It has to be pointed out that the Bn{N) vectors are bidimensional column vectors and they can 
be decomposed in orthonormal basis {61,62} as 



2.NL0 sr^ An{x) 



n'. 



BniN) = eii?(i)(iV) + e2B^^\N) = S+(iV) + B-{N) 



(1.182) 



Also, we can re-arrange the Bn+i{N) relation into the form 

Bn+i{N) = (Ro - l) Bn{N) + RiR^A^iN). 



(1.183) 



The last step to follow in order to construct the truncated solution involves a projection of the 
recursion relations (|1.183p in the basis of the projectors e±, and in the basis {ei, 62}. To do this, 
we separate the equations as 



B 



n+l 



e+r^ + e_r_ 
++ 1 £>+- 



11 



+ B^' + B'+ + B' 



+ 



R^ + R-y^ + R 



-+ 
1 



R^ 



[e+Tl + e.rl)Ao, 



(1.184) 



where we have used the notation (not to be confused with the previous one) 

e+B«(iV)ei=4+ 
e+B^^\N)e2 = Bt~ 
e.B^^\N)e,=B-+ 
e.B^^\N)e2 = B--. 

Then, we can split the relation (11.1831) in four recursion relations 



(1.185) 



R++ 



B 
B 



n+l 



n+l 



1)4 

1)4 
1)4 
1)4 



+ Rt+rlAo, 
+ Rt-r'lAo, 
+ R^+rlAo, 
+ R^-r'lAo . 



Finally, it is an easy task to verify that the solutions of the recursion relations at NLO 

by 



(1.186) 
are given 



B++ = 


VI - {r+ - 


1)" 


4+^, 


Bn~ = 


[r'l - (r_ - 


1)" 


4"^o, 


4" = 


[-rl + (r+ 


-1] 


J r+ — r_ — 1 


4+ = 




-1) 


J j-_ _ j.^ _ \ 



Aq, 

Aq, 



(1.187) 



where we have expressed the nth iterate in terms of the initial conditions, and we have taken 
Bq = 0. Summing over all the projections, we arrive at the following expression for the NLO 



truncated solution of the singlet parton distributions 



f{N, as) = Y,— An + as (B++ + B'- + 3'+ + B, 



n\ L 

n=0 

which can be easily exponentiated to give 



f{N,as) = c+Aq 



a. 



+ 



a. 



Us { e+Rie+ ( — ) - e+^ie+ ( — ) 
"0/ V"o/ 



(r+-l) 



+ 



.ao/ V«o/ 



fl.l88) 



1 



(r+ — r_ — 1) 
1 

(r_ — r_(_ — 1) 



-e+Rie- ( — + e+RiC- — 
\aoJ \aoJ 

-e_i?ie+ I — + e_i2ie+ — 
,ao/ V«o/ 



+ 



^0 



(1.189) 



Finally, organizing the various pieces we obtain exactly the solution in eq. (11.1771) . We have 
therefore shown that the logarithmic ansatz coincides with the solution of the singlet DGLAP 
equation at NLO known from the previous literature and reported in the previous section. It 
is intuitively obvious that we can build with this approach truncated solutions of higher orders 
improving on the standard solution (|1.177p known from moment space, and we can do this with 
any accuracy. However, before discussing this point in one of the following sections, we want to 
show how the same strategy works at NNLO. 



1.9.4 Truncated Solution at NNLO 

At this point, to complete our investigation, we need to discuss the generalization of the procedure 
illustrated above to the NNLO case. As usual, we start from a truncated version of eq. (jl.lSSp . 
that at NNLO can be written as 



df{N,as) _ 1 



Rq + asRi + agR2 



(1.190) 



where 



f? =-- £!^,Rih Rob2 
^ vr i 27r/3o 4 16n 



(1.191) 



whose solution is expected to be of the form [19] 



f{N,as)= I + asUi{N) + a'MN) L{as,ao,N) 1 + aoUi{N) + a'oU2{N) f{N,ao), 



-1 ^ 



(1.192) 



where 



-Ro) Ui — Ui — Ri, 

Ro, 1/2] = -R2 - RiUi + 2U2. (1.193) 
Using the projectors e+,e_ in the it subspaces, one can remove the commutators, obtaining 

R^ ^ 



R^ R^ + R2 



R^ R^ + R2 



r_ — r-i- — 1 



R^ R^ 



1 



r+ — r_ — 2 
1 



-Ri Ri ~ R2 ~\~ 



-Ri ~^ — R2 ~l" 



r_ — r_)_ — 2 

and the formal solution from Mellin space can be simplified to 

f{N,as) 



R^^ R^ 
r+ — r_ — 1 

R^ R^ ^ 
r_ — r+ — 1 



(1.194) 



L + agUiL — aoLUi 

+alU2L - asaoUiLUi + a^L (uf - U2 ) /(iV, ao) • 



(1.195) 



At this point we introduce our (1-st truncated) logarithmic ansatz that is expected to reproduce 
p.l95p . Now it includes also an infinite set of new coefficients C„, similar to the non-singlet 
NNLO case 



f{N,as) = J2 



n=0 



An + OgBn + Q^Cn 



(1.196) 



Inserting the NNLO logarithmic ansatz into (I1.190p . we obtain in moment space the recursion 
relations 





= RoA 


Bn+1 


= [Ro 


Cn+1 


= [Ro 







2]C„, 



+ 



R2 + T~Rl 



RqAq 



Bn + Bn+l ) + 



Rq + Ri 



Br, 



(1.197) 



whose solution has to coincide with (|1.195p . Also in this before, we use the e± projectors 

and notice that the structure of the recursion relations for the coefficients An and Bn remain the 
same as in NLO. Therefore, the solutions of the recursion relations for An and Bn are still given 
by (|1.181|] and (I1.187|) . We then have to find only an explicit solution of the relations for the 
new coefficients Cn+i{N). 

These relations can be solved in terms of Aq, Bq and Co with the help of p.l87p . Finally, 
taking Bq = and Co = after a lengthy computation we obtain the explicit solutions for the 
projected components 



C+' 



1 



2 (r+ — r_ — — r+ — 1) 

[2(r_ - 1)" - (r+ - 2r - (r+ - 2)"r+ - rl + rl+^ + r_ ((r+ - 2)" - rl)] Aq 

[rl - 2(r+ - 1)" + (r+ - 2)"] Ao 
+lRt-' [rl - (r+ - 2)"] Ao, 



1 R^ 



a 



2 (r+ — r_ — l)(r_ — r+ — 1) 
[2(r+ - 1)" - (r_ - 2)" - (r_ - 2)"r_ - rl + rl+^ + r+ ((r_ - 2)" - rl)] Aq 

+^R--R-- [r^ - 2(r_ - 1)" + (r_ - 2)"] Ao 
+^R,- [r!! - (r_ - 2)"] Ao, 

^ '%+-(r!!-(r+-2r)(l+r_-r+) 



" 2 + r2 + r_(3 - 2r+) - 3r+ +r2 

(2(r_ - 1)" + (r_ - 1)'V_ - r1+^ 
-(r+ - 2)" + r!!(r+ - 1) - (r_ - l)"r+) 

(r!! + (r+ - 2)" + r_(r+ - 2)" 
-2(r+ - 1)" - r_(r+ - 1)" - r+(r+ - 2)" + r+(r+ - 1)")] Ao, 

C-+ = 5 \ J \Ro^{r. - r+ - 1) ((r_ - 2)" - r") 

2 + r2 +3r+ +r2 - r_(3 + 2r+) L ^ ^ + J\\ ) +) 

-R{-k{+ (-(r_ - 2)^^ + 2(r„ - 1)" + ((r_ - 2)" - (r_ - l)")r_ 

-(r_ - 2)"r+ + (r_ - l)"r+ - r!^) 

((r_ - 2)" - 2(r+ - 1)'^ + r_(r+ - l)'^ 

-(r+ - + r!^ - r_r!^ + rl+^)] Ao- 

(1.198) 

Re-inserting these solutions into the NNLO ansatz (|1.196|1 and after exponentiation one can 
show explicitely that the logarithmic solution so obtained coincides with p.l95p . Details can be 
found in appendix A. In the practical implementations of these solutions, there are two obvious 
strategies that can be followed. One consists in the implementation of the recursion relations 
as we have done in various cases above: a sufficiently large number of iterates will converge to 



the truncated solution (|1.195p . This is obtained by implementing eqs. (I1.197|) and incorporating 
them into (I1.196p . A second method consists in the direct computation of (11.195(1 in x-space 
which becomes 



f{x,as) 



L + asUi (g) L - aoL Ui 

+aj,U2 ®L- OsaoUi (g) L ® C/i + ajl ® ( C/i (g) C/i - C/2 



/(x,ao) 



(1.199) 



and with L now replaced by its operatorial (0) form 



L 



\n=0 / 



(1.200) 



An implementation of (|1.199|) would reduce its numerical evaluation to that of a sequence of LO 
solutions built around "artificial" initial conditions given by Ui (8> f{x,ao), Ui Ui f{x,ao) 
and so on. 



1.10 Higher order logarithmic approximation of the NNLO sin- 
glet solution 

The procedure studied in the previous section can be generalized and applied to obtain solutions 
that retain higher order logarithmic contributions in the NLO/NNLO singlet cases. The same 
procedure is also the one that has been implemented in all the existing codes for the singlet: 
one has to truncate the equation and then try to reach the exact solution by a sufficiently high 
number of iterates. On the other hand, x-space (non brute force) implementations are, from this 
respect, still lagging since they are only based on the Rossi-Storrow formulation [3l[5], which we 
have analized thoroughly and largely extended in this section. 

Therefore, the only way at our disposal to reach from x-space the exact solution is by using 
higher order truncates. This is of practical relevance since our algorithm allows to perform 
separate checks between truncated solutions of arbitrary high orders built either from Mellin or 
from x-space 0. We recall, if not obvious, that in the analysis of hadronic processes the two 
criteria of using either truncated or exact solutions are both acceptable. 

To summarize: since in the singlet case is not possible to write down a solution of the DGLAP 
equation in a closed exponential form because of the non-commutativity of the operators Ri, the 
best thing we can do is to arrange the singlet DGLAP equation in the truncated form, as in 
eqs. p.lTOp and (ll.l9Qp . The truncated vector solutions (|1.195|1 and p.l76p are equivalent to 
those obtained using the vector recursion relations at NLO/NNLO. 

The current benchmarks available at NNLO are limited to exact solutions and do not involve comparisons 
between truncated solutions 



Now, working at NNLO, we will show how the basic NNLO solution can be improved and 
the higher truncates identified. Clearly, it is important to show explicitely that these truncates, 
generated after solving the recursion relations, can be rewritten exactly in the form previously 
known from Mellin space. We are going to show here that this is in fact the case, although 
some of the explicit expressions for the higher order coefficient functions Cn,D'^s will be given 
explicitely only in part. The expressions are in fact slightly lengthy [f]. Therefore, here we will 
just outline the procedure and illustrate the proof up to the second truncate of the NNLO singlet 
only for the sake of clarity. 

The exact singlet NNLO equation in Mellin space is given by 



da. 



(^) p(0)(iV) + (^)'P«(iV) + m'P^'HN) 



f{N,as) 



(1.201) 



where we have introduced the singlet kernels. After a Taylor expansion of P 
up to it becomes 



NNLO 



(iV,a,)//3^^^0(a.) 



da. 



1 

as 



Rn 



(47r 



Ri + asRi + agR2 



f{N,as), 



(1.202) 



which is the truncated equation of order . The Ri {i = 0, 1, 2) operators are listed below 



Ri 

R2 



p(l) 



+ 



+ 



The formal solution of this equation can be written as [20|.[23| 

62 



f{N,as 



exp 



da[ 



' a' 



Ro 



Ri + a'gRi + a^R2 



(47r)2 

U{N, as)L{as, ao)U~\N, ao)f{N, ao), 



p(0)_ 



^1 /3A 

— a„ R2 



(1.203) 



f{N,ao) 
(1.204) 



where the Tq, operator acts on the exponential similarly to a time-ordered product, but this time 
in the space of the couplings. Again, expanding the U{N,as) and U~^{N,ao) operators in the 



They will be included in a file that will be made available in the same distribution of our code, which is in 
preparation. 



formal solution around = and uq = 0, we have 
f{N,as) 



L + ttgUiL — aoLUi + alU2L — agaoUiLUi + a^L (tli —IJ2 



-alL (Uf - U1U2 - U2U1 + U3)] m, ao). 



(1.205) 



Inserting the expanded solution into eq. (|1.205p and equating the various power of a<j we arrive 
at the following chain of commutation relations 



Rq,Ui — Ui — Ri, 



Ro, U2 
Ro, U3 



-R2- R1U1 + 2U2, 

Ri + ^R2-RiU2-R2Ui + 3U3- 



(47r)2"^ ' (47r 

Removing the commutators by the e± projectors one obtains 



(1.206) 



U2^ 



Ri^, 
Ri , 



r+ — r_ — 1 

Ri^ 

r_ — r+ — 1 
R^^ R^^ + R^ 

R^ R^ + R2 
1 



R^ R^ ^ 
r_ — r+ — 1 

R^ ^ Ri 
r+ — r_ — 1 



r+ — r_ — 2 
1 

r_ — r-i- — 2 

1 
3 
1 
3 



~Ri Ri 



-R\ ^ Ri^ 



R+- + 

R2^ + 



R^^ R^ 
r_|. — r_ — 1 

r_ — r-L — 1 



RP 



R^ 



1 



r+ — r_ — 3 
1 



(4ir) 1 



r_ — r I 



(4^^2 -2 
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1 ^2 



12 ^1 



(47r) 1 (47r) 



2 -^2 ^ ~ -^1 ^^2^^ ~ -^1 f^2 ^ ~ -^2 ^^1^^ 



R^Ur 



Rru{ 



As one can see, the tl^ operator is expressed in terms of the kernels P^''), p(^) and P^'^\ One can 
prove that by imposing a higher order ansatz in Melhn space of the form 



U 



n=0 



fiN, a,) = V — AniN) + asBn{N) + alCn{N) + a^DniN) 



the solution (|1.205p is generated. The vector recursion relations in this case become 

2 



Cn+1 



/3o 

-Br, 



1 



27r2/3o 



/?2 



(47r)2/3o 



A 



n+l 



1 

7r/3o 



47r/3o 47r/3o Po 



1 



27r2/5( 



P2 



P2 



Pi ^ 



27r/?o 



(4vr)/3o /3o 



fl.208) 



(1.209) 



Applying the properties of the e± operators and decomposing into the {ei,e2} basis, we project 
out the -B^^_^, C^_^_l ■ ■ ■ components of these relations, and imposing the initial conditions = 
Co = -Do = together with 



fiN,ao)=Ao, 



(1.210) 



we obtain the explicit form of f{N,as). In order to construct the f{N,as) solution, the four 
lb projections of Dn-\-i must be solved with respect to ^o(-^); since the other projections B^, 
are already known. A direct computation shows that the structure of the solution can be 



organized as follows in terms of the components of Ri 

D++{N) = Wi(i?++3, r+, r_, AT, Aq) + W2{Rt~ Rf^ ,r+,r-, N, Aq) 

+Ws{Ri-R^-R^+,r+,r_,N, Ao) + Wi{R++ R+- ,r+,r., N, Ao) 
+W5(i?++i?2^+, r+, r_, iV, Ao) + We{Rt^ Rf^ ,r+,r., N, Ao) 
+ W7 {Rt-R2^: r+ , r_ , iV, lo ) + Ws {RfRl^ ,r+,r^,N,Ao) 
+ W9 {R++ ,r+,r-,N,Ao) + Wio + ,r+,r^,N,Ao), 

D--{N) = Wi(4~ ^ r+, r-,N, Ao) + W2(i?r+i^i^"^r"> ''-^ ^, ^) 

+W3(^r^^i"^-R]^"> '^-> ^> ^0) + W4{R~R^+R+-,r+,r-,N, Aq) 
+W5{Ri-R^-,r+,r-,N, Ao) + W6(^2 --Rf", ^, ^) 
+W7(^r''^^", ^) + W8(4-+^+-, r+, r_, iV, lo) 

+ W9 {R^- ,r+,r-,N,Ao) + Wio (^2 ~ , r+ , r_ , iV, Iq) , 

D^- (N) = Zi {R++R+'R^- , r+, r„ , iV, Ao) + Z2 , r+ , r_ , A^, Iq) 
+Z3(i?|-i?r"^r" > r+ , r_ , iV, lo) + Z4 (i^^- Ri+Rt-, r+,r^, N, Ao) 
+Z5 (i?^- i22"- , r+,r_,N,Ao) + Ze(Rt- R^' , r+ , r_ , iV, ) 
+Z7(i?^+«^- , r+, r_ , iV, lo) + ZsiRt^Rf- ,r+,r^,N, Ao) 
+Z9 ,r+,r-,N,Ao) + Zio{Rt~,r+,r-,N,Ao), 

D-+{N) = Z,{R~R^+Rt+,r+, r^,N, Ao) + Z2{RT~R~R^+, r+,r-,N, Ao) 
+Z3{Rj+Rt+Rt+,r+, r-,N, Ao) + Zi{R^+R+-R^+,r+,r^,N, Ao) 
+Z5(^J-+i2++, r+, r_, iV, Ao) + Zg{R^+R++ ,r+, r^,N, Ao) 
+ZriRrR2^, r+, r-,N, Ao) + Z^{iq-R{+, r+, r_, iV, Iq) 
+Z9(^]-+, r+, r_, A^, Ao) + Zio(^2 ^> ^> ^0), 

(1.211) 

where we used the notation Rl^"^ = R^^ R^^ R'^'^ . The expressions of the functions W and Z's 
can be extracted by a symbohc manipulation of the coefficients D:^^. We have included one of 
the projections for completeness in an appendix for the interested reader. 

In x-space, the structure of the (2nd) truncated (or af) NNLO solution can be expressed as 



a sequence of convolution products of the form 

f^^aff{x,as) = L{x) + asUi{x) ® L{x) - aoL{x) Ui{x) + a'iU2{x) (g> L{x) 
-asaoJJiix) ® L{x) ® tji{x) + alL{x) ® {Ui{x) (g) Ui{x) - U2{x) 
+alU3{x) L{x) + asalUi{x) L{x) (^Ui{x) Ui{x) - (J2{x) 

—0^00^72(2;) ® L{x) (8) Ui{x) — aoL(x)(8' 
Uiix) Ui{x) (g> Ui{x) - Ui{x) (g> U2{x) - U2{x) ® Ui{x) + U^ix) 



® /(x,ao), 
(1.212) 



and is reproduced by the logarithmic expansion 



U 



n=0 



/(x, a^) = — An{x) + asBn{x) + a^C„(x) + alDn{x) 
^-^ n! 



with the initial condition 



(1.213) 



/(x,ao) = Aq{x). 



(1.214) 



The study of higher order truncates is performed numerically with the implementation of the 
generalized recursion relations (ll.QOp given in the previous sections. 



1.11 Comparison with existing programs 



In this section we present a numerical test of our solution algorithm. For this aim, we compare 
the results of a computer program that implements our method with the results of QCD-Pegasus 
[2], a PDF evolution program based on Mellin-space inversion, which has been used by the QCD 
Working Group to set some benchmark results [SUES]- In the following, we refer to these results 
as to the benchmark. 

In the tables we are going to show in this chapter, we set the renormalization and factorization 
scales to be equal, and we adopt the fixed flavor number scheme. The final evolution scale is 
ji'^p = lO^GeV^. We limit ourselves to this case because it is enough to test the reliability of 
our method. A more lenghty and detailed analysis, which will take into account many other 
cases with renormalization scale dependence and variable flavor number scheme, is planned to 
be presented for the future. 

As in the published benchmarks, we start the evolution at //|.q = 2GeV^, where the test 
input distributions, regardless of the perturbative order, are parametrized by the following toy 



model 



XUy{x) 

xg{x) 
xd{x) 
xu{x) 



5.107200j;°-'^(1 -x)^ 
3.064320j;°-'^(1 -x)"^ 
1.700000x"°-^(l - j;)^ 
0.1939875x"°-i(l -x) 
(1 — x)xd{x) 



6 



xs{x) = xs{x) 



^.2x{u + d){x) 



(1.215) 



and the running coupling has the value 



as{li%o = 2GeV2) = 0.35. 



(1.216) 



We remind that Qv = Q — Q, Q+ = Q + Q, L± = d± u. Our results are obtained using the exact 
solution method for the nonsinglet and the LO singlet, and the K-th truncate method for the 
NLO and NNLO singlet, with k = 10. In each entry in the tables, the first number is our result 
and the second is the difference between our results and the benchmark. 

In Table 11.31 we compare our results at leading order with the results reported in Table 2 of 
[24| . The agreement is excellent for any value of x for the nonsinglet (xu^ and xdy); regarding 
the singlet (xg column), the agreement is excellent except at very high x: we have a sizeable 
difference at x = 0.9. In Table 11.41 we analyze the next-to-leading order evolution; the results 
for the proposed benchmarks are reported in Table 3 of [24]. The agreement is very good for 
any value of x for the nonsinglet; for the singlet the agreement is good, except at very high x 



Moving to the NNLO case (Table [LSl the benchmarks are shown in Table 14 of [25]), some 
comments are in order. We don't solve the non-singlet equation as in PEGASUS, since (I1.125|l 
admits an exact solution (|1.127p . which in PEGASUS is obtained only by iteration of truncated 
solutions. Our implementation is based on the exact solution presented in this thesis. The 
discrepancy between our results and PEGASUS are of the order of few percent, and they become 
large in the gluon case at x = 0.9, as for the lower orders. Another comments should be made 
for the sea asimmetry of the s quark, which is nonvanishing at NNLO. In this case we have a 
sizeable relative discrepancy for any value of x (reported in the column xs„), but it is evident 
that this asymmetry is quite small, especially if compared with the column xs+, whose entries are 
several orders of magnitude larger. This means that xs and xs should be comparable and their 
difference therefore pretty small. We have to notice that xs^ cannot be computed directly by a 
single evolution equation. Indeed, it is computed by a difference between very close numbers, a 
procedure that amplifies the relative error. 



(x = 0.9). 



LO, uf =4, fij,^ ^1 = 10-1 GeV' 


X 






xL- 


2xL+ 


XS+ 




xg 


10-' 


5.7722 • lO-'' 
0.0000-10-5 


3.4343 - lO-'' 
0.0000-10-5 


7.6527 - 10- ' 
0.0000 - 10-^ 


9.9465 - 10+1 
0.0000 - 10+1 


4.8642 - 10+1 
0.0000 - 10+1 


4.7914 - 10+1 
0.0000 - 10+1 


1.3162 - 10+3 
0.0000 - 10+3 


10-*^ 


3.3373 • 10-1 
0.0000- lo-'' 


1.9800 - lO-'' 
0.0000-10-1 


5.0137 - 10-*^ 
0.0000 - lo-*' 


5.0259 - 10+1 
0.0000 - 10+1 


2.4263 - 10+1 
0.0000 - 10+1 


2.3685 - 10+1 
0.0000 - 10+1 


6.0008 - 10+2 
0.0000 - 10+2 


lO-'"" 


1.8724 • 10"^ 
0.0000 • 10-3 


1.1065 - 10-3 
0.0000 - 10-3 


3.1696 - 10-5 
0.0000-10-5 


2.4378 - 10+1 
0.0000 - 10+1 


1.1501 - 10+1 
0.0000 - 10+1 


1.1042 - 10+1 
0.0000 - 10+1 


2.5419 - 10+2 
0.0000 - 10+2 


10-1 


1.0057- 10-^ 
0.0000-10-2 


5.9076 - 10-3 
0.0000 - 10-3 


1.9071 - lO-"' 
0.0000 - 10-1 


1.1323 - 10+1 
0.0000 - 10+1 


5.1164 - 10+° 
0.0000 - 10+° 


4.7530 - 10+° 
0.0000 - 10+° 


9.7371 - 10+1 
0.0000 - 10+1 


10-^ 


5.0392 - 10-2 
0.0000-10-2 


2.9296 - 10-2 
0.0000-10-2 


1.0618 - 10-3 
0.0000 - 10-3 


5.0324 - 10+1 
0.0000 - 10+° 


2.0918 - 10+° 
0.0000 - 10+° 


1.8089 - 10+° 
0.0000 - 10+° 


3.2078 - 10+1 
0.0000 - 10+1 


10-^ 


2.1955- 10-1 
0.0000-10-1 


1.2433- 10-1 
0.0000-10-1 


4.9731 - 10-3 
0.0000 - 10-3 


2.0433 - 10+° 
0.0000 - 10+° 


7.2814- 10-1 
0.0000-10-1 


5.3247- 10-1 
0.0000-10-1 


8.0546 - 10+° 
0.0000 - 10+° 


0.1 


5.7267-10-1 
0.0000-10-1 


2.8413- 10-1 
0.0000-10-1 


1.0470-10-2 
0.0000-10-2 


4.0832- 10-1 
0.0000-10-1 


1.1698-10-1 
0.0000-10-1 


5.8864-10-2 
0.0000 - 10-2 


8.8766-10-1 
0.0000-10-1 


0.3 


3.7925-10-1 
0.0000-10-1 


1.4186- 10-1 
0.0000-10-1 


3.3029 - 10-3 
0.0000 - 10-3 


4.0165- 10-2 
0.0000 - 10-2 


1.0516- 10-2 
0.0000-10-2 


4.1380- 10-3 
+0.0001 - 10-3 


8.2676-10-2 
0.0000-10-2 


0.5 


1.3476-10-1 
0.0000-10-1 


3.5364-10-2 
0.0000-10-2 


4.2815-10-1 
0.0000 - lo-i 


2.8624- 10-3 
0.0000 - 10-3 


7.3137-10-1 
-0.0001 - 10-1 


2.6481 - lO-'i 
0.0000 - 10-"^ 


7.9242 - 10-3 
+0.0002-10-3 


0.7 


2.3123- 10-2 
0.0000-10-2 


3.5943-10-3 
0.0000-10-3 


1.5868-10-5 
0.0000-10-5 


6.8970-10-5 
+0.0009 - 10-5 


1.7730-10-5 
+0.0005 - 10-5 


6.5593 - 10-° 
+0.0044 - 10-° 


3.7301 - 10-1 
-0.0010-10-1 


0.9 


4.3443 - 10-1 
0.0000 - 10-1 


2.2287- 10-5 
0.0000-10-5 


1.1042-10-** 
0.0000 - 10-8 


3.3030 - 10-» 
-0.3263- 10-8 


8.5607- lO-'^ 
-0.1631 - 10-8 


3.2577- 10-*^ 
-1.6316- 10-9 


1.3887-10-° 
+0.2969-10-° 



Table 1.3: Comparison between our and the benchmark results at LO. In each entry, the first number is our result and the second one is the 
difference between our result and the benchmark. The benchmark values are reported in Table 2 of [21] 



NLO, n/ = 4, /4 = /4 = 10* GeV^ 


X 


XUy 


xdv 


xL^ 


2xL+ 


xs^ 


xc^ 


xg 




1.0620- 10-4 
+0.0004-10-4 


6.2353- lO-'^ 
+0.0025 - 10-'^ 


4.2455 - lO-t* 
+0.0015-10-6 


1.3710-10+2 
+0.0112- 10+2 


6.7469- 10+1 
+0.0556- 10+1 


6.6750 - 10+1 
+0.0555 - 10+1 


1.1517- 10+3 
+0.0034-10+3 


10-" 


5.4196- 10-4 
+0.0019-10-4 


3.1730- 10-* 
+0.0011 - 10-4 


1.9247- lO-'' 
+0.0006 - 10-5 


6.8896 - 10+1 
+0.0500-10+1 


3.3592-10+1 
+0.0250-10+1 


3.3021 - 10+1 
+0.0250 - 10+1 


5.4048 - 10+2 
+0.0137-10+2 


10-^ 


2.6878-10-^ 
+0.0008 - 10-3 


1.5682-10-'^ 
+0.0005 - 10-3 


8.3598-10-5 
+0.0023-10-5 


3.2936-10+1 
+0.0208-10+1 


1.5788-10+1 
+0.0103-10+1 


1.5335-10+1 
+0.0104-10+1 


2.3578 - 10+2 
+0.0050-10+2 


10-* 


1.2844- 10-^ 
+0.0003 - 10-2 


7.4576 - 10-3 
+0.0018-10-3 


3.4919- 10-4 
+0.0008 - 10-4 


1.4824-10+1 
+0.0078 - 10+1 


6.8744 - 10+" 
+0.0389 - 10+" 


6.5156 - 10+" 
+0.0387- 10+" 


9.3026- 10+1 
+0.0154- 10+1 


10-^ 


5.7937-10-2 
+0.0011 - 10-2 


3.3343 - 10-2 
+0.0006 - 10-2 


1.4164- 10-3 
+0.0002 - 10-3 


6.1899- 10+" 
+0.0251 - 10+° 


2.6783 - 10+" 
+0.0124- 10+" 


2.4001 - 10+" 
+0.0123- 10+" 


3.1540- 10+1 
+0.0038 - 10+1 


10-^ 


2.3029- 10-^ 
+0.0003-10-1 


1.2930- 10-1 
+0.0002-10-1 


5.3258 - 10-3 
+0.0007-10-3 


2.2587 - 10+" 
+0.0060 - 10+" 


8.4518-10-1 
+0.0298-10-1 


6.5540-10-1 
+0.0294- 10-1 


8.1120- 10+" 
+0.0054- 10+" 


0.1 


5.5456- 10-1 
+0.0004-10-1 


2.7338-10-1 
+0.0002-10-1 


1.0012- 10-2 
+0.0001 - 10-2 


3.9392-10-1 
+0.0056-10-1 


1.1517- 10-1 
+0.0028-10-1 


6.0619-10-2 
+0.0268-10-2 


8.9872-10-1 
+0.0005-10-1 


0.3 


3.5395- 10-1 
+0.0002 - 10-1 


1.3158-10-1 
0.0000-10-1 


3.0363-10-3 
+0.0001 - 10-3 


3.5884-10-2 
+0.0036 - 10-2 


9.2210- 10-3 
+0.0180-10-3 


3.4066- 10-3 
+0.0176-10-3 


8.3415-10-2 
-0.0036 - 10-2 


0.5 


1.2271 - 10-1 
0.0000-10-1 


3.1968-10-2 
+0.0001 - 10-2 


3.8266- 10-* 
+0.0001 - 10-4 


2.4149- 10-3 
+0.0023 - 10-3 


5.8539-10-* 
+0.0115-10-4 


1.7068- 10-* 
+0.0113- 10-* 


8.0412- 10-3 
-0.0061 - 10-3 


0.7 


2.0429- 10-2 
0.0000-10-2 


3.1474- 10-3 
+0.0001 - 10-3 


1.3701 - 10-5 
0.0000-10-5 


5.3703 - 10-5 
+0.0081 - 10-5 


1.2432 - 10-5 
+0.0039 - 10-5 


2.8201 - 10-" 
+0.0394- 10-" 


3.8654- 10-* 
-0.0067- 10-* 


0.9 


3.6097-10"'' 
+0.0001 - 10-4 


1.8317- 10-^ 
0.0000-10-5 


8.9176-10-9 
-0.0054- 10-9 


1.6196- 10-** 
-0.4724-10-^ 


1.6717- 10-9 
-2.3673- 10-9 


-2.6084- 10-9 
-2.3681 - 10-9 


1.8308-10"'' 
+0.6181 - 10-'^ 



Table 1.4: Same as in Table [L3] in the NLO case. The benchmark values are re] 



ported in 



Table 



3 of [24] 



NNLO, uf ^ A, ^ fi% = lO^GeV^ 








xL^ 


2xL+ 




XS+ 


XC+ 


xg 


10-' 


1.4069 • 10--* 
-0.1218- 10-* 


9.0435 - 10-5 
-0.1201 - IQ-'' 


5.5759- 10-5 
-0.1259- 10-5 


1.4184-10+2 
+0.0994 - 10+2 


1.9569- 10-5 
-1.1868-10-5 


6.9844- 10+1 
+0.4967-10+1 


6.9127-10+1 
+0.4966 - 10+1 


1.0519- 10+3 
+0.0543 - 10+3 


10-" 


6.5756 • lO-'^ 
-0.3420- lO-'' 


4.0826- 10-1 
-0.3458 - 10-1 


2.4722-10-5 
-0.0688-10-5 


7.1794-10+1 
+0.3295 - 10+1 


5.8862-10-5 
-3.5417-10-5 


3.5043-10+1 
+0.1646- 10+1 


3.4474 - 10+1 
+0.1646-10+1 


5.1093-10+^ 
+0.1969-10+2 


10-5 


3.0260 - 10-^ 
-0.0721 - 10-3 


1.8199- 10-3 
-0.0775 - 10-3 


1.0393-10-1 
-0.0326- lO"'' 


3.4373 - 10+1 
+0.0902-10+1 


1.4230- 10-1 
-0.8560-10-1 


1.6509- 10+1 
+0.0450-10+1 


1.6057-10+1 
+0.0450 - 10+1 


2.2899 - 10+2 
+0.0602 - 10+2 


10--* 


1.3656-10-^ 
-0.0066-10-2 


8.0052-10-3 
-0.0967-10-3 


4.1299- 10-1 
-0.1259- 10-1 


1.5403 - 10+1 
+0.0199-10+1 


2.2837- 10-1 
-1.3807-10-1 


7.1661 - 10+" 
+0.0991 - 10+" 


6.8085 - 10+" 
+0.0988 - 10+" 


9.2125- 10+1 
+0.1457-10+1 


10-^ 


5.9360-10-^ 
+0.0200-10-2 


3.4135-10-^ 
+0.0085 - 10-2 


1.5650 - 10-3 
-0.0358-10-3 


6.3657-10+^ 
+0.0427 - 10+° 


8.9572-10-5 
-0.5522- 10-1 


2.7684 - 10+" 
+0.0210- 10+" 


2.4913 - 10+" 
+0.0209 - 10+" 


3.1592- 10+1 
+0.0243- 10+1 


10-^ 


2.3139- 10-1 
+0.0061 - 10-1 


1.2958-10-1 
+0.0039-10-1 


5.5064- 10-3 
-0.0624- 10-3 


2.2868 - 10+1 
+0.0116- 10+" 


-3.5702-10-1 
+2.1611 - 10-1 


8.6094-10-1 
+0.0592-10-1 


6.7224- 10-1 
+0.0601-10-1 


8.1503 - 10+" 
+0.0122- 10+" 


0.1 


5.5125-10-1 
-0.0052-10-1 


2.7142- 10-1 
-0.0023-10-1 


9.9834-10-3 
-0.0040- 10-2 


3.9119-10-1 
+0.0100-10-1 


-1.9045- 10-1 
+1.1582- 10-1 


1.1453- 10-1 
+0.0067-10-1 


6.0520-10-2 
+0.0747-10-2 


8.9909-10-1 
-0.0654- 10-1 


0.3 


3.5017-10-1 
-0.0054-10-1 


1.3005-10-1 
-0.0020-10-1 


3.0025 - 10-3 
-0.0073-10-3 


3.4975 - 10-2 
-0.0383-10-2 


-1.9830-10-5 
+1.2061 - 10-5 


8.8758 - 10-3 
-0.1722- 10-3 


3.1421 - 10-3 
-0.1640- 10-3 


8.3041 - 10-2 
-0.1145- 10-2 


0.5 


1.2099- 10-1 
-0.0018-10-1 


3.1485-10-^ 
-0.0043 - 10-2 


3.7667-10-1 
-0.0075 - 10-1 


2.1876- 10-3 
-0.1991 - 10-3 


-1.6924- 10-5 
+1.0291 - 10-5 


4.8155-10-1 
-0.9810-10-1 


7.4120- 10-5 
-0.9758-10-1 


7.9784-10-3 
-0.1342- 10-3 


0.7 


2.0052- 10-2 
-0.0025-10-2 


3.0849-10-3 
-0.0037-10-3 


1.3411 - 10-5 
-0.0023- 10-5 


1.5984-10-5 
-3.8260 - 10-5 


-6.2854- 10-** 
+0.3821 - lO-'^ 


-6.1500- 10-5 
-1.9084- 10-5 


-1.5545- 10-5 
-1.9075-10-5 


3.8226- 10-1 
-0.0722- 10-1 


0.9 


3.5078 - 10-1 
-0.0033-10-1 


1.7767-10-5 
-0.0016-10-5 


8.6326 - 10-** 
-0.0184- 10-3 


-6.4293 - 10-^ 
-6.6986-10-^ 


-9.1828- 10-11 
+0.5579-10-1" 


-3.2776- 10-' 
-3.3487-10-^ 


-3.3190- 10-' 
-3.3487-10-^ 


1.9280- 10-5 
+0.7144- 10-5 



Table 1.5: Same as in Table [L3] in the NNLO case. The benchmark values are reported in Table 14 of [25] 



1.12 Future objectives 

We have shown that logarithmic expansions identified in x-space and implemented in this space 
carry the same information as the solution of evolution equations in Mellin space. This has been 
obtained by the introduction of new and generalized expansions that we expect to be very useful 
in order to establish benchmarks for the evolution of the pdf's at the LHC. Our analysis has 
been presented up to NNLO. We have also shown how exact expansions can be derived. We have 
presented analytical proofs of the equivalence and clarified the role of previous similar analysis 
which were quite limited in their reach. We have also presented a numerical comparison of our 
results against those obtained using PEGASUS, for a specific setting. The overall agreement, 
as we have seen, is very good down to very small x-values. One future objective will be a more 
detailed analysis based on the numerical implementation of our results with various comparisons 
between our approach and other approaches with different set up. 

There are several issues which are still unclear in this area and concern the role of the 
NNLO effects in the evolution and in the hard scatterings, the role of the theoretical errors in 
the determination of the pdf's, whether they dominate over the NNLO effects or not, and the 
impact of the choices of various truncations in the determination of the numerical solution of the 
pdf's, along the lines of our work. Similar analysis can be performed by other methods, but we 
think that it is important, in the search for precise determination of cross sections at the LHC, 
to state clearly which algorithm is implemented and what accuracy is retained, with a particular 
attention to the issues connected to the resummation of the perturbative expansion |26j . Our 
work, here, has been limited to a (fixed order) NNLO analysis. We hope that our analysis has 
shown that x-space approaches have a very solid base and provide a simple view on the structure 
of the solutions of the DGLAP equations, valid to all orders. 

1.13 Appendix A. Derivation of the recursion relations at NNLO 

As an illustration we have included here a derivation of the recursion relations for the first 
truncated ansatz of 0{a'^) that appears at NNLO. 

Inserting the NNLO truncated ansatz for the solution into the DGLAP equation we get at 
the left-hand-side of the defining equation 



Note that the first sum starts at n = 1, because the n = term in (I1.217P does not have a 




n=0 



(1.217) 



dependence. Sending re ^ ?i + 1 in the first sum, using the three-loop expansion of the beta 
function (|1.2p and neglecting all the terms of order or higher, the previous formula becomes 



n=0 



re! 



47r 



f3l 2 



167r2 



* a! 



647r3 



^ n! [ Att ' 167r2 



^ re! V 47r^ 167r2^^ 



A, 



re! 



-a. 



re! 



47r 



(1.218) 



At this point we use the NNLO expansion of the kernels. We get at the right-hand-side of the 
defining equation 



E 

n=0 



L" (a 



re! 27r 



+ 



+ 



87r3 
4^ 



{x) + 



4^ 



a, 
27r 



(X) 

(x) 
(x) 



(1.219) 



Equating p.218p and p.219p term by term and grouping the terms proportional respectively to 
Os, and we get the three desired recursion relations p.61|) . Setting Q = Qo in (I1.58P we 
get 



fix, Ql) = Aoix) + as{Ql)Bo{x) + (asiQ^))' Co{x). 
We have seen that the initial conditions should be chosen as 



(1.220) 



Bo{x) = Coix) = 0, fix, Ql) = Aoix) 



fl.221) 



in order to reproduce the moments of the truncated solution of the DGLAP equation. 



1.14 Appendix B. NNLO singlet truncated solution 



Putting all the projections of the coefficients Cn into the NNLO singlet ansatz we get 

An + as (b++ + B-- + B-+ + S+ 



^ jn 



n=0 



+- 



(1.222) 



and exponentiating 



/(iV,a,) = e+Ao ('^y +e_Io ( -] ' 



asle+Rie+(^^ - e+Rie+ 



V "0 / 



(^+-1) 



+ 



e-Rie- 



— e-RiC- 



(r_-l) 



+ 
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1 
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This expression is equivalent to that obtained from eq. (11.195(1 . Projecting over all the it 
components, and formally introducing the quantities q{N, as)~^^ ,q{N, a^)"' , . . . , we can write 
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Using Aq = f{N, ao) = q{N, uq) it is simple to obtain that 
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For example we can check the first of the relations above, gives 
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Factorizing {as/aoY+ and expanding the power of the previous expression becomes 
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1.15 Appendix C. Calculation of 

An explicit calculation of the vector coefficient D^^ of the k = 4 (4th truncated) solution of the 
NNLO singlet equation has been done in this section. Since the expressions of the coefficients 
, D~~^ and D~~ have a structure similar to , we omit them and give only the explicit 



form of this one 
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40 6i7rE^+ r_ r+ -4^2 r.^ r+ + 327r^ Rt"'^rJr+ - 166i7rE^+ r_V++ 
240.^i.-3,_ ^Z|±a)%, _ 9e.^^-3,_3 (_|±I^y _ 
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Chapter 2 



Applications: Solving the x-space 
Evolution Equations for Transversity at 
NLO 

2.1 Introduction 

One of the most fascinating aspects of the structure of the nucleon is the study of the distribution 
of spin among its constituents, a topic of remarkable conceptual complexity which has gained 
a lot of attention in recent years. This study is entirely based on the classification and on the 
phenomenological modeling of all the leading-twist parton distributions, used as building blocks 
for further investigations in hadronic physics. 

There are various theoretical ways to gather information on these non-local matrix elements. 
One among the various possibilities is to discover sum rules connecting moments of these distribu- 
tions to other fundamental observables. Another possibility is to discover bounds - or inequalities 
- among them and use these results in the process of their modeling. There are various bounds 
that can be studied, particularly in the context of the new generalized parton dynamics typical of 
the skewed distributions [271 [28]. All these relations can be analized in perturbation theory and 
studied using the Renormalization Group (RG), although a complete description of their per- 
turbative dynamics is still missing. This study, we believe, may require considerable theoretical 
effort since it involves a global understanding both of the (older) forward (DGLAP) dynamics 
and of the generalized new dynamics encoded in the skewed distributions. 

In this context, a program aimed at the study of various bounds in perturbation theory using 
primarily a parton dynamics in x-space has been outlined [29]. This requires accurate algorithms 
to solve the equations up to next-to-leading order (NLO). Also, underlying this type of description 
is, in many cases, a probabilistic approach [30] which has some interesting consequences worth of 
a closer look . In fact, the DGLAP equation, viewed as a probabilistic process, can be rewritten 
in a master form which is at the root of some interesting formal developements. In particular, a 
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wide set of results, available from the theory of stochastic processes, find their way in the study 
of the evolution. In a recent work [29] it has been proposed a Kramers-Moyal expansion of the 
DGLAP equation as an alternative way to describe the dynamics of parton evolution. Here, this 
analysis will be extended to the case of the non-forward evolution. 

With these objectives in mind, in this study we test x-space algorithms up to NLO developed 
in the previous chapter, and verify their accuracy using a stringent test: Soffer's inequality. As 
usual, we are bound to work with specific models of initial conditions. The implementations on 
which our analysis are based are general, with a varying flavour number n/ at any threshold of 
intermediate quark mass in the evolution. Here, we address Soffer's inequality using an approach 
based on the notion of "superdistributions" [31], which are constructs designed to have a simple 
(positive) evolution thanks to the existence of an underlying master form [30l[29]. The original 
motivation for using such a master form (also termed kinetic or probabilistic) to prove positivity 
has been presented in [3l], while further extensions of these arguments have been presented 
in [29]. In a final section we propose the extension of the evolution algorithm to the case of 
the skewed distributions, and illustrate its implementation in the non-singlet case. As for the 
forward case, numerical tests of the inequality are performed for two different models. We show 
that even starting from a saturated inequality at the lowest evolution scale, the various models 
differ significantly even for a moderate final factorization scale of Q = 100 GeV. Finally, we 
illustrate in another application the evolution of the tensor charge and show that, in the models 
considered, differences in the prediction of the tensor charge are large. 



2.2 Prelude to x-space: A Simple Proof of Positivity of hi to NLO 

There are some nice features of the parton dynamics, at least in the leading logarithmic approx- 
imation (LO), when viewed in x-space, once a suitable "master form" of the parton evolution 
equations is identified. 

The existence of such a master form, as firstly shown by Teryaev, is a special feature of the 
evolution equation itself. The topic has been addressed before in LO [M] and reanalized in more 
detail in [29] where, starting from a kinetic interpretation of the evolution, a differential equation 
obtained from the Kramers-Moyal expansion of the DGLAP equation has also been proposed. 

The arguments of refs. [HUE^ are built around a form of the evolution equation which has a 
simple kinetic interpretation and is written in terms of transition probabilities constructed from 
the kernels. 

The strategy used, at least in leading order, to demonstrate the positivity of the LO evolution 
for special combinations of parton distributions Q± [31], to be defined below, or the NLO evo- 
lution for /ii, which we are going to address, is based on some results of ref.[3T], briefly reviewed 
here, in order to be self-contained. 



A master equation is typically given by 

—f{x,T)= j dx' {w{x\x')f{x' ,t) — w{x'\x)f{x,T)) dx' 

and if through some manipulations, a DGLAP equation 

dq{x,Q^) _ c^y ^2 



(2.1) 



dlog(Q2) y 
with kernels -P(x), is rewritten in such a way to resemble eq. (|2.1 



P{x/y)q{y,Q'), (2.2) 



^qix, r) = C dyP ( A ^J^lllll _ T (1) ^ 
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with a (positive) transition probability 
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w{x\y) = ^P(x/y)^^^^ (2.4) 
ivr y 

then positivity of the evolution is established. 

For equations of non-singlet type, such as those evolving = q — q, the valence quark 
distribution, or hi, the transverse spin distribution, this rewriting of the equation is possible, 
at least in LO. NLO proofs are, in general, impossible to construct by this method, since the 
kernels turn out, in many cases, to be negative. The only possible proof, in these cases, is 
just a numerical one, for suitable (positive) boundary conditions observed by the initial form of 
the parton distributions. Positivity of the evolution is then a result of an unobvious interplay 
between the various contributions to the kernels in various regions in x-space. 

In order to discuss the probabilistic version of the DGLAP equation it is convenient to 
separate the bulk contributions of the kernels (a; < 1) from the edge point contributions at 
X = 1. For this purpose we recall that the structure of the kernels is, in general, given by 

P{z) = P{z)-6{l-z) [ P{z)dz, (2.5) 

Jo 

where the bulk contributions {z < 1) and the edge point contributions (~ 5{z — 1)) have been 
explicitely separated. We focus on the transverse spin distributions as an example. With these 
prerequisites, proving the LO and NLO positivity of the transverse spin distributions is quite 
straightforward, but requires a numerical inspection of the transverse kernels. Since the evo- 
lutions for ATq^^^ = h\ are purely non-singlet, diagonality in flavour of the subtraction terms 
(~ w{y\x)q{x,T)) is satisfied, while the edge-point subtractions can be tested to be positive 
numerically. We illustrate the explicit construction of the master equation for hi in LO, since 
extensions to NLO of this construction are rather straighforward. 

In this case the LO kernel is given by 



^tPII}{x) = Cf 



2 + ^5(1 -x) 



(2.6) 



and by some simple manipulations we can rewrite the corresponding evolution equation in a 
suitable master form. That this is possible is an elementary fact since the subtraction terms 
can be written as integrals of a positive function. For instance, a possibility is to choose the 
transition probabilities 

wi[x\y] = — % -2]6{y>x)9{y<l) 
y \^-x/y J 

W2[y\x] = ^ (^—^ - 0(y > -x)0(y < 0) 



(2.7) 



which reproduce the evolution equation for hi in master form 



dhi 



dywi{x\y)hi{y,T) - dyw2{y\x)hi{x,T). (2.8) 
JQ 



The NLO proof of positivity is also rather straightforward. For this purpose we have analized 
numerically the behaviour of the NLO kernels both in their bulk region and at the edge-point. 
We show in Table 1 of Appendix B results for the edge point contributions to NLO for both of 
the IS.tP'^^ components, which are numerically the same. There we have organized these terms 
in the form ~ C5{1 — x) with 

C = - \og{l- k)A + B, (2.9) 

with A and B being numerical coefficients depending on the number of flavours included in the 
kernels. The (diverging) logarithmic contribution (~ dz/{l — z)) have been regulated by 
a cutoff. This divergence in the convolution cancels when these terms are combined with the 
divergence at x = 1 of the first term of the master equation (12. 8p for all the relevant components 
containing "+" distributions. As for the bulk contributions (x < 1), positivity up to NLO of the 
transverse kernels is shown numerically in Fig. (|2.ip . All the conditions of positivity are therefore 
satisfied and therefore the 1^t±<1 distributions evolve positively up to NLO. The existence of a 
master form of the equation is then guaranteed. 

Notice that the NLO positivity of IS.T±q implies positivity of the nucleon tensor charge [32] 

5q = j dx (h\{x) - hf'''^'^ (2.10) 

for each separate flavour for positive initial conditions. As we have just shown, this proof of 
positivity is very short, as far as one can check numerically that both components of eq. (|2.8p are 
positive. 



2.3 SofFer's inequality 



Numerical tests of Soffer's inequality can be performed either in moment space or, as we are 
going to illustrate in the next section, directly in x-space, using suitable algorithms to capture 
the perturbative nature of the evolution. We recall that Soffer's inequality 

\hi{x)\<q+ix) (2.11) 

sets a bound on the transverse spin distribution hi (x) in terms of the components of the positive 
helicity component of the quarks, for a given flavour. An original proof of Soffer's inequality in 
LO has been discussed in ref.[35, while in [31] an alternative proof was presented, based on a 
kinetic interpretation of the evolution equations. 

We recall that hi, also denoted by the symbol 

Arqix, Q2) ^ q^x, Q2) _ qi(x, Q^), (2.12) 

has the property of being purely non-singlet and of appearing at leading twist. It is identifiable 
in transversely polarized hadron-hadron collisions and not in Deep Inelastic Scattering (from 
now on we will omit sometime the x-dependence in the kernels and in the distributions when 
obvious). In the following we will use interchangeably the notations hi = hi and Ayg to denote 
the transverse asymmetries. We introduce also the combinations 

Ariq + q) = hl + hl 
Arq^-^ = Ariq - q) = h\-h\ 

ATq^+') = ^AT{q^ + q^) 

i 

(2.13) 

where we sum over the flavor index (i), and we have introduced singlet and non-singlet contri- 
butions for distributions of fixed helicities 

9^ ^ = Yl (^+* + ^+*) 

i 

(-) - - V- 

q\. = q+i- q+i = S. 

(2.14) 

In our analysis we solve all the equations in the helicity basis and reconstruct the various helicities 
after separating singlet and non-singlet sectors. We mention that the non-singlet sector is now 
given by a set of 2 equations, each involving zt helicities and the singlet sector is given by a 
4- by- 4 matrix. 

In the singlet sector we have 



+P^^^«)G+ + P^^0G_). (2.15) 



while the non-singlet (valence) analogue of this equation is also easy to write down 



(-)/ 



Above, i is the flavor index, (it) indicate q:^q components and the lower subsctipt it stands for 
the helicity. 

Similarly to the unpolarized case the flavour reconstruction is done by adding two additional 
equations for each flavour in the helicity it 

= - ^^'i' (2.17) 



whose evolution is given by 



dX-i{x) _ ^S(jyNS^^, I T3NS 



dt 2tt 

(2.18) 



The reconstruction of the various contributions in flavour space for the two helicities is finally 
done using the linear combinations 

q±i = \{q''±^ +X±i + ^q'i^y (2.19) 



We will be needing these equations below when we present a proof of positivity up to LO, and 



we will thereafter proceed with a NLO implementation of these and other evolution equations. 
For this we will be needing some more notations. 

We recall that the following relations are also true to all orders 

P(x) = i(P++(x)+P+_(x)) 

between polarized and unpolarized (P) kernels and 

P++(x) = P__(x), P_+(x) = P+_(x) (2.20) 

relating unpolarized kernels to longitudinally polarized ones. Generically, the kernels of various 
type are expanded up to NLO as 

P(x) = ^pW(x) + (g)'p«(x), (2.21) 

and specifically, in the transverse case we have 



ATP^iL = ATP«±ATP4-^ (2-22) 

(2.23) 



with the corresponding evolution equations 

d 



^Tq±{Q') = ATPgq,±{as{Q')) ATq±{Q') • (2.24) 



dlnQ2 

We also recall that the kernels in the helcity basis in LO are given by 

p(0) _ p(0) _ p(0) 

^NS±,++ ~ ^qq,++ ~ ^qq 

p{0) _ p{0) _ 

^qq,+- — ^qq-+ ~ ^ 



i2 



^qg,++ — "-f-^ 

Pqg,+~ = Pqg,-+ =nf{x-l 

Pgq,++ = Pgq, 

Pgll^ = ^S!++ = ^^((r^ + ^-^-^-^')^^°'^^^-^^ 

P^l]+. = nJsx + ^-S-xA . (2.25) 



An inequality, such as Soffer's inequality, can be stated as positivity condition for suitable 
linear combinations of parton distributions [3T] and this condition can be analized - as we have 



just shown for the hi case - in a most direct way using the master form. 
For this purpose consider the Hnear valence combinations 

Q+ = Q+ + hi 

Q = q+- hi 

(2.26) 

which are termed "superdistributions" in ref.[3T]. Notice that a proof of positivity of the Q 
distributions is equivalent to verify Soffer's inequality. However, given the mixing of singlet 
and non-singlet sectors, the analysis of the master form is, in this case, more complex. As we 
have just mentioned, what can spoil the proof of positivity, in general, is the negativity of the 
kernels to higher order. We anticipate here the result that we will illustrate below where we 
show that a LO proof of the positivity of the evolution for Q can be established using kinetic 
arguments, being the kernels are positive at this order. However we find that the NLO kernels 
do not satisfy this condition. In any case, let's see how the identification of such master form 
proceeds in general. We find useful to illustrate the result using the separation between singlet 
and non-singlet sectors. In this case we introduce the combinations 



with h^{^^ = Arg(±). 

Differentiating these two linear combinations (|2.27ll we get 



(2.27) 



dlog(Q2) 

(2.28) 



which can be rewritten as 



(ilog(Q2) 2 V ++ ' T '2 

dlog(Q2) - 2 1^++"^^ +2[P+++Pt) Q- +P+-q- 



(2.29) 



with P*- ) = P^^ being the non-singlet (NS) kernel. 
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Figure 2.1: Plot of the transverse kernels. 



At this point we define the linear combinations 



P2± = \iP++^PT) (2-30) 



and rewrite the equations above as 



(2.31) 

where we have reintroduced i as a flavour index. From this form of the equations it is easy to 
establish the leading order positivity of the evolution, after checking the positivity of the kernel 
and the existence of a master form. 

The second non-singlet sector is defined via the variables 

Xi± = - ^qll^ (2.32) 
which evolve as non-singlets and the two additional distributions 

Q^i,± = Xi+ ± (2.33) 



Also in this case we introduce the kernels 




) 



(2.34) 



to obtain the evolutions 



dQxi+ 



diog(g2) 



cilog(Q2) 



(2.35) 



For the singlet sector, we simply define Q 



,(+) 



q^^\ and the corresponding evolution is 



+ 



similar to the singlet equation of the helicity basis. Using the equations above, the distributions 
Qj± are then reconstructed as 



and result positive for any flavour if the addends are positive as well. However, as we have just 
mentioned, positivity of all the kernels introduced above is easy to check numerically to LO, 
together with their diagonality in flavour which guarantees the existence of a master form. 

As an example, consider the LO evolution of Q±. The proof of positivity is a simple conse- 
quence of the structure of eq. (12. Sip . In fact the edge-point contributions appear only in 
i.e. they are diagonal in the evolution of Q±. The inhomogenous terms on the right hand side 
of (I2.3ip . proportional to q- are are harmless, since the kernel has no edge-point contri- 
butions. Therefore under 1) diagonality in flavour of the subtraction terms and 2) positivity of 
first and second term (transition probabilities) we can have positivity of the evolution. A refined 
arguments to support this claim has been presented in |29| . 

This construction is not valid to NLO. In fact, while the features of flavour diagonality of the 
master equation are satisfied, the transition probabilities w{x,y) are not positive in the whole 
x,y range. The existence of a crossing from positive to negative values in P^^ can, in fact, be 
established quite easily using a numerical analysis. We illustrate in Figs. (|2.2|) and (|2.3p plots 
of the Q kernels at LO and NLO, showing that, at NLO, the requirement of positivity of some 
components is violated. The limitations of this sort of proofs -based on kinetic arguments- are 
strictly linked to the positivity of the transition probabilities once a master form of the equation 
is identified. 




(2.36) 



^Based on the article published in JHEP 0311:059, (2003) 



2.4 The numerical investigation 



We have seen that NLO proofs of positivity, can be -at least partially- obtained only for suitable 
sets of boundary conditions. To this purpose, we choose to investigate the numerical behaviour 
of the solution using x-space based algorithms which need to be tested up to NLO. 

Our study validates a method which can be used to solve evolution equations with accuracy in 
leading and in next-to-leading order. The method is entirely based on an expansion [6j used in the 
context of spin physics [31] and in supersymmetry [H]. An interesting feature of the expansion, 
once combined with Soffer's inequality, is to generate an infinite set of relations among the scale 
invariant coefficients which characterize it. 

In this approach, the NLO expansion of the distributions in the DGLAP equation is the one 
studied in the previous section and it is given by 

where, to simplify the notation, we assume a short-hand matrix notation for all the convolution 
products. Therefore f{x,Q'^) stands for a vector having as components any of the helicities of 
the various flavours {Q±,q±,G±). The ansatz implies a tower of recursion relations once the 
running coupling is kept into account and implies that (see Eqns. 11.331) 

An+i{x) = An{x) (2.38) 

to leading order and 

= -i?n(x)+(^P(°)®A„(x)) 
ilTfJo Po 

(2.39) 

to NLO, relations which are solved with the initial condition Bq{x) = 0. The initial conditions 
for the coefficients Aq{x) and Bq{x) are specified with q{x, Qq) as a leading order ansatz for the 
initial distribution 

Ao{x) = 6{l - x) q{x, Ql) = qo{x) (2.40) 
which also requires Bq{x) = 0, since we have to satisfy the boundary condition 



Ao{x) + aoBo{x) = qo{x). 



(2.41) 



If we introduce Rossi's expansion for hi, and the linear combinations Q± (in short form) 

hi ~ (^^^r) 

(2.42) 

we easily get the inequalities 

(-1)"(^^+ + A:j) >0 (2.43) 

and 

(-1)"(^^+-^;:) >0 (2.44) 

valid to leading order,which we can check numerically. Notice that the signature factor has to 
be included due to the alternation in sign of the expansion. To next to leading order we obtain 

{-ir+'{Al+{x) + aiQ')BlHx)) <{-ir[A';,{x) +a{Q')B^{x)) < (-1)^ {Al{x) + a{Q')B?^Hx)) , 

(2.45) 

valid for n > 1, obtained after identification of the corresponding logarithmic powers log (a((5^)) 
at any Q. In general, one can assume a saturation of the inequality at the initial evolution scale 

Q_(x, Ql) = hi{x, Ql) - ^q+{x, Ql) = 0. (2.46) 

This initial condition has been evolved in Q solving the equations for the Q± distributions to 
NLO. 



2.5 Relations among moments 

In this section we elaborate on the relation between the coefficients of the recursive expansion 
as defined above and the standard solution of the evolution equations in the space of Mellin 
moments. We will show that the two sohitions can be related in an unobvious way. 

Of our concern here is the relation between the Mellin moments of the coefficients appearing 
in the expansion 

A{N) = C dxx^-^A{x) 
Jo 

B{N) = f dxx^-^B{x) 
Jo 

(2.47) 



and those of the distributions 



ATq±{N,Q^)= [ dxx^~^ATq±{x,Q^)). (2.48) 
^0 



For this purpose we recall that the general (non-singlet) solution to NLO for the latter moments 
is given by 

/ry (n^\\-^^TP!,?(N)/l3o 

ATq±{N,Q^) = K{QlQ\N)l^^] ATq±{N,Ql) 
with the input distributions Axq^iQ^) at the input scale Qq and where we have set 



ArP^lliN) - ^^ArP^^iiN) 



(2.49) 



In the expressions above we have introduced the corresponding moments for the LO and NLO 
kernels {AtP^T , ^TP^l^f)- 

We can easily get the relation between the moments of the coefficients of the non-singlet 
2;-space expansion and those of the parton distributions at any Q, as expressed by eq. (|2.49p 

An{N) + asBn{N) = ATq±{N, Ql)K{Qo, Q, N) (^Z^^^i^^^ " . (2.50) 

As a check of this expression, notice that the initial condition is easily obtained from (|2.50|) 
setting Q — > Qo,n — > 0, thereby obtaining 

<^(iV) + asB^'iN) = ATq±iN, Ql) (2.51) 

which can be solved with A^'^{N) = Arg±(A^,Qo) and B^'^{N) = 0. 

It is then evident that the expansion (12.37P involves a resummation of the logarithmic con- 
tributions, as shown in eq. ()2.5Q|) . 



2.6 An Example: The Evolution of the Transverse Spin Distri- 
butions 

LO and NLO recursion relations for the coefficients of the expansion can be worked out quite 
easily, although the numerical implementation of these equations is far from being obvious. 
Things are somehow simpler to illustrate in the case of simple non-singlet evolutions, such as 
those involving transverse spin distributions, as we are going to show below. Some details and 
definitions can be found in the appendix. 
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Figure 2.2: Plot of the LO kernels for the Q distributions 



For the first recursion relation (eq. (|2.38|1 ) we have 



^ dy yA^{y)-xAi{a 

y y -X 



+ At{x) log(l-x) 



(2.52) 



As we move to NLO, it is convenient to summarize the structure of the transverse kernel 
/^TP^q^^\x) as 



ArPg^'(^)(x) = Kf{x)5{l -x) + Kf{x)S2{x) + Kf{x) log(x) 



+Kf{x) log2(x) + Kf{x) log(x) log(l - x) + Kf{x) 



(1 - x)_ 



+ Kf{x). (2.53) 



Hence, for the (+) case we have 



AtP+'(i)(x) AUx) = K+A+ix) + [Kt{z)S2{z) + K+{z) \og{z) 

Jx y 

+ log^{z)K+{z) + log(z) log(l - z)K^{z)] A+{y) + 



1 dyyA+{y)-xA+ix) 



y -X 



+ A+(x)log(l-x)|+i^+ / ^At{y), (2.54) 

J J X y 



where z = x/y. For the (— ) case we get a similar expression. 

Now we are ready to write down the expression for the B^_^^{x) coefficient to NLO, similarly 
to eq. (|2.39|1 . So we get (for the (+) case) 



B^^,{x) = -Bt{x) + ^^\2Cp 



1 dy yA+{y) - xA+{x) 



X y 



y-x 



+ A+{x)\og{l - x) 



+ 



-2Cf ^At\y)^ + - ^K+AUx) + ^ [K^{z)S2iz)+ 



+ K+{z) log(z) + log\z)K+{z) + log(z) log(l - z)K+iz)] 



1 



A^iy) 



^ dy yA-^iy) - xA+{x) 



Cf 



X y y-x 
dy yBtiy) - xB^{x) 



+ yl+(a;)log(l -x) 



47r/3o 

Jx y 



y 



y-x 



+ B^(x)log(l-x) 



+ 



x) . 



(2.55) 



As we have already mentioned, the implementation of these recursion relations require particular 
numerical care, since, as n increases, numerical instabilities tend to add up unless high accuracy 
is used in the computation of the integrals. In particular we use finite element expansions to 
extract analitically the logarithms in the convolution (see the discussion in Appendix A). NLO 
plots of the coefficients An{x) + a{Q'^)Bn{x) are shown in figs. (|2.4|2.5p for a specific set of 
initial conditions (GRSV, as discussed below). As the index n increases, the number of nodes 
also increases. A stable implementation can be reached for several thousands of grid-points and 
up to n 10. Notice that the asymptotic expansion is suppressed by n\ and that additional 
contributions (n > 10) are insignificant even at large (> 200 GeV) final evolution scales Q. 



2.7 Nonforward Extensions 



In this section we finally discuss the nonforward extension of the evolution algorithm. In the case 
of nonforward distributions a second scaling parameter C, controls the asymmetry between the 
initial and the final nucleon momentum in the deeply virtual limit of nucleon Compton scattering. 
The solution of the evolution equations, in this case, are known in operatorial form. Single and 
double parton distributions are obtained sandwiching the operatorial solution with 4 possible 
types of initial/final states < p\...\p >, < p\...\0 >, < p'\...\p >, corresponding, respectively, to the 
case of diagonal parton distributions, distribution amplitudes and, in the latter case, skewed and 
double parton distributions [28]. Here we will simply analize the non-singlet case and discuss the 
extension of the forward algorithm to this more general case. Therefore, given the off-forward 
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Figure 2.3: Plot of the NLO kernels for the Q distributions, showing a negative behaviour at 
large x 

distributions Hq{x,S), in Ji's notation, we set up the expansion 

which is the natural extension of the forward algorithm discussed in the previous sections. We 
recall that in the light-cone gauge H{x,^) is defined as 

H,{x,iA')) = ^e--^^^-(P'|V;,(0,^,Ox)^7+^,(0,^,Ox)|P) (2.57) 

with A = P' - P, P+ = 1/2{P + P) [2^ (symmetric choice) and = 1/2 A+. 

This distribution describes for a; > ^ and x < — ^ the DGLAP-type region for the quark 
and the antiquark distribution respectively, and the ERBL [33] (see also [BB] for an overview) 
distribution amplitude for — ^ < x < ^. In the following we will omit the A dependence from 

Again, once we insert the ansatz (I2.56P into the evolution equations we obtain an infinite set 
of recursion relations which we can solve numerically. In LO, it is rather simple to relate the 
Gegenbauer moments of the skewed distributions and those of the generalized scaling coefficients 
An. We recall that in the nonforward evolution, the multiplicatively renormalizable operators 
appearing in the light cone expansion are given in terms of Gegenbauer polynomials [28] . The 
Gegenbauer moments of the coefficients An of our expansion (|2.56p can be easily related to those 
of the off-forward distribution 




Cn{^, Q') = C Cl'\z/i)H{z, e, Q'')dz. (2.58) 
The evolution of these moments is rather simple 



with 




(2.60) 



being the non-singlet anomalous dimensions. If we define the Gegenbauer moments of our ex- 
pansion 



A^^\i,Q^) = C fcl'\z/OH{zi,Q^)dz (2.61) 



we can relate the moments of the two expansions as 



k 



Ar{i) = Cn{C,Qt)[f^j ■ (2.62) 

Notice that expansions similar to (|2.56|) hold also for other choices of kinematical variables, such 
as those defining the non- forward distributions [28], where the t-channel longitudinal momentum 
exchange A+ is related to the longitudinal momentum of the incoming nucleon as A = C,P. We 
recall that Hg{x.(,) as defined in [27] can be mapped into two independent distributions Tq{X, Q 
and J^q{X, Q through the mappings [37] 

(1 + 



(1+0 

^ = C/(2-C) 
1 

1-7/2^ 
1 



(2.63) 



in which the interval — 1 < x < 1 is split into two coverings, partially overlapping (for — ^ < a; < ^, 
or ERBL region) in terms of the two variables — ^ < xi < 1 (0 < Xi < 1) and — 1 < X2 < 
(0 < X2 < 1). In this new parameterization, the momentum fraction carried by the emitted 
quark is X, as in the case of ordinary distributions, where it is parametrized by Bjorken x. For 
definitess, we focus here on the DGLAP-like (X > C) region of the non-singlet evolution. The 



non-singlet kernel is given in this case by (x = X) 



a 



PAxX) = -Cf 



vr 



1 



y-x 



1 + 



XX 

yy' 



6{x — y) dz 



1 + 
1 - z 



(2.64) 



we introduce a LO ansatz 



k=0 



(2.65) 



and insert it into the evolution of this region to obtain the very simple recursion relations 



/3o 
2 

To 



Cf 



UyyAn{yX)-xAn{XX) 2^ fUy{X - Q {yAn{X,C) - XAn{y,0) 



X y 

CFAn{X,0 



y-X 



Po^^ Jx yiy-C) 



y-x 



3 ^^ {l-X)\l-x/C) 
2 1 - C 



(2.66) 



The recursion relations can be easily reduced to a weighted sum of contributions in which C is a 
spectator parameter. Here we will not make a complete implementation, but we will illustrate 
in an appendix the general strategy to be followed. There we show a very accurate analytical 
method to evaluate the logarithms generated by the expansion without having to rely on brute- 
force computations. 



2.8 Positivity of the non-singlet Evolution 



Positivity of the non-singlet evolution is a simple consequence of the master-form associated to 
the non- forward kernel (|2.64p . As we have already emphasized above, positivity of the initial 
conditions are sufficient to guarantee a positivity of the solution at any scale Q. The master- 
form of the equation allows to reinterpret the parton dynamics as a random walk biased toward 
small-x values as r = log(Q^) increases. 

In the non-forward case the identification of a transition probability for the random walk 
[29] associated with the evolution of the parton distribution is obtained via the non-forward 
transition probability 



W(^{x\y) 
w'^{y\x) 



^Cf- 
vr y 



1 



X 



1 + ^ ~ 

y y-C 



IT — y) 



and the corresponding master equation is given by 



Oiy > x) 



(2.67) 



dJ^a 



dyw(^{x\y)Tq{y,Q, 



dy w[{y\x)J^q{x,C, 



(2.68) 



dT Jx Jo 
that can be re-expressed in a form which is a simple generalization of the formula for the forward 



evolution 



dJ^q 



nl px 

dlogQ2 ^ J dyw(;{x\y)Tg{yX,T-) - dy w'^{y\x)Tg{x, C,t) 

j-a(x) p~x 

= - dywc_{x + y\x)*Tq{x,C,T) + I + y|x)J^g(x, C, r), (2.69) 

Jo Jo 

where a Moyal-like product appears 

Wi;{x + y\x)*J='q{x,C,'r) = Wi;{x + y\x)e J='q{x,C,T) (2.70) 

and a{x) = x — 1. A Kramers-Moyal expansion of the equation allows to generate a differential 
equation of infinite order with a parametric dependence on ( 



dTq 



Ci log Q2 



t-—x 

dywc_{x + y\x)Tq{x,C,T) + / dy w'(-{x + y\x)Tq{xX^T) 

cx{x) Jo 

dy^—^d,^{w^ix + y\x)J^q{xX,r)). (2.71) 



n=l 



We define 

/"O r—x 

aoix,C) = / dywQ{x + y\x)J='q{x,C,r) + dyw'^{x + y\x)J^q{x,C,T) 

Ja(x) Jo 
j-a{x) 

an{x,C) = / dyy''w(;{x + y\x)J='q{x,C,r) 
Jo 

a„(x,C) = / dyy''d^''{w(^{x + y\x)Fq{xX,T)) n=l,2,... (2.72) 
Jo 

If we arrest the expansion at the first two terms (n = 1, 2) we are able to derive an approximate 
equation describing the dynamics of partons for non-diagonal transitions. The procedure is a 
slight generalization of the method presented in |29| . to which we refer for further details. For 
this purpose we use the identities 



ai(x,C) = dxai{xX) — a{x)dxa{x)w(^{x + a{x)\x)!Fq{xX^T) 
02(2;, C) = d1a2{x,C,) -'^oi{x){dxa{x))'^WQ{x + a{x)\x)J^q{x,C,,T) 
-a{xfdxa{x)dx {w^^ix + a{x)\x)J^q{x,C,T)) 

-a^{x)dxa{x)dx {w^{x + y\x)J^q{xX,T)) \y=a{x)- (2-73) 



which allow to compute the first few coefficients of the expansion. Using these relations, the 
Fokker-Planck approximation to this equation can be worked out explicitely. We omit details 



on the derivation which is unobvious since particular care is needed to regulate the (canceling) 
divergences and just quote the result. 

A lengthy computation gives 



dr 



Ot„ ( Xo_3 Xo_i \ 

a ( Xi_3 Xi-i \ ^ . N , ^0-3 p,2-t: I a \ 



(2.74) 



where we have defined 

- ((-1 + xf (17x3 - (3 + 40 + SxC (3 + 50 - 3x2 (3 + 7^))^ 

"^"■-^ = 12^^ 

-29x^ - 3 + x^ (-1 + + 2C - 2x (1 + 30 + (12 + 230 
= 3^5 

1 3 , (1-x) 
a;o,o = 4+-^-- + 21og^ 

ZX^ X X 



Xl _l 



((-1 + 6x - 15x2 + 14x3) (x - C)) 
3x2 



1 5x , 23x^ 7C 3C 5xC 



2, , 131x3C 5(2 , e , ,0 .2 39x2(2 ^3 g^^a 



-ISx^C + — — ^ - ^ + ^ - — + 13^C 7^ - 3C + + 



12 2 4x2 X ^ 4 32;2 3 



((-1 + x)2 (x - C)^ (3 + 23x2 + 4C - 2x (7 + 80)) 



X2,-3 = ^ (2.75) 

This equation and all the equations obtained by arresting the Kramers-Moyal expansion to 
higher order provide a complementary description of the non-forward dynamics in the DGLAP 
region, at least in the non-singlet case. Moving to higher order is straightforward although the 
results are slightly lengthier. A full-fledged study of these equations is under way and we expect 
that the DGLAP dynamics is reobtained - directly from these equations - as the order of the 
approximation increases. 



2.9 Model Comparisons, Saturation and the Tensor Charge 

In this last section we discuss some implementations of our methods to the standard (forward) 
evolution by doing a NLO model comparisons both in the analysis of Sofl'^er's inequality and for 
the evolution of the tensor charge. We have selected two models, motivated quite independently 
and we have compared the predicted evolution of the SofFer bound at an accessable flnal evolution 
scale around 100 GeV for the light quarks and around 200 GeV for the heavier generations. At 
this point we recall that in order to generate suitable initial conditions for the analysis of Soffer's 
inequality, one needs an ansatz in order to quantify the difference between its left-hand side and 



right-hand side at its initial value. 

The well known strategy to build reasonable initial conditions for the transverse spin distri- 
bution consists in generating polarized distributions (starting from the unpolarized ones) and 
then saturate the inequality at some lowest scale, which is the approach we have followed for all 
the models that we have implemented. 

Following Ref. [3^ (GRSV model), we have used as input distributions - in the unpolarized 
case - the formulas in Ref. [39], calculated to NLO in the MS scheme at a scale Qq = 0.40 GeV^ 

x{u-u){x,Ql) = 0.632x°-^^(l-x)2-°^(l + 18.2x) 
x{d-d){x,Ql) = 0.624(1 -x)^-°x(n-n)(x,Qg) 
x(d-u){x,Ql) = 0.20x°-*^3(l-3;)i2.4(^_-^3 3^_^gQQ^^ 

x{u + d){x,Ql) = 1.24x°-2°(l-a;)^-^(l-2.3Vx + 5.7x) 

xgix,Ql) = 20.80x^-^(1 - x)''-^ (2.76) 

and xqi{x, Qq) = xql{x, Ql) = for qi = s, c, b, t. 

We have then related the unpolarized input distribution to the longitudinally polarized ones 
by as in Ref. [38] 

xAu{x,Ql) = im^x^-^'^{l-xf-^'^xu{x,Ql) 

xAd{x,Ql) = -0.669x°-'^^xd(x,Qg) 

xAu{x,Ql) = -0.272x°-^^xIi(x,Q^) 

xAd(x,Qo) = xAu{x,Qq) 

xAg{x,Ql) = lA19x^-'^^{l-xf-^^xg{x,Ql) (2.77) 

and xAqi{x, Qq) = xAgi(x, Qq) = for qi = s, c, b, t. 

Following [10], we assume the saturation of Soffer inequality: 

xA^,.(x, Ql) = (2.78) 

and study the impact of the different evolutions on both sides of Soffer's inequality at various 
final evolution scales Q. 

In the implementation of the second model (GGR model) we have used as input distributions 
in the unpolarized case the CTEQ4 parametrization [41], calculated to NLO in the MS scheme 
at a scale Qq = 1.0 GeV 



x{u-u){x,Ql) = 1.344x°-50i(l-x)3-689(l + 6.402x°-^^3) 

x{d-d){x,Ql) = 0.64x°-^°i(l-x)^-247(i + 2.69x0-333) 

xs{x, Ql) = xs{x, Ql) = 0.064j;-0-^^3(i _ a;)8.04i(^ ^ 6.112x) 

x{d-u){x,Ql) = 0.071x°-50^(l-x)S-°^^(l + 30.0x) 

x{u + d){x,Ql) = 0.255x-0-i^3(i _^)8.04i(^^g_-^^2x) 

xg{x,Ql) = 1.123x-0-206(i_a;)4.673(;L^4 269a;i-508) (2.79) 

and xqi{x,Qf^) = xgi(x,(5o) = fo'^ = c, 6, t and we have related the unpolarized input 
distribution to the longitudinally polarized ones by the relations |l2] 



xAm(x, Qo) 
xAn(x, (5o) 
xAd(x,Qg) 
xAd(x,Qg) 
xAs(x, Qo) = xAs(x, Qq) 



xriu{x)xu{x,Ql) 
cos ^D(a;, Qo 

cos6d{x,Q1) - 
xrjs{x)xs{x,Ql) 



x{u — u) — -x{d — d) 



{x, Ql) + xAu(x, Qo) 



-x{d - (i)(x,Qo) 



+ xAd{x,Ql) 



(2.80) 



and xAqi{x, Qq) = xAqi{x, Qq) = for qi = c, b, t. 

A so-called "spin dilution factor" as defined in [12], which appears in the equations above is 
given by 

2a,(Q2) (i-x)2n-i 



cos 6d{x,Qq) 



1 + 



(2.81) 



In this second (GGR) model, in regard to the initial conditions for the gluons, we have made use 
of two different options, characterized by a parameter 77 dependent on the corresponding option. 
The first option, that we will denote by GGRl, assumes that gluons are moderately polarized 



xAg{x,Ql) = x-xg(x,Ql) 
Vu{x) = Vdix) = -2.49 + 2.8^/^ 
7]s{x) = -1.67 + 2.1^^, 

while the second option (GGR2) assumes that gluons are not polarized 

: 



(2.82) 



xAg{x,Ql) 
r]u{x) = r]d{x) 

7]s{x) 



-3.03 + 3.0^/^ 
-2.71 2.9^/^. 



(2.83) 



We have plotted both ratios At/ and differences {xf~^ — xAxf) for various flavours as a func- 
tion of X. For the up quark, while the two models GGRl and GGR2 are practically overlapping, 
the difference between the GGR and the GRSV models in the the ratio Atu/u~^ is only slightly 
remarked in the intermediate x region (0.1 — 0.5). In any case, it is just at the few percent level 
(Fig. (|2.6|) ). while the inequality is satisfied with a ratio between the plus helicity distribution 
and transverse around 10 percent from the saturation value, and above. There is a wider gap in 
the inequality at small x, region characterized by larger transverse distribution, with values up 
to 40 percent from saturation. A similar trend is noticed for the x-behaviour of the inequality in 
the case of the down quark (Fig. 12. 7|) . In this latter case the GGR and the GRSV model show a 
more remarked difference, especially for intermediate x-values. An interesting features appears 
in the corresponding plot for the strange quark (Fig. (12. 8]) ). showing a much wider gap in the 
inequality (50 percent and higher) compared to the other quarks. Here we have plotted results 
for the two GGR models (GGRl and GGR2). Differently from the case of the other quarks, in 
this case we observe a wider gap between Ihs and rhs at larger x values, increasing as x — > 1. In 
figs. (12. 91) and (|2.10l) we plot the differences (x/"*" —xAt/) for strange and charm and for bottom 
and top quarks respectively, which show a much more reduced evolution from the saturation 
value up to the final corresponding evolving scales (100 and 200 GeV). As a final application 
we finally discuss the behaviour of the tensor charge of the up quark for the two models as a 
function of the final evolution scale Q. We recall that like the isoscalar and the isovector axial 
vector charges defined from the forward matrix element of the nucleon, the nucleon tensor charge 
is defined from the matrix element of the tensor current 

(P^tIVS^t'^'tsA^IP, St) = 25q''{Ql) {P^S!^ - P'^S!^) (2.84) 

where (5"g(Qo) denotes the flavour (a) contribution to the nucleon tensor charge at a scale Qo 
and St is the transverse spin. 

In fig. (12. lip we plot the evolution of the tensor charge for the models we have taken in 
exam. At the lowest evolution scales the charge is, in these models, above 1 and decreases 
slightly as the factorization scale Q increases. We have performed an evolution up to 200 GeV as 
an illustration of this behaviour. There are substantial differences between these models, as one 
can easily observe, which are around 20 percent. From the analysis of these differences at various 
factorization scales we can connect low energy dynamics to observables at higher energy, thereby 
distinguishing between the various models. Inclusion of the correct evolution, up to subleading 
order is, in general, essential. 

2.10 Conclusions 

We have illustrated the use of x-space based algorithms for the solution of evolution equations in 
the leading and in the next-to-leading approximation and we have provided some applications of 
the method both in the analysis of Soffer's inequality and in the investigation of other relations. 



such as the evolution of the proton tensor charge, for various models. The evolution has been 
implemented using a suitable base, relevant for an analysis of positivity in LO, using kinetic 
arguments. The same kinetic argument has been used to prove the positivity of the evolution 
of hi and of the tensor charge up to NLO. In our implementations we have completely relied 
on recursion relations without any reference to Mellin moments. We have provided several 
illustrations of the recursive algorithm and extended it to the non-forward evolution up to NLO. 
Building on previous work for the forward evolution, we have presented a master-form of the 
non-singlet evolution of the skewed distributions, a simple proof of positivity and a related 
Kramers Moyal expansion, valid in the DGLAP region of the skewed evolution for any value 
of the asymmetry parameter (^. We hope to return with a complete study of the nonforward 
evolution and related issues not discussed here in the near future. 



2.11 Appendix A. Weighted Sums 

In this appendix we briefly illustrate the reduction of recursion relations to analytic expressions 
based on finite element decompositions of the corresponding integrals. The method allows to 
write in analytic forms the most dangerous integrals thereby eliminating possible sources of insta- 
bilities in the implementation of the recursion relations. The method uses a linear interpolation 
formula for the coefficients An, which, in principle can also be extended to higher (quadratic) 
order. However, enugh accuracy can be achieved by increasing the grid points in the discretiza- 
tion. Notice that using this method we can reach any accuracy since we have closed formulas for 
the integrals. In practice these and similar equations are introduced analytically as functions in 
the numerical integration procedures. 

Below, we will use a simplified notation {X = x for simplicity). 

We define P{x, Q = xP{x, Q and A{x, Q = xA{x) and the convolution products 



The integration interval in y at any fixed x- value is partitioned in an array of increasing points 
ordered from left to right {xo,xi,X2, ■■■,Xn,Xn+i) with xq = x and x^+i = 1 being the upper 
edge of the integration region. One constructs a rescaled array {x,x/xn, ■■■,x/x2,x/xi,l). We 
define Si = x/xi, and s„+i = x < Sn < s„_i < ...si < sq = 1. We get 




(2.85) 




(2.86) 



At this point we introduce the linear interpolation 




(2.87) 



and perform the integration on each subinterval with a change of variable y— > x/y and replace 
the integral J{x,Q with its discrete approximation Jn{x) to get 

Mx,Q = A{xo)-^ f'^P{y,Q{y-s,) 

- •si Jsi y 



^ /-s -1 rl 

*i J Si y 



(2.88) 



with the condition A{xjsi+i,Q = 0- Introducing the coefficients W(x,x,Q and W{xi,x,C,), the 
integral is cast in the form 



Jn{x, C) = W{x, X, C)A(x, O + Y. ^(^i' C)^(^i, C) 



(2.89) 



where 



1 - Si J SI y 



W{xi,x,C) = 



Si - Si+i Jg.^^ y 

Si p-' dy 
Si-i - Si Js y 



^{y-Si+,)P{y,0 

(y-Si_i)P(y,C). 



(2.90) 



For instance, after some manipulations we get 



I dy yAniy,C) xAniX,0 ^ ino(^)^„(^, iJn,{x) - JnU{x)) An{xi)-ln{l-x)An{x,0 

Jx y y — Tr^ 

(2.91) 



where 
lo{x) 

i',{x) 

Jl(x) 



1 



1 

Si—\ Si 
1 



log(si) + log(l - Si) 



1 - Si 



log ( T- — — 1 + Si-i log 



1-s. 



i-1 



1 - Sj+1 Si 

Si 



Si-1 



Si-1 



1 - s. 



i-1 



1 - Si 



1-si 



log Si. 



i = 2,3,. .N 



(2.92) 



These functions, as shown here, and similar ones, are computed once and for all the kernels and 
allow to obtain very fast and extremely accurate implementations for any ^. 
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Figure 2.4: Coefficients An{x) + as{Q^)Bn, with n = 0, . . . , 4 for a final scale Q = 100 GeV for 
the quark up. 
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Uf 


A 


B 


3 


12.5302 


12.1739 


4 


10.9569 


10.6924 


5 


9.3836 


9.2109 


6 


7.8103 


7.7249 



Table 1. Coefficients A and B for various flavour, to NLO for ATPqq,±. 



2.13 Appendix C 



Here we define some notations in regard to the recursion relations used for the NLO evolution 
of the transverse spin distributions. 
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Figure 2.5: Coefficients An{x) + as{Q^)Bn, with n = 0, . . . , 4 for a final scale Q = 100 GeV for 
the quark down. 



For the (+) case we have these expressions 

K+{x) = ^CF(-2n/(3 + 47r2) + iVc(51 + 447r^ - 216C(3)) + 9Cf(3 - iir^ + 48C(3)) 



^(^)^ ^^-^-^^-;^^- (2.93) 
Kti.) = ^H9C.-lliVc + 2n,). 



Kti^) = ^ (2.95) 

= ^ (2.96) 

K+{x) = -^CF(10n/ + Arc(-67 + Sn^)) (2.97) 

K+{x) = ^CF(10n/ + Nc{-67 + Sir^)) (2.98) 



(2.99) 




Figure 2.6: Test of Soffer's inequality for quark up at Q = 100 GeV for different models. 



and for the (— ) case we have 
1 



CF{-2nf{3 + 47r^) + Nc{51 + A4tt^ - 216C(3)) + 9Cf(3 - Att^ + 48C(3)) 



72 

2Cf{+2Cf - Nc)x 
1 + x 

CF{9CF-llNc + 2nf)x 
3(x - 1) 

CfNcx 
1 — X 

1-x 

—CFilOuf + Arc(-67 + Stt^)) 
9 

-^CF(10n/ - 18Cf(x - 1) + Nc{-76 + Stt^ + 9x)) 



(2.100) 

(2.101) 

(2.102) 

(2.103) 

(2.104) 

(2.105) 
(2.106) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 

Figure 2.7: Test of SofFer's inequality for quark down at Q = 100 GeV for different models 
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Figure 2.9: Soffer's inequality for strange and charm in the GRSV model. 
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Figure 2.10: Soffer's inequality for bottom and top in the GRSV model. 
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Figure 2.11: Tensor charge qt as a function of Q for up and down quark for the GRSV and GGR 
models. 



Chapter 3 



Double Transverse- Spin Asymmetries in 
Drell-Yan Processes with Antiprotons 



3.1 Introduction 

In this chapter we will use what we learned about the evolution of the transversely polarized 
parton distributions to make predictions about the NLO asymmetry in DY processes with trans- 
versely polarized antiprotons. 

The experiments with antiproton beams planned for the next decade in the High- Energy Stor- 
age Ring at GSI will provide a variety of perturbative and non-perturbative tests of QCD |44j. 
In particular, the possible availability of transversely polarised antiprotons opens the way to 
direct investigation of transversity, which is currently one of the main goals of high-energy spin 
physics [45]. The quark transversity [i.e. transverse polarisation) distributions AtQ were first 
introduced and studied in the context of transversely polarised Drell-Yan (DY) production [46j; 
this is indeed the cleanest process probing these quantities. In fact, whereas in semi-inclusive 
deep-inelastic scattering transversity couples to another unknown quantity, the Collins frag- 
mentation function [47], rendering the extraction of AtQ a not straightforward task, the DY 
double-spin asymmetry 

^TT - rf^TT _ dan - ^ y-^-^^ 

only contains combinations of transversity distributions. At leading order, for instance, for the 
process p^p^ — > I'^l^ X one has 

'''^^ Egel[q{xi,M^)q{x2,M^) + q{xi,M^)q{x2,M^)] ' ^ ' > 

where M is the invariant mass of the lepton pair, q{x,M'^) is the unpolarised distribution func- 
tion, and aTT is the spin asymmetry of the QED elementary process qq £~^i~ . In the dilepton 



^ Based on the article [l^ 
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centre-of-mass frame, integrating over the production angle 0, one has 



(3.3) 



where (p is the angle between the dilepton direction and the plane defined by the collision and 
polarisation axes. 

Measurement oip^p^ DY is planned at RHIC |l8j. It turns out, however, that A^^{pp) is 
rather small at such energies [lH IMl El], no more than a few percent (similar values are found 
for double transverse-spin asymmetries in prompt-photon production [52| and single-inclusive 
hadron production [53]). The reason is twofold: 1) Af^^{pp) depends on antiquark transversity 
distributions, which are most likely to be smaller than valence transversity distributions; 2) 
RHIC kinematics {^/^ = 200 GeV, M < lOGeV and xiX2 = M'^/s < 3 x lO'^) probes the 
low-x region, where QCD evolution suppresses Axqix, IVP) as compared to the unpolarised 
distribution q{x,M'^) [Ml ES]. The problem may be circumvented by studying transversely 
polarised proton-aniiproton DY production at more moderate energies. In this case a much 
larger asymmetry is expected [MUMIET] since A^^{pp) is dominated by valence distributions at 
medium x. The PAX collaboration has proposed the study of p^p^ Drell-Yan production in the 
High-Energy Storage Ring (HESR) at GSI, in the kinematic region 30GeV2 < s < 200 GeV^ 
2 GeV < M < 10 GeV and xiX2 > 0.1 [58]. An antiproton polariser for the PAX experiment is 
currently under study [59]: the aim is to achieve a polarisation of 30-40%, which would render 
the measurement of A^^{pp) very promising. 

Leading-order predictions for the pp asymmetry at moderate s were presented in [56]. It was 
also suggested there to access transversity in the J/ip resonance production region, where the 
production rate is much higher. The purpose of this study is to extend the calculations of [56] 



to next-to-leading order (NLO) in QCDB This is a necessary check of the previous conclusions, 
given the moderate values of s in which we are interested. We shall see that the NLO corrections 
are actually rather small and double transverse-spin asymmetries are confirmed to be of order 
20-40%. 

3.2 The kinematics 

The kinematic variables describing the Drell-Yan process are (1 and 2 denote the colliding 
hadrons) : 



with r = M'^/s. We denote by xi and X2 the longitudinal momentum fractions of the incident 
partons. At leading order, and ^2 coincide with xi and X2, respectively. The QCD factorisation 
formula for the transversely polarised cross-section for the proton-antiproton Drell-Yan process 

^The results presented here were communicated at the QCD-PAC meeting at GSI (March 2005) and reported 
by one of us (M.G.) at the Int. Workshop "Transversity 2005" (Como, September 2005) ifiOj. 




6 = V^ey 



6 = V^e-y 




(3.4) 



IS 



+ A,,-(x„;.^)A,,-(x„;.^)]^^^, (3.5) 



where /j. is the factorisation scale and we take the quark (antiquark) distributions of the 
antiproton equal to the antiquark (quark) distributions of the proton. Note that, since gluons 
cannot be transversely polarised (there is no such thing as a gluon transversity distribution for 
a spin one-half object like the proton), only quark- antiquark annihilation subprocesses (with 
their radiative corrections) contribute to cZAtct. In Eq. (|3.5I) we use the fact that antiquark 
distributions in antiprotons equal quark distributions in protons, and viceversa. At NLO, i.e. at 
order as, the hard-scattering cross-section dAyo"^^^, taking the diagrams of Fig. [SH] into account, 
is given by [50] 



003(299) 
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,,,3,1, (l-6)(l-6) ^ , '^(xi-gi) 
5(xi - Ci)6{x2 - 6) ( 4 + 2 In ^ J + (^^3^ 

+ [1^2], 



(3.6) 



where we have taken the factorisation scale fi equal to the renormalisation scale. In our calcula- 
tions we set // = M. 



(a) (b) (c) 

Figure 3.1: Elementary processes contributing to the transverse Drell-Yan cross-section at NLO: 
(a, b) virtual-gluon corrections and (c) real-gluon emission. 



The unpolarised Drell-Yan differential cross-section can be found, for instance, in [6T]; besides 
the diagrams of Fig. 13.11 it also includes the contribution of quark-gluon scattering processes. 



3.3 Drell-Yan Asymmetries 



To compute the Drell-Yan asymmetries we need an assumption for the transversity distributions, 
which as yet are completely unknown. We might suppose, for instance, that transversity equals 
helicity at some low scale, as suggested by confinement models [Hj (this is exactly true in the 
non-relativistic limit). Thus, one possibility is 

ATg(x,/x2) = A(7(x,^2)^ (37) 

where typically fJ-o ^ 1 GeV. Another possible assumption for AtQ is the saturation of Soffer's 
inequality [i^2j, namely 

|ATg(x,/xo)| = i[g(x,/io) + Ag(x,/io)], (3.8) 
which represents an upper bound on the transversity distributions. 

Since Eqs. (|3.7I) and (|3.8p make sense only at very low scales, in practical calculations one 
has to resort to radiatively generated helicity and number densities, such as those provided by 
the GRV fits The GRV starting scale is indeed (at NLO) = 0.40 GeV^. We should 

however bear in mind that in the GRV parametrisation there is a sizeable gluon contribution 
to the nucleon's helicity already at the input scale {Ag is of order 0.5). On the other hand, 
as already mentioned, gluons do not contribute to the nucleon's transversity. Thus, use of 
Eq. (13.7(1 with the GRV parametrisation may lead to substantially underestimating the quark 
transversity distributions and hence is a sort of "minimal bound" for transversity. Incidentally, 
the experimental verification or otherwise of the theoretical predictions of Att based on the 
low-scale constraints (13.71 13.81) would represent an indirect test of the "valence glue" hypothesis 
behind the GRV fits. Note too that, although the assumption (|3.7I) may, in principle, violate the 
Soffer inequality, we have explicitly checked that this is not the case with all the distributions 
we use. 

After setting the initial condition (|3.7p or (|3.8p . all distributions are evolved at NLO according 
to the appropriate DGLAP equations (for transversity, see [6^; the numerical codes we use to 
solve the DGLAP equations are described in [3j). The u sector of transversity is displayed in 
Fig. 13.21 for the minimal bound (|3.7p and for the Soffer bound (|3.8p . 

The transverse Drell-Yan asymmetry A^^ /axT, integrated over M between 2 GeV and 3 GeV 
{i.e. below the J ftp resonance region), for various values of s is shown in Fig. 13.31 As can be 
seen, the asymmetry is of order of 30% for s = 30 GeV^ (fixed-target option) and decreases 
by a factor two for a centre-of-mass energy typical of the colhder mode {s = 200 GeV^). The 
corresponding asymmetry obtained by saturating the Soffer bound, that is by using Eq. (|3.8|) for 
the input distributions, is displayed in Fig. 13.41 As expected, it is systematically larger, rising 
to over 50% for fixed-target kinematics. 

Above the J/ip peak Ai^^ /utt appears as shown in Fig. 13.51 where we present the results 
obtained with the minimal bound (|3.7p . Comparing Figs. l3.3l and l3.5l we see that the asymmetry 
increases at larger M (recall though that the cross-section falls rapidly with growing M). 



The importance of NLO QCD corrections may be appreciated from Fig. 13.61 where one 
sees that the NLO effects hardly modify the asymmetry since the K factors of the transversely 
polarised and unpolarised cross-sections are similar to each other and therefore cancel out in the 
ratio. As for the dependence on the factorisation scale fi (we recall that the results presented in 
all figures are obtained setting ^ = M), we have repeated the calculations with two other choices 
{fj, = 2M and /j. = M/2) and found no sensible differences. 

A caveat is in order at this point. The GSI kinematics is dominated by the domain of large 
r and large z = t/xiX2, where real-gluon emission is suppressed and where there are powers of 
large logarithms of the form ln(l — z), which need to be resummed to all orders in [65]. It 
turns out that the effects of threshold resummation on the asymmetry A^^ in the regime we 
are considering, although not irrelevant, are rather small (about 10%) if somewhat dependent 
on the infrared cutoff for soft-gluon emission. 

The feasibility of the Att measurement at GSI has been thoroughly investigated by the PAX 
Collaboration (see App. F of [58]). In collider mode, with a luminosity of 5 • lO^'^ cm~^ s~^, a 
proton polarisation of 80%, an antiproton polarisation of 30% and considering dimuon invariant 
masses down to M = 2 GeV, after one year's data taking one expects a few hundred events per 
day and a statistical accuracy on Att of 10-20%. 

3.4 The J/iIj region 

Before concluding, we briefly comment on the possibility of accessing transversity via J/i/j pro- 
duction in pp scattering. It is known that the dilepton production rate around M = 3 GeV, i.e. 
at the J/ip peak, is two orders of magnitude higher than in the region M ~ 4 GeV. Thus, with 
a luminosity of 5 • 10'^'^ cm~^ s~^, one expects a number of pp J /il) £^£^ events of order 10^ 
per year at GSI collider energies. This renders the measurement of Att in the J/^/'-resonance 
region extremely advantageous from a statistical point of view. 

As explained in [56], if J/^p formation is dominated by the qq annihilation channel, at leading 
order the double transverse-spin asymmetry at the J/ip peak has the same structure as the 
asymmetry for Drell-Yan continuum production, since the J/ip is a vector particle and the 
qqJ/ip coupling has the same helicity structure as the qq^* coupling. The GERN SPS data [HE] 
show that the pp cross-section for J/ip production at s = 80 GeV^ is about ten times larger than 
the corresponding pp cross-section, which is a strong indication that the gg-fusion mechanism is 
indeed dominant. Therefore, at the s values of interest here [s < 200 GeV^) dilepton production 
in the J/ip resonance region can be described in a manner analogous to Drell-Yan continuum 
production, with the elementary subprocess — > 7* — > £~i~^ replaced hy qq ^ J/ip ^ [67]. 
Using this model, which successfully accounts for the SPS J/ip production data at moderate 
values of s, it was found in [56] that the transverse asymmetry at the J/tp peak is of the order 
of 25-30%. 

At next-to-leading order, due to QGD radiative corrections, one cannot use a point-like 



qqJ/ip coupling, and therefore it is not possible to extend in a straightforward way the model 
used to evaluate Att at leading order. Were NLO effects not dominant, as is the case for 
continuum production, one could still expect the J/i/j asymmetry to be quite sizeable, but this 
is no more than an educated guess. What we wish to emphasise, however, is the importance 
of experimentally investigating the J/^p double transverse-spin asymmetry, which can shed light 
both on the transversity content of the nucleon and on the mechanism of J/xp formation (since 
gluon-initiated hard processes do not contribute to the transversely polarised scattering, the 
study of Att in the J/ip resonance region may give information on the relative weight of gluon 
and quark-antiquark subprocesses in J ftp production). 



3.5 Conclusions 



In conclusion, experiments with polarised antiprotons at GSI will represent a unique opportunity 
to investigate the transverse polarisation structure of hadrons. The present chapter, which 
confirms the results of [56], shows that the double transverse-spin asymmetries are large enough 
to be experimentally measured and therefore represent the most promising observables to directly 
access the quark transversity distributions. 
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Figure 3.2: The u and u transversity distributions, as obtained from the GRV parametrisation 
and Eq. top panel, or Eq. (IT8]l . bottom: xAtu at //^ = = 0.40 GeV^ (dashed curve) 

and = OGeV^ (sohd curve); xAtu at /i^ = ^2 ^ 0.40 GeV^ (dotted curve) and = 9GeV^ 
(dot-dashed curve). Note that the u transversity distributions have been multiphed by a factor 
of 10. 
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Figure 3.3: The NLO double transverse-spin asymmetry ATT{y)/o-TT, integrated between M = 
2 GeV and M = 3 GeV, for various values of s; the minimal bound (|3.7p is used for the input 
distributions. 
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Figure 3.4: As Fig. 13.31 but with input distributions corresponding to the Soffer bound (|3.8I) . 
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Figure 3.5: The NLO double transverse-spin asymmetry ATT{y)/o-TT, integrated between M = 
4GeV and M = 7GeV, for various values of s; the minimal bound (|3.7p is used for the input 
distributions. 
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Figure 3.6: NLO vs. LO double transverse-spin asymmetry ATTiy)/o-TT at M = 4 GeV for 
s = 45GeV^ and s = 200 GeV^; the minimal bound (|3.7I1 is used for the input distributions. 



Chapter 4 



On the Scale Variation of the Total 
Cross Section for Higgs Production at 
the LHC and at the Tevatron 

The validity of the mechanism of mass generation in the Standard Model will be tested at the 
new collider, the LHC. For this we require precision studies in the Higgs sector to confirm its 
existence. This program involves a rather complex analysis of the QCD backgrounds with the 
corresponding radiative corrections fully taken into account. Studies of these corrections for 
specific processes have been performed by various groups, to an accuracy which has reached the 
next-to-next-to-leading order (NNLO) level in a^, the QCD coupling constant. The quantification 
of the impact of these corrections requires the determination of the hard scattering partonic 
cross sections up to order , together with the DGLAP kernels controlling the evolution of the 
parton distributions determined at the same perturbative order. Therefore, the study of the 
evolution of the parton distributions, using the three-loop results on the anomalous dimensions 
|17| . is critical for the success of this program. Originally NNLO predictions for some particular 
processes such as total cross sections [6^ have been obtained using the approximate expressions 
for these kernels |69| . The completion of the exact computation of the NNLO DGLAP kernels 
motivates more detailed studies of the same observables based on these exact kernels and the 
investigation of the factorization (fJ-F) and renormalization (fJ^R) scale dependences of the result, 
which are still missing. In this work we are going to reanalyze these issues from a broader 
perspective. Our analysis is here exemplified in the case of the total cross sections at the LHC 
(pp) and at the Tevatron (pp) for Higgs production using the hard scatterings computed in [68] 
and their dependence on the factorization and renormalization scales. Our study is based on the 
exact and well defined NNLO computations of the hard scatterings for this process and we have 
not taken into account any threshold re summation since this involves further approximations. 
These effects have been considered in [70]. The DGLAP equation is solved directly in x-space 
using a method which is briefly illustrated below and which is accurate up to order a^. Our 
input distributions at a small scale will be specified below. We also analyse the corresponding 
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Figure 4.1: The leading order diagram for Higgs production by gluon fusion 




Figure 4.2: A typical NLO diagram for Higgs production by gluon fusion 

K-factors and the region of stability of the perturbative expansion by studying their variation 
under changes in all the relevant scales. It is shown that the NNLO corrections are sizeable while 
the region of reduced scale dependence is near the value fip = with around the same 
value but slightly higher. 



4.1 Higgs production at LHC 



The Higgs field, being responsible for the mechanism of mass generation, can be radiated off by 
any massive state and its coupling is proportional to the mass of the same state. At the LHC 
one of the golden plated modes to search for the Higgs is its production via the mechanism of 
gluon fusion. The leading order contribution is shown in Fig. (|4.1|) which shows that dependence 
of the amplitude is through the quark loop. Most of the contribution comes from the top quark, 
since this is the heaviest quark and has the largest coupling to the Higgs field. NLO and NNLO 
corrections have been computed in the last few years by various groups [71], [72]. A typical 
NLO correction is shown in Fig. (|4.2p . In the infinite mass limit of the quark mass in the loop 
(see ^73j for a review), an effective description of the process is obtained in leading order by the 



Lagrangian density 
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(4.1) 



being the contribution of Nh heavy fermion loops to the QCD beta function. This effective 
Lagrangian can be used to compute the radiative corrections in the gluon sector. A discussion of 
the NNLO approach to the computation of the gluon fusion contributions to Higgs production 
has been presented in [68], to which we refer for more details. We recall that in this chapter we 
presented a study for both scalar and pseudoscalar Higgs production, the pseudoscalar appearing 
in 2-Higgs doublets models. The diagonalization of the mass matrix for the Higgs at the minimum 
introduces scalar and pseudoscalar interactions between the various Higgs and the quarks, as 
shown from the structure of the operator O2 below in eq. (|4.3|) . In the large top-quark mass limit 
the Feynman rules for scalar Higgs production (H) can be derived from the effective Lagrangian 
density [75], [74], [76], 



(4.2) 



£H ^ = Gn ^^(x) 0{x) with 0{x) = C^^ix) C^^^^ix) 
whereas the production of a pseudo-scalar Higgs [T^, (A) is obtained from 
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(4.3) 



where $^(x) and ^^{x) represent the scalar and pseudo-scalar fields respectively and Uf denotes 
the number of light flavours. Ga'^ is the field strength of QCD and the quark fields are denoted 
by Qi. We refer the reader to [68] for further details. 

Using the effective Lagrangian one can calculate the total cross section of the reaction 



(4.4) 



where Hi and H2 denote the incoming hadrons and X represents an inclusive hadronic state 
and B denotes the scalar or the psudoscalar particle produced in the reaction. The total cross 



section is given by 

0"tot = o /jY2 ^ I) X] dxi dx2 fa{xi,fl^) fb{x2,H^) 

^ ' a,b=q,q,g '' ^ Jx/xi 

fx m? 

X B , 

\XiX2 /U^ 

2 

with x='^ , 5 = (Pi + P2)^ , (4.5) 

where the factor 1/{N'^ — 1) is due to the average over colour. The parton distributions /a(y, /i^) 
{a,b = q,q,g) depend on the mass factorization/renormahzation scale fi. Aq^ b denotes the 
partonic hard scattering coefficient computed with NNLO accuracy. 



4.2 The NNLO Evolution 

We summarize the main features of the NNLO DGLAP evolution. As usual we introduce singlet 
(+) and non-singlet (— ) parton distributions 

Uf 

^=m±c!^, «(^) = E^r^ (4-6) 

1=1 

whose evolution is determined by the corresponding equations 



for the singlet combination and a scalar one for the non-singlet case 
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dlogQ2^* 
The convolution product is defined by 



(4.7) 



I \x,Q^) = PNs{x,as{Q^))0qi{x,Q^). (4.8) 



[a b] {x) = t (^) b{y) = ^a{y)b ) . (4.9) 



We recall that the perturbative expansion, up to NNLO, of the kernels is 

P{x, a,) = a,P^^\x) + a2p(i)(x) + alP^^\x) + .... (4.10) 
where = Qs/(47r). In order to solve the evolution equations directly in 2:-space [6], [78], (see 



[3] for an NLO implementation of the method), we assume solutions of the form [79] 



^^ '^^ - n! a,(Q2)+«^(^)l. log 



+a?(Q^) log" ^ (4.11) 

for each parton distribution /, where Qq defines the initial evolution scale. The ansatz is in- 
troduced into the evolution equations and used to derive recurrence relations for its unknown 
coefficients A„,i?„,C„, involving polylogarithmic functions [HOlEl] which are then implemented 
numerically. 

This ansatz corresponds to a solution of the DGLAP equation accurate up to order (trun- 
cated solution). It can be shown that [79] this ansatz reproduces the solution of the DGLAP 
equation in (Mellin) moment space obtained with the same accuracy in o^. Modifications of 
this ansatz also allow to obtain the so-called "exact" solutions of the equations for the moments 
[2J. These second solutions include higher order terms in and can be identified only in the 
non-singlet case. Exact approaches also include an exact solution of the renormalization group 
equation for the /3-function, which embodies the effects of the coefficients (3o, (3i and f32 to higher 
order in o^. The term "exact" is, however, a misnomer since the accuracy of the solution is 
limited to the knowledge of the first three contributions to the expansion in the beta function 
and in the kernels. It can be shown both for exact and for the truncated solutions that solving 
the equations by an ansatz in x-space is completely equivalent to searching for the solution in 
moment space, since in moment space the recursion relations can be solved exactly [79] . 

A numerical comparison of our approach with that of [2] has been presented in where it is 
shown that at LO and NLO there is excellent agreement, while at NNLO there are discrepancies 
of a few percent (mainly in the Singlet case, and at very small and large x values). For a more 
detailed discussion we refer the reader to section 11 of [79]. 



4.3 Renormalization scale dependence 

For a better determination of the dependence of the perturbative cross section on the scales of a 
certain process it is important to keep these scales independent and study the behaviour of the 
corresponding hadronic cross section under their variation. In our case the two relevant scales 
are the factorization scale /ip and the renormalization scale fiR which can be both included in 
the evolution by a rearrangement of the evolution kernels up to NNLO. 

The study of the dependence of the solution upon the various scales is then performed in 
great generality and includes also the logarithmic contributions log{nF / fiR) coming from the 
hard scatterings given in [68] . where, however, only the specific point = I^r = rnn was 
considered. The separation of the scales should then appear not only in the hard scatterings but 
also in the evolution equations. This issue has been addressed in and can be reconsidered 



also from x-space [79] using the x-space logarithmic ansatz ()4.11|) . 

The scale dependence of the parton distribution functions is then expressed by a generalized 
DGLAP equation 

d 

ain/i^ = PiAx^hI'.^i'r) ® fj{x,njp,nji) , (4.12) 

where ;Ui? is now a generic factorization scale. 

Generally speaking, both the kernels and the PDF's have a dependence on the scales and 
liR, and formally, a comparison between these scales is always possible up to a fixed order by 
using the renormalization group equations for the running coupling a^- 

The renormalization scale dependence of the ansatz (14. lip that solves (|4.12|) is obtained quite 
straightforwardly by a Taylor expansion of the running coupling UsifJ.'j^) in terms of a<j(/i|j) [TS] 



(4.13) 



where the ^'jp dependence is included in the factor L = hi(/i|,//i|j), and the coefficients of the 
/3-function, (the Pi) are listed below [82], [EH], [H] 

/3o = yiVc - ^Tf, 

f3i = ^N^-^Ncnf-2CFnf, 

= ^iV3 + 2C^T, - '^CNcT, - Hi^iV^T, + fc.Tf + 'fNcTf. 

(4.14) 

As usual we have set 

= ^^ = ^iV^ = 3' ^/ = ^«"/=2"/' (4.15) 

where Nc is the number of colors and is the number of active flavors. This number is 
varied as we step into a region characterized by an evolution scale ^ larger than a specific quark 
mass (^ > nig). Also the NNLO matching conditions across flavor thresholds [HS], [M] are 
implemented. 

Since the perturbative expansion of eq. (|4.10|) contains powers of as(//|,) which can be related 
to the value of as(;U^) by (I4.13|) . from 

/^7^^«(x,;x^) = ^ (^)'^'/^(f)(x), (4.16) 

k=0 \ / 



substituting eq. (|4.13|] into (I4.16|) . we obtain the corresponding expression of the kernels orga- 



nized in powers of as(/i|j) up to NNLO, and it reads [2] 
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The implementation of the method in x-space is quite straightforward and allows us to 
perform a separate study of the predictions in terms oi fxp and ixr. 



4.4 Numerical Results 

The use of the NNLO evolution of the parton distributions together with the results of [68] 
allows us to provide accurate predictions for the total cross section for Higgs production. Here 
we summarize and discuss our numerical results. We use as inital conditions at low scales the 
sets of distributions given by MRST [88] and Alekhin [89]. Our final plots refer to center-of-mass 
energies which are reachable at the LHC, with 14 TeV being the largest one achievable in a not 
so distant future, and at the Tevatron, where we have selected the corresponding value as 2 TeV. 
We have also taken the Higgs mass ruH as a parameter in the prediction, with an interval of 
variability which goes from a light to a heavy Higgs (100 GeV to 300 GeV). Therefore ^p-, fJ'R 
and rriH are studied chosing various combinations of their possible values in the determination of 
total cross sections at leading {(Tlo), next-to-leading {(Jnlo), and next-to-next-to-leading order 
{ctnnlo)- We present both standard two-dimensional plots and also some three dimensional 
plots in order to characterize in detail the structure of the region of stability of the perturbative 
expansion. We have also evaluated the K-factors for the total Higgs cross section at NLO, defined 
by 

and at NNLO 

j^NNLO ^ <^NNLO _ ^^_^g^ 
(yNLO 

The study of the K-factors has been performed first by keeping the three scales equal {fJ-F = 
Ai_R = it^h) and then letting them vary around the typical value mn- A second set of studies 
has been performed by taking typical values of mn and varying the value of the renormalization 
scale. 



^ Based on the article [87l 



4.4.1 The errors on the cross sections 

We present in Figs. (|4.3|1 and (|4.4|1 the LO, NLO and NNLO results for the total Higgs cross 
sections at the LHC {VS = 14 TeV) and at the Tevatron {^/S = 2 TeV), with the corresponding 
errors, by setting the condition = fip. We have chosen to compute these only for two figures, 
as an illustration of the size of the errors on the parton distribution functions compared to the 
best fits, since these are smaller than the variation induced when moving from one perturbative 
order to the next. The numerical determination of the errors is computationally very intensive 
and has been performed on a cluster. 

In the figures on the left the cross sections obtained using MRST input are represented by 
a solid line, and the ones obtained using Alekhin's input, by a dashed line. In the figures on 
the right we present a plot of the difference between the values of the MRST cross sections and 
Alekhin's cross sections for each perturbative order, with the respective errors. The calculation of 
the error bands has been done following the usual theory of the linear propagation of the errors. 
Starting from the errors on the PDFs known in the literature (see [89], [90]), we have generated 
different sets of cross sections. Then, the error on the cross section has been calculated using 
the formula 



where is the k-th cross section belonging to a certain set, and is the number of free 
parameters, which is 15 for MRST and 17 for Alekhin. 

The PDFs with the related error analysis are available at all orders for the Alekhin input but 
only at NLO for the MRST's input (Figs. g31)-c and (lOl-c). 

When in Figs. (I4.3|l . (|4.4p . and (|4.5|) we plot more than one line for a single set the lower line 
is the minimal value (best fit minus error) and the upper line is the maximal value (best fit plus 
error) . 

The LO cross sections increase by a factor of approximately 100 as we change the energy 
from 2 TeV Figs. (14.41) . to 14 TeV Figs. (|4.3p . and sharply decrease as we raise the mass of the 
Higgs boson. At 14 TeV the range of variation of aLO is between 30 and 5 pb, with the highest 
value reached for mn = 100 GeV. 

In the same figures we compare LO, NLO and NNLO cross sections at these two typical 
energies. It is quite evident that the role of the NLO corrections is to increase by a factor of 
approximately 2 the LO cross section bringing the interval of variation of (TjqLo between 60 and 
10 pb, for an increasing value of niH- NNLO corrections at 14 TeV increase these values by an 
additional 10 per cent compared to the NLO prediction, with a growth which is more pronounced 
for the set proposed by Alekhin. 

Comparing the results computed using Alekhin and MRST's inputs, for the LHC case {y/S = 
14 TeV) we observe that at LO (see Fig. (|4.3|) -b) the two sets give results which are compatible 




(4.20) 



within the error bands for rriH < 150 GeV, while for larger values of mn we observe only small 
differences between the two. At NLO, Fig. (|4.3|) -d. where the error analysis is available for both 
sets, the results are compatible within the error bands for mn < 190 GeV, and we have small 
differences for larger values of the Higgs mass. For the NNLO case. Fig. (I4.3ll -f. we notice that 
there are sensible differences among the two sets. 

By a similar inspection of Figs. (|4.4l) -b. (I4.4p -d and (I4.4l) -f. we notice that at the Tevatron 
energy of = 2 TeV, the two PDFs sets give quite different predictions at all the three orders. 

The numerical values of the total cross sections and the K-factors as a function of the center 
of mass energy with the respective errors have also been reported in Table 14.11 in the case of 
Alekhin's inputs. 



4.4.2 K-factors 

A precise indication on the impact of the NLO /NNLO corrections and the stability of the pertur- 
bative expansion comes from a study of the K-factors K^^'-' and K^^^^ , defined above. From 
the plots in Figure US] the different behaviour of the predictions derived from the two models for 
the parton distributions is quite evident. At 14 TeV the NLO K-factors from both models are 
large, as expected, since the LO prediction are strongly scale dependent. The increase of (Jnlo 
compared to the LO predictions is between 65 and 90 per cent. 

In Fig. (|4.5p -a one can observe that the impact of the NLO corrections to the LO result 
predicted by both sets increases for an increasing tuh, with the corrections predicted by Alekhin 
being the largest ones. The trend of the MRST model in the NNLO vs NLO case Fig. ()4.5l) -b is 
similar, for an increasing Higgs mass, ranging from 1.16 to 1.22, while Alekhin X^^^'-^-factor is 
approximately constant around a value of 1.21. 

The evaluation of the overall impact of this growth on the size of these corrections should, 
however, also keep into consideration the fact that these corrections are enhanced in a region 
where the cross section is sharply decreasing (Figs. (|4.3p and (|4.4I) ). 

Moving to Tevatron energy we notice that all the K-factors are larger than in the LHC case. 
For a MRST input we continue to observe a growth of K'^^^ and K^^^O increasing mu, 

while for the Alekhin case this trend is slower for K^^^ , and it is even reversed for x^^^^ . 
Unlike the LHC situation, the Alekhin's K-factors are smaller than MRST ones at Tevatron 
energies. 



4.4.3 Renormalization/factorization scale dependence 

Now we turn to an analysis of the dependence of our results on /xj? and fiR. In Figs. 14.61 we 
perform this study by computing o" as a function of the Higgs mass for an incoming energy of 2 
and 14 TeV and choose 

lAi = 2l^F I^F = 2mH- (4.21) 



We have seen that for a typical Higgs mass around 100 GeV and 14 TeV of energy (for niH = 
/^F = (J'r) the cross section doubles when we move from LO to NNLO, and a similar trend 
is also apparent if we fix the relations among the scales as in (|4.2ip . In this case, however, 
the impact of the NLO and NNLO corrections is smaller, a trend which is apparently uniform 
over the whole range of the Higgs mass explored. For niH = 100 GeV the scalar cross section 
o'NNLO is around 58 pb for coincident scales, while a different choice, such as (|4.21|) lowers 
it to approximately 45 pb. At Tevatron energies the variations of the cross section with the 
changes of the various scales are also sizeable. In this case for mn = 100 GeV the LO, NLO 
and NNLO predictions (0.6,1.2 and 1.6 pb respectively) change approximately by 10-20 per 
cent if we include variations of the other scales as well. A parallel view of this trend comes 
from the study of the dependence of the K-factors. This study is presented in Figs. 14.71 The 
interval of variation of the K-factors is substantially the same as for coincident scales, though 
the trends of the two models [88] and [89] is structurally quite different at NLO and at NNLO, 
with several cross-overs among the corresponding curves taking place for mn around 200 GeV. 
Another important point is that the values of K^^^'-^ are, of course smaller than K^^^ over 
all the regions explored, signaling an overall stability of the perturbative expansion. We show in 
Figs. 11-14 3-dimensional plots of the cross section and the corresponding K-factors as functions 
of the factorization and renormalization scales. Notice that as we move from LO to higher 
orders the curvature of the corresponding surfaces for the cross sections change from negative to 
positive, showing the presence of a plateau when the scales are approximately equal. 



4.4.4 Stability and the Choice of the Scales 

The issue of determining the best of values of mn, fJ'F and fiR in the prediction of the total 
cross section is a rather important one for Higgs searches at LHC. We have therefore detailed in 
Figs. 14.81 and 14.91 a study of the behaviour of our results varying the renormalization scale fiR at 
a fixed value of the ratio between and mn. In these figures we have chosen two values for the 
ratio between these two scales. Apart from the LO behaviour of the scalar cross section, which 
is clearly strongly dependent on the variation of both scales (see Figs. 14.81 (a) and (d)), and does 
not show any sign of stability since the cross section can be drastically lowered by a different 
choice of hf, both the NLO and the NNLO predictions show instead a clear region of local 
stability for fiji > fip but not too far away from the "coincidence region" fip = t^R = fnu- This 
can be illustrated more simply using Fig. 14.81 (b) as an example, where we have set the incoming 
energy of the p-p collision at 14 TeV. In this case, for instance, we have chosen mu = fJ-F = 100 
GeV (C = 1), and it is clear from the plots that a plateau is present in the region of ~ 130 
GeV. Similar trends are also clearly visible at NNLO, though the region of the plateau for the 
scalar cross section is slightly wider. Also in this case it is found that the condition /i^ > fJ-F 
generates a reduced scale dependence. Away from this region the predictions show a systematic 
scale dependence, as shown also for the choice of C = 1/2 in the remaining figures. In Figs. 14. 91 we 
repeat the same study for the K-factors, relaxing the condition on the coincidence of all the scales 
and plotting the variations of K^^^ and K^^^^ in terms of jiR. In the case mn = ^J^f = 100 



GeV the plateau is reached for ^ir ~ 150 GeV for K'^^^ and ^ir ~ 200 GeV for k'^nlu_ ^j^g 
first case the NLO corrections amount to an increase by 100 per cent compared to the LO result, 
while the NNLO corrections modify the NLO estimates by about 20 per cent (MRST). Similar 
results are obtained also iox hf = 50 GeV. In this case, at the plateau, the NLO corrections are 
still approximately 100 per cent compared to the LO result and the NNLO corrections increase 
this value by around 15 per cent (MRST). 



4.4.5 Energy Dependence 

The energy dependence of the NNLO predictions for the total cross sections and the corre- 
sponding K-factors at the LHC are shown in Figs. I4.15ll4.17l where we have varied the ratio 
C = fip/^H and k = fij^/ fip in order to illustrate the variation of the results. The cross sec- 
tions increase sharply with energy and the impact of the NNLO corrections is significant. The 
K-factors, in most of the configurations chosen, vary between 1 and 2.2. We have chosen the 
MRST input. The behaviour of the K-factors is influenced significantly by the choice of the ratio 
(fc) between ^r and ij,f- In particular, in Figs. 14.161 the NNLO K-factors increase with ^/S for 
k = 2 , the center of mass energy, which is not found for other choices of scales. The case k = 1/2 
is close in behaviour to the coincident case fij^ = Hp. The overall stability of the K-factors is 
clearly obtained with the choice k = 1. We have finally included in Tables 2-10 our numerical 
predictions in order to make them available to the experimental collaborations. 



4.5 Conclusions 



A study of the NNLO corrections to the cross section for Higgs production has been presented. 
We have implemented the exact three-loop splitting functions in our own parton evolution code. 
We used as initial conditions (at small scales) the boundary values of Martin, Roberts, Thorne 
and Stirling and of Alekhin. This study shows that the impact of these corrections are important 
for the discovery of the Higgs and for a reconstruction of its mass. The condition of stability of 
the perturbative expansion is also quite evident from these studies and suggests that the optimal 
choice to fix the arbitrary scales of the theory are near the coincidence point, with fiR in the 
region of a plateau. The determination of the plateau has been performed by introducing in 
the perturbative expansion and in the evolution a new independent scale (/^j?), whose variation 
allows to accurately characterize the properties of the expansion in a direct way. 

While this thesis was being completed several authors have presented studies of the total 
Higgs cross section based on threshold resummation of soft or soft-plus- virtual logarithms, see 
[91],[92],[93], [94]- Our work is based on exact NNLO partonic cross sections. Another relevant 
paper which recently appeared is [95] . 





Figure 4.3: Cross sections for the scalar Higgs production at the LHC with fi^ = jx^p = mjj. 
When available (Alekhin at all orders and MRST at NLO) the error bands are shown. See sec. 
5 for comments. 
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Figure 4.4: Like in Fig. (14. 3 1) but for the Tevatron. In Figs, (a), (c) and (e) the error bands are 
so tiny that the Hnes look superimposed. 
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Figure 4.5: K-factors for the scalar Higgs at NNLO/NLO and NLO/LO with jLt| = 
When available (Alekhin at all orders and MRST at NLO) the error bands are shovi^n. 
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Figure 4.6: Cross sections for the scalar Higgs production at the LHC and Tevatron with = 
(l/2)/x|n and fip = 2mif 
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Figure 4.7: K-factors for the scalar Higgs at NNLO/NLO and NLO/LO with /xfj = {l/2)fij, and 




Figure 4.8: Cross sections for the scalar Higgs production at the LHC as a function of /ir, with 
IJ.F = CruH and ruH = 100 GeV 
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Figure 4.9: K-factors for the scalar Higgs production at the LHC, NNLO/NLO and NLO/LO as 
a function of fiR, with jj^f = Ctuh and niH = 100 GeV 
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Figure 4.10: Three-dimensional graphs for the LO cross sections for the scalar Higgs production 
at the LHC, as a function of jiR and with with a fixed value of mn 
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Figure 4.11: Three-dimensional graphs for the NLO cross sections for the scalar Higgs production 
at the LHC, as a function of jiR and with with a fixed value of mn 
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Figure 4.12: Three-dimensional graphs for the NNLO cross sections for the scalar Higgs produc- 
tion at the LHC, as a function of /ir and with /if with a fixed value of mn 




Figure 4.13: Three-dimensional graphs for the K-factor (Ti\jlo/^lo for the scalar Higgs produc- 
tion at the LHC, as a function of fiR and fip and at a fixed value of mn- 
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Figure 4.14: Three-dimensional graphs for the i^-factor (Tnnlo/'^nlo for the scalar Higgs pro- 
duction at the LHC, as a function of iir and fip and at a fixed value of mn- 
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Figure 4.15: Cross sections and if-factors for the scalar Higgs production at the LHC as a 
function of y/S with fip = Cmn, with fij^ = kfi^ and niH = 114 GeV. MRST inputs have been 
used. 
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Figure 4.16: Cross sections and if-factors for the scalar Higgs production at the LHC as a 
function of y/S with fip = Cmn, with fij^ = kfi^ and niH = 114 GeV. MRST inputs have been 
used. 
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Figure 4.17: Cross sections and if- factors for the scalar Higgs production at the LHC as a 
function of y/S with fip = Cmn, with fij^ = kfi^ and niH = 114 GeV. MRST inputs have been 
used. 
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8.5 


9.75 ± 0.19 


18.04 ±0.18 


22.1 ± 0.3 


1.85 ± 0.04 


1.224 ± 0.020 


9.0 


10.8 ±0.2 


19.9 ±0.2 


24.4 ± 0.3 


1.84 ±0.04 


1.22 ± 0.02 


9.5 


12.0 ±0.3 


21.9 ±0.3 


26.7 ±0.4 


1.83 ±0.04 


1.22 ±0.02 


10.0 


13.1 ±0.3 


23.9 ±0.3 


29.1 ±0.5 


1.82 ±0.05 


1.22 ± 0.03 


10.5 


14.3 ± 0.4 


25.9 ±0.4 


31.5 ±0.6 


1.81 ± 0.05 


1.21 ±0.03 


11.0 


15.6 ±0.4 


28.1 ±0.4 


34.0 ± 0.7 


1.80 ±0.06 


1.21 ±0.03 


11.5 


16.8 ±0.5 


30.2 ±0.5 


36.6 ±0.8 


1.80 ±0.06 


1.21 ±0.03 


12.0 


18.1 ±0.5 


32.4 ±0.6 


39.2 ±0.9 


1.79 ± 0.06 


1.21 ±0.04 


12.5 


19.4 ±0.6 


34.6 ± 0.7 


41.8 ± 1.1 


1.78 ± 0.07 


1.21 ± 0.04 


13.0 


20.8 ±0.7 


36.9 ±0.8 


44.5 ± 1.2 


1.77 ±0.07 


1.21 ±0.04 


13.5 


22.2 ±0.8 


39.2 ±0.9 


47.2 ± 1.4 


1.77 ±0.07 


1.20 ±0.04 


14.0 


23.6 ±0.9 


41.6 ± 1.0 


50.0 ± 1.5 


1.76 ± 0.08 


1.20 ±0.05 



Table 4.1: Values of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of VS with //j? = mn, with fj,^ = //^ and uih = 114 GeV for Alekhin, with errors. 



fa 




<^NLO 


C^NNLO 


j^N LO 


j^NN LO 


2.0 


0.4155 


0.8670 


1.242 


2.087 


1.433 


2.5 


0.7410 


1.521 


2.084 


2.053 


1.370 


3.0 


1.153 


2.335 


3.093 


2.025 


1.325 


3.5 


1.645 


3.292 


4.248 


2.001 


1.290 


4.0 


2.212 


4.380 


5.529 


1.980 


1.262 


4.5 


2.847 


5.587 


6.924 


1.962 


1.239 


5.0 


3.547 


6.903 


8.419 


1.946 


1.220 


5.5 


4.308 


8.318 


10.01 


1.931 


1.203 


6.0 


5.125 


9.826 


11.68 


1.917 


1.189 


6.5 


5.995 


11.42 


13.42 


1.905 


1.175 


7.0 


6.916 


13.09 


15.23 


1.893 


1.163 


7.5 


7.885 


14.84 


17.11 


1.882 


1.153 


8.0 


8.899 


16.66 


19.05 


1.872 


1.143 


8.5 


9.956 


18.55 


21.04 


1.863 


1.134 


y.u 


1 1 C\K 


on /lo 


no 


1 St^/l 
i.o04 


1 1 07 


9.5 


12.19 


22.50 


25.18 


1.846 


1.119 


10.0 


13.37 


24.56 


27.31 


1.837 


1.112 


10.5 


14.58 


26.68 


29.49 


1.830 


1.105 


11.0 


15.83 


28.85 


31.71 


1.822 


1.099 


11.5 


17.11 


31.06 


33.97 


1.815 


1.094 


12.0 


18.42 


33.32 


36.26 


1.809 


1.088 


12.5 


19.76 


35.62 


38.59 


1.803 


1.083 


13.0 


21.13 


37.97 


40.95 


1.797 


1.078 


13.5 


22.53 


40.36 


43.33 


1.791 


1.074 


14.0 


23.96 


42.78 


45.75 


1.785 


1.069 



Table 4.2: \ alues of the cross sections and i^-factors for the scalar Higgs production at the LHC 
as a function of ^/S with hf = niH, with jj,^ = jij^ and ruH = 114 GeV. MRST inputs have been 
used. 



fa 




ctnlo 


O'NNLO 






2.0 


0.3029 


0.6871 


1.081 


2.268 


1.573 


2.5 


0.5508 


1.221 


1.830 


2.217 


1.499 


3.0 


0.8700 


1.893 


2.733 


2.176 


i AAA 

1.444 


3.5 


1.256 


2.690 


3.771 


2.142 


1.402 


4.0 


1.706 


3.602 


4.928 


2.111 


1.368 


4.5 


2.216 


4.619 


6.191 


2.084 


1.340 


5.0 


2.781 


5.733 


7.549 


2.061 


1.317 


5.5 


3.400 


6.937 


8.993 


2.040 


1.296 


6.0 


4.069 


8.225 


10.52 


2.021 


1.279 


6.5 


4.786 


9.590 


12.11 


2.004 


1.263 


7.0 


5.548 


11.03 


13.77 


1.988 


1.248 


7.5 


6.354 


12.53 


15.49 


1.972 


1.236 


8.0 


7.201 


14.11 


17.27 


1.959 


1.224 


8.5 


8.088 


15.74 


19.10 


1.946 


1.213 
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1 7 /IQ 
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1 OC\A 


9.5 


9.974 


19.17 


22.90 


1.922 


1.195 


10.0 


10.97 


20.97 


24.87 


1.912 


1.186 


10.5 


12.00 


22.82 


26.88 


1.902 


1.178 


11.0 


13.06 


24.71 


28.93 


1.892 


1.171 


11.5 


14.16 


26.65 


31.01 


1.882 


1.164 


12.0 


15.28 


28.63 


33.13 


1.874 


1.157 


12.5 


16.44 


30.65 


35.28 


1.864 


1.151 


13.0 


17.62 


32.71 


37.47 


1.856 


1.146 


13.5 


18.83 


34.81 


39.68 


1.849 


1.140 


14.0 


20.07 


36.95 


41.92 


1.841 


1.135 



Table 4.3: \ alues of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of with jip = "^ruH, with jip = /j,^ and niH = 114 GeV. MRST inputs have 
been used. 
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VO 






O'NNLO 




j^NNLO 


2.0 


0.5817 


1.098 


1.396 
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888 


1.271 
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1.014 


1.904 


2.328 


1 


878 


1.223 


3.0 


1.551 


2.896 


O A Ar\ 

3.440 


1 


867 


1.188 


3.5 


2.182 


4.054 


4.707 


1 


858 


1.161 


4.0 


2.899 


5.362 


6.111 


1 


850 


1.140 


4.5 


3.695 


6.804 


7.635 


1 


841 


1.122 


5.0 


4.563 


8.369 


9.266 


1 


834 


1.107 


5.5 


5.498 


10.05 


10.99 


1 


828 


1.094 


6.0 


6.496 


11.83 


12.81 


1 


821 


1.083 


6.5 


7.552 


13.70 


14.71 


1 


814 


1.074 


7.0 


8.662 


15.67 


16.67 


1 


809 


1.064 


7.5 


9.824 


17.71 


18.71 


1 


803 


1.056 


8.0 


11.03 


19.84 


20.81 


1 


799 


1.049 


8.5 


12.29 


22.03 


22.97 


-1 

1 


793 


1.043 


y.u 


io.oy 




oc; 1 fi 
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1 


788 


1 nQfi 

l.UoD 


9.5 


14.93 


26.63 


27.44 


1 


784 


1.030 


10.0 


16.31 


29.02 


29.75 


1 


779 


1.025 


10.5 


17.72 


31.46 


32.10 


1 


775 


1.020 


11.0 


19.18 


33.96 


34.50 


1 


771 


1.016 


11.5 


20.66 


36.51 
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Table 4.4: \ alues of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of y/S with //j? = (l/2)mff , with ij!p = ij?^ and mn = 114 GeV. MRST inputs have 
been used. 
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Table 4.5: \ alues of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of with jip = "niH, with = 2//|j and uih = 114 GeV. MRST inputs have 
been used. 
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Table 4.6: \ alues of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of \fS with fip = "^"mH, with /j^p = 2iJ?^ and uih = 114 GeV. MRST inputs have 
been used. 
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Table 4.7: \ alues of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of y/S with jip = (l/2)mH, with n'p = 2//^ and tuh = 114 GeV. MRST inputs 
have been used. 



fa 
V o 




<^NLO 


(^NNLO 






2.0 


0.4899 


1.276 


1.635 


2.605 


1.281 


2.5 


0.8641 


2.083 


2.567 


2.411 


1.232 


3.0 


1.333 


3.022 


3.622 


2.267 


1.199 


3.5 


1.890 


4.069 


4.775 


2.153 


1.174 


4.0 


2.526 


5.206 


6.009 


2.061 


1.154 


4.5 


3.237 


6.419 


7.312 


1.983 


1.139 


5.0 


4.016 


7.698 


8.673 


1.917 


1.127 


5.5 


4.858 


9.034 


10.08 


1.860 


1.116 


6.0 


5.761 


10.42 


11.54 


1.809 


1.107 


6.5 


6.719 


11.85 


13.03 


1.764 


1.100 


7.0 


7.730 


13.32 


14.55 


1.723 


1.092 


7.5 


8.790 


14.82 


16.11 


1.686 


1.087 


8.0 


9.898 


16.35 


17.69 


1.652 


1.082 


8.5 


11.05 


17.91 


19.30 


1.621 


1.078 


y.u 


10 0/1 


1 Q t^n 
iy.ou 


on 


1 KQ'i 

i.oyo 


1 n7Q 

l.U / o 


9.5 


13.48 


21.11 


22.57 


1.566 


1.069 


10.0 


14.75 


22.74 


24.24 


1.542 


1.066 


10.5 


16.06 


24.39 


25.92 


1.519 


1.063 


11.0 


17.41 


26.06 


27.61 


1.497 


1.059 


11.5 


18.79 


27.74 


29.32 


1.476 


1.057 


12.0 


20.20 


29.44 


31.05 


1.457 


1.055 


12.5 


21.64 


31.15 


32.78 


1.439 


1.052 


13.0 


23.11 


32.87 


34.53 


1.422 


1.051 


13.5 


24.62 


34.60 


36.28 


1.405 


1.049 


14.0 


26.15 


36.35 


38.04 


1.390 


1.046 



Table 4.8: \ alues of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of y/S with = mn, with = (1/2)//^ and uih = 114 GeV. MRST inputs have 
been used. 
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Table 4.9: \ alues of the cross sections and K-factors for the scalar Higgs production at the LHC 
as a function of y/S with jip = 2mff, with n'p = (l/2)//|j and tuh = 114 GeV. MRST inputs 
have been used. 
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Table 4.10: Values of the cross sections and ii'-factors for the scalar Higgs production at the 
LHC as a function of with hf = {l/2)mH, with = (l/2)//| and mn = 114 GeV. MRST 
inputs have been used. 



Chapter 5 



Deeply Virtual Neutrino Scattering 
(DVNS) 

5.1 Introduction 

In this second part of the thesis, after a brief summary of the generalities of the Deeply Virtual 
Compton Scattering (DVCS) process, we discuss its generalization to the case of neutral currents 
and the phenomenological implications concerning the Physics of Neutrino-Nucleon interactions 
in the few GeV's range of energies. 

Exclusive processes mediated by the weak force are an area of investigation which may gather 
a wide interest in the forthcoming years due to the various experimental proposals to detect neu- 
trino oscillations at intermediate energy using neutrino factories and superbeams [10]. These 
proposals require a study of the neutrino-nucleon interaction over a wide range of energy start- 
ing from the elastic/quasi-elastic domain up to the deep inelastic scattering (DIS) region (see 
[M],[II],|nZ],|nE],|M] for an overall overview). However, the discussion of the neutrino nucleon 
interaction has, so far, been confined either to the DIS region or to the form factor/nucleon 
resonance region, while the intermediate energy region, at this time, remains unexplored also 
theoretically. Clearly, to achieve a "continuos" description of the underlying strong interaction 
dynamics, from the resonant to the perturbative regime, will require considerable effort, since it 
is experimentally and theoretically difficult to disentangle a perturbative from a non-perturbative 
dynamics at intermediate energy, which appear to be superimposed. This is best exemplified - at 
least in the case of electromagnetic processes, such as Compton scattering - in the dependence of 
the intermediate energy description on the momentum transfer [36j . In this respect, the interac- 
tion of neutrinos with the constituents of the nucleon is no different, once the partonic structure 
of the target is resolved. From our viewpoint, the presence of such a gap in our knowledge well 
justifies any attempt to improve the current situation. 

We recall that DVCS has been extensively studied in the last few years for electromagnetic 
interactions. The extension of DVCS to the case of neutral currents is presented here, while the 
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Figure 5.1: Leading hand-bag diagrams for the process 

charged current version will be presented in the next chapter. Also, in this chapter we will just 
focus on the lowest order contributions to the process, named by us Deeply Virtual Neutrino 
Scattering (DVNS) in order to distinguish it from standard DVCS, while a renormalization group 
analysis of the factorized amplitude, which requires an inclusion of the modifications induced by 
the evolution will also be presented elsewhere. The application of the formalism that we develop 
here also needs a separate study of the isoscalar cross sections together with a detailed analysis of 
the various experimental constraints at neutrino factories in order to be applicable at forthcoming 
experiments. 



5.2 The Generalized Bjorken Region and DVCS 

Compton scattering has been investigated in the near past by several groups, since the original 
work by Ji and Radyushkin [1001 llOll I102| . Previous work on the generalized Bjorken region, 
which includes DVCS and predates the "DVCS period" can be found in |103| . 

A pictorial description of the process we are going to illustrate is given in Fig 15.11 where a 
neutrino of momentum I scatters off a nucleon of momentum Pi by an interaction with a neutral 
current; from the final state a photon and a nucleon emerge, of momenta q2 and P2 respectively, 
while the momenta of the final lepton is The process is described in terms of new constructs 
of the parton model termed generalized parton distributions (GPD) or also non-forward (off- 
forward) parton distributions. We recall that the regime for the study of GPD's is characterized 
by a deep virtuality of the exchanged photon in the initial interaction {e+p ^ e+p + j) ( 
2 GeV^), with the final state photon kept on-shell; large energy of the hadronic system {W'^ > 6 
GeV^) above the resonance domain and small momentum transfers \t\ < 1 GeV^. The process 
suffers of a severe Bethe-Heitler (BH) background, with photon emission taking place from the 
lepton. Therefore, in the relevant region, characterized by large and small t, the dominant 
Bethe-Heitler background (~ 1/t) and the 1/Q behaviour of the DVCS scattering amplitude 
render the analysis quite complex. From the experimental viewpoint a dedicated study of the 
interference BH-VCS is required in order to explore the generalized Bjorken region, and this is 
done by measuring asymmetries. Opting for a symmetric choice for the defining momenta, we 
use as independent variables the average of the hadron and gauge bosons momenta 



with —A = P2 — Pi being the momentum transfer. Clearly 

p . A = 0, t = A^ p2 ^ - ^ (5.2) 

and M is the nucleon mass. There are two scaling variables which are identified in the process, 
since 3 scalar products can grow large in the generalized Bjorken limit: q'^ , A ■ q, P ■ q. 

The momentum transfer t = A^ is a small parameter in the process. Momentum asymmetries 
between the initial and the final state nucleon are measured by two scaling parameters, ^ and r/, 
related to ratios of the former invariants 



2P-q ' 2P-q 

where ^ is a variable of Bjorken type, expressed in terms of average momenta rather than nucleon 
and Z-boson momenta. The standard Bjorken variable x = —q\/{2Pi ■ qi) is trivially related to 
^ in the t = hmit. In the DIS hmit (Pi = P2) = and x = ^, while in the DVCS hmit r] = ^ 
and x = 2^/(1 + ^), as one can easily deduce from the relations 



.f=^l + ^J.- + 4, qi=[l-'jjt + ^^. (5.4) 

We introduce also the inelasticity parameter y = Pi ■ I /{Pi ■ qi) which measures the fraction of 
the total energy that is transferred to the final state photon. Notice also that ^ = measures 
the ratio between the plus component of the momentum transfer and the average momentum. 
A second scaling variable, related to ^ is C = ^^/-Pi^, which coincides with Bjorken x {x = C) 
when t = 0. ^, therefore, parametrizes the large component of the momentum transfer A, which 
can be generically described as 

A = 2^P + A (5.5) 
where all the components of A are 0{\/t) |1Q4| . 



5.3 DIS versus DVNS 



In the study of ordinary DIS scattering of neutrinos on nucleons (see Fig. 15. 2|) . the relevant 
current correlator is obtained from the T-product of two neutral currents acting on a forward 
nucleon state of momentum Pi 



= u{l')r{-l + As\n^ew+l5)u{l) 



(5.6) 



where 9w is the Weinberg angle, I and I' are the initial and final-state lepton. The relevant 
correlator is given by 

T^AqIi^) = I j d^ze^'^-^(Pi|r(4(OJ^(0)|Pi) (5.7) 

with V = E — E' being the energy transfered to the nucleon, Pi is the initial-state nucleon 
4-momenta and gi = / — Z' is the momentum transfered. The hadronic tensor W^j^u is related to 
the imaginary part of this correlator by the optical theorem. We recall that for an inclusive elec- 
troweak process mediated by neutral currents the hadronic tensor (for unpolarized scatterings) 
is identified in terms of 3 independent structure functions at leading twist 

Wf,^, = l^-g^^ + jWi{v,Q) + le^^^XaqiPi (5.8) 

where transversality of the current is obvious since Pi^ = — g^Pi ■ qi/qf- 

The analysis at higher twists is far more involved and the total number of structure func- 
tions appearing is 14 if we include polarization effects. These are fixed by the requirements of 
Lorenz covariance and time reversal invariance, neglecting small CP-violating effects from the 
CKM matrix. Their number can be reduced to 8 if current conservation is imposed, which is 
equivalent to requiring that contributions proportional to non-vanishing current quark masses 
can be dropped (see also |128j ). We recall that the DIS limit is performed by the identifications 

MWiiQ^u) = Fi{x,Q^) 
vW2{Q\u) = F2{x,Q^) 

uWsiQ^u) = F3{x,Q^), (5.9) 

in terms of the standard structure functions Fi, F2 and F3. There are various ways to express 
the neutrino-nucleon DIS cross section, either in terms of and the energy transfer, in which 
the scattering angle 9 is integrated over, or as a triple cross section in {Q'^,1',6), or yet in terms 
of the Bjorken variable x, inelasticity y and the scattering angle {x,y,6). This last case is close 
to the kinematical setup of our study. In this case the differential Born cross section in DIS is 
given by 

i^-'4^"'iQ"^'-rwr. (5.10) 

i 

the index i denotes the different current contributions, ( i = \ ,\jZ\,\Z\'^ for the neutral 
current) and a denotes the fine structure constant. By 9 we indicate the azimuthal angle of the 
final-state lepton, while y = {Pi ■ qi)/{l • Pi) is the inelasticity parameter, and = —qf- The 
factors r]i(Q'^) denote the ratios of the corresponding propagator terms to the photon propagator 




Figure 5.2: Leading diagram for a generic DIS process 



squared, 



2\/27ra Q^ + Mf 

(5.11) 

where Gp is the Fermi constant and Mz is the mass of the Z boson while the leptonic tensor 
has the form 

^l^u = E [^(^'' ^'hM + all^Hk, A)] * uik', X'h,{g'^ + 5if75)n(A;, A). (5.12) 

A' 

In the expression above A and A', denote the initial and final-state helicity of the leptons. The 
indices ii and 12 refer to the sum appearing in eq. (|5.10h 



9v = 1' 9l = 0, 

= -^ + 2sm^ 9w, 9a = h 



In the case of neutrino/nucleon interaction mediated by the neutral current the dominant diagram 
for this process appears in Fig. 15.21 A similar diagram, with the obvious modifications, describes 
also charged current exchanges. We recall that the unpolarized cross section is expressed in terms 
of Fi and F2, since F^ disappears in this special case, and in particular, after an integration over 
the scattering angle of the final state neutrino one obtains 



2-5^(<7^)%I^^I(Q2) (2X/F1 + 2(1 - . - ^)F2^ . (5.14) 



dxdy 

where S = 2MEi, is the nucleon-neutrino center of mass energy. 
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Figure 5.3: The cross section of a neutrino process of DIS-type at x ^ 0.1 with neutral current 
at ultrahigh energy 



We also recall that in this case the cross section in the parton model is given by 
(fa G%ME^ ( Ml 



'^{q^TMi) [^Ax^Q') + ^fM')i^-yf] (5-15) 



where q^{x,Q'^) and (f{x,Q^) are linear combinations of parton distributions 



with 



+ 



+ 



Uy{x,Q'^) + dy{x,Q'^) u{x,Q'^) + d{x,Q 



+ 



u{x,Q^)+d{x,Q'^ 



Uy{x,Q'^) + dy{x,Q'^) u{x,Q'^) + d{x,Q 



+ 



u{x,Q'^) + d{x,Q 



2\1 



(5.16) 



Lu = 1 — - sin 9w ■ 



Ru = —- sin 9w , 



Ld = -1 + 3 sin Ow 

2 2 
-Rd = 2 sm 9w 



(5.17) 



and we have identified the sea contributions Us and ds with u and d rispectively. 

Let's now move to the nonforward case. Here, when a real photon is present in the final state, 



the relevant correlator is given by 

T,,{qlv)=^ I d'ze'^-^P-^\T{4{-z/2)J^{z/2)\P + ^). (5.18) 

The dominant diagrams for this process appears in Fig l5.ll We impose Ward identities on both 
indices, which is equivalent to requiring that terms proportional to the quark masses in d ■ Jz 
are neglected. This approximation is analogous to the one performed in the forward case in 
order to reduce the structure functions from 8 to 3 (in the absence of any polarization), imposing 
symmetric trasversality conditions on the weak currents 

[P''-^^T^. = (p/^ + ^) T^, = 0. (5.19) 

The leading twist contribution to DVNS is obtained by performing a collinear expansion of the 
loop momentum of the hand-bag diagram and neglecting terms of order 0(A^/Q^) and M^/Q^. 
Transversality is satisfied at this order. Violation of transversality condition in the hand-bag 
approximation is analogous to the DVCS case, where it has been pointed out that one has 
to include systematically "kinematical" twist-3 operators, which appear as total derivatives of 
twist-2 operators |104| |105| in order to restore it. 

For the parametrization of the hand-bag diagram (Fig. 15. ip we use the light-cone decompo- 
sition in terms of 2 four- vectors (n, n), where 



= A(l, 0,0,1) 
= ^(1,0,0,-1) 
n = n = , h ■ n = 1. 



At the same time we set 



M 

= (l+^)n^ + (l-e)-n'^-^ 

"/Tf 

P,^ = {l+^)n^ + {l+^)-n'^ + ^ 

p = {k ■ n) h^' + {k ■ h) n^" + k^^ 
Is? 

jf = -— (5.20) 



with P"^ = M . We will also use the notation —q\ = for the invariant mass of the virtual 
Z boson and we will denote by q the average gauge bosons momenta respectively. 



After a collinear expansion of the loop momentum we obtain 



57^ [D, - 7^) rEl^'gdM^k)] (5.21) 
where we have used the following notations 



o^, gd = o^' 9 = 7— — 
o o 4 cos I 



8 4 
C4 = 1 - - sin^ Ow, ^ ~ 3 



(5.22) 

^ -a4 +qi 



{k — aA + qif' + ie 



= ^ +^('-°) ,5.23) 

where the constant a (a is a free parameter) ranges between and 1. The M matrix is the 
quark density matrix and is given by 



M^lik) = j d%e^^-y{P'\i^^^{-ay)4\{l-a)y)\P). (5.24) 

The index i = u,d runs on flavours. Using a Sudakov decomposition of the internal loop we 
can rewrite T^'^ and T^'^ as 



(5.25) 



(2vr) 

Tr {~gr {Uy - 7') pEr9uM^{k)+ 

~gr [Dy - 7') pEl^'gdM^k)] (5.26) 



to which we wiU refer as the direct and the exchange diagram respectively. It is also convenient 
to introduce two new linear combinations T^^ = + = Tj^ + T'^ which will turn useful 
in order to separate Vector (V) and axial vector parts (A) of the expansion 



~gguUyTr{[^^f>Dl^ + 7^/'s7l + 
~ggdDyTr[Yl^PDl'' + I^PeY] M\k)] 



(5.27) 



55,rr {[7Ti57^ + t'^/'etI 7'M"(^)} + 

-ggaTr { + T^/'etI l^'K^k)] . (5.28) 

with Ta including the vector parts iV xV + ^4 x ^) and Tb the axial-vector parts 
iy X A + AxV). After some algebraic manipulations we finally obtain 



rpl-LV 



T.99^C^ [ ^e*^4(W + W-5^na(z)(P'|V' 



IP) 



+ze''-/^n«n,/3(.)(P'|V^« (-^) 7V ^« (^) l^')} 



i=u,d 



dXdz 



(2vr) 



2 / 



An 



+ze^'^'*^n,n^/3(z)(P'|V 



"2"; 



IP) 



(5.29) 



where Ci = Uy, and 



a{z) = + , I3{z) = : (5.30) 

are the "first order" propagators appearing in the factorization of the amphtude. We recall, if 
not obvious, that differently from DIS, DVCS undergoes factorization directly at amplitude level 
[106]. 

The parameterizations of the non- forward light cone correlators in terms of GPD's is of the 
form given by Ji at leading twist |100j 

H{z, e, ^')U{P')rU{P) + E{z, e, A')U{P')^^U{P) + 



H{z,C,A^)U{P')r7'U{P)+E{z,C,A^)U{P')-L—U{P) + (5.31) 

which have been expanded in terms of functions H, E, H, E |127| and the ellipses are meant 
to denote the higher-twist contributions. It is interesting to observe that the amplitude is still 
described by the same light-cone correlators as in the electromagnetic case (vector, axial vector) 
but now parity is not conserved. 



5.4 Operatorial analysis 

The operatorial structure of the T-order product of one electroweak current and one electro- 
magnetic current is relevant in order to identify the independent amplitudes appearing in the 
correlator at leading twist and the study is presented here. We will identify four operatorial 
structures. For this purpose let's start from the Fourier transform of the correlator of the two 
currents 

T^, = iJ <fxe"'^P2\T (jJ(x/2) (-x/2)) , (5.32) 
where for the neutral and electromagnetic currents we have the following expressions 

J^^«(-x/2) = -^^^J-x/2)7^(5fy + 5fA7')V'.(-^/2) +?,(-x/2)7'^(5|y +5|47')V'd(-x/2), 
I cos 0\Y 

r'^{x/2) = ^Ax/2)Y (-^e) + ^ui^mY (^e) Mx/2) . (5.33) 

^ Based on the article published in JHEP 0502, 038 (2005) 



By simple calculations one obtains 



(P2|r(JJ(x/2)jf (-x/2))|Pi) = 

(P2|V^Ja;/2)g„7,5(x)7^((7fy + <7fA7')V'n(-a;/2) - 
^d{x/'2)9dluS{x)-i^,{g^y + g^A-i^)^IJd{-x/2) + 
^u{x/2)-i^{g^y + g^j^-i^)S{-x)gulu'4)u{x/2) - 
lPa{x/2)j^{gly + 5|475)S(-x)5d7.Vd(a;/2)|Pi 



(5.34) 



The coefficients used in eqs. ()5.33l \5^Mh gy and g^, are 



9dV 



— I — sin 
2 3 



- + 2 sin Ow 



and 



9u — 2 ) 



1 

9d-^ 



9uA 



9dA 



1 

~2 
1 

2 ^ 



(5.35) 



(5.36) 



are the absolute values of the charges of the up and down quarks in units of the electron charge. 
The function S{x) denotes the free quark propagator 



S{x) 



27r2(x2 -ie)2' 



(5.37) 



After some standard identities for the 7's products 



S^aul3 = {9lJ.agvP + 9iyagiil3 — 9iJ.ugaf3) 



(5.38) 



we rewrite the correlators as 



9u9uV I J^iaufS^u ~ 1'^tJLOivl3'^u 



q r)5/3 

>J iirvijf-i^ti 



a'? 



—9d9dV [Sfiauf^O^ — ie^aulsO^f) + 9dgdA [S^aulsO^f — i^/iaulsO'^ 



\Pl)- 

(5.39) 



The x-dependence of the operators in the former equations was suppressed. 



Whence the relevant operators are denoted by 



6f(x/2,-x/2) = V^,(x/2)7'3V'„(-x/2) +V^,(-x/2)7'3^„(x/2), 
Of (x/2,-x/2) = V^,(x/2)757^^,(-x/2) -V^,(-x/2)7S'3Va(a;/2), 
Of(x/2,-x/2) = V^,(x/2)7''V^.(-.x/2) -?,(-x/2)7''^a(x/2), 
Of (x/2, -x/2) = V^„(x/2)757/^Va(-x/2) +liia{-^/2)^''i^^a{x/2) , 

(5.40) 

where a is a flavour index. 



5.5 Phases of the Amphtude 



The numerical computation of the cross section requires a prescription for a correct handling 
of the singularities in the integration region at z = The best way to proceed is to work 
out explicitly the structure of the factorization formula of the amplitude using the Feynman 
prescription for going around the singularities, thereby isolating a principal value integral (P.V., 
which is real) and an imaginary contribution coming from the d function term. A P.V. integral 
is expressed in terms of "plus" distributions and of logarithmic terms, as illustrated below. The 
expression of the factorization formula of the process in the parton model, in which a{z) and 
(3{z) appear as factors in the coefficient functions, is then given by 

Mfi = 4iqi)D{qi)e''*iqi-A) 



a{z) 



H^{z, A2)i7(P2)# C/(Pi) + E^{z, A^)f/(P2 



fjL'—^V 



2M 



U{Pi 



+ 



P{z)ie^^-^han(, [H^iz, ^, A^)U{P2)i l'U{P,) + E-{z, A'')U{P2h' (A • n) C/(Pi)' 
-^g9dDv J dz{u^ d}- 



+ 



^99u j\dz{-n^'rf -h^ni' + gn 



a{z) 



H^{z, A^)U{P2)i l^'UiPi) + E^{z, A^)U{P2 



i7^A • n 
' 2M 



U{Pi) 



+ 



(5{z)ie>''"^Pne,nfi 

' J dz {u d}^ . 



H^iz, A2)t/(P2)yi UiPi) + E^iz, t A^)U{P2) " (Pi] 



2M 



:99d 



(5.41) 



To handle the singularity on the path of integration in the factorization formula, as we have 
already mentioned, we use the Feynman (ie) prescription, thereby generating imaginary parts. 
In particular, any standard integral containing imaginary parts is then separated into real and 
imaginary contributions as 

dz ^^^] = PV C dz'^^ T iT^Ti±0 (5.42) 
for a real coefficient T{z). We then rewrite the P.V. integral in terms of "plus" distributions 



= dzQ{z)H{z) + J'^ dzQ{z)H{z) + H{0 log (^^ 



(5.43) 



where 



Q{z) = e{-i<z<o ^ 



Q{z) = e{^<z<i) ^ 



(5.44) 



and the integrals are discretized using finite elements methods, in order to have high numer- 
ical accuracy. This last point is illustrated in Appendix C, where the computations are done 
analytically on a grid and then the grid spacing is sent to zero. 

We can now proceed and compute the cross section. We define the scalar amplitude 

Mfi = J^^,^{q,)D{q{)T>^'^e*^{q, - A) (5.45) 

where D{qi) is the Zq propagator in the Feynman gauge and Jif,p{qi) is the leptonic current and 
we have introduced the polarization vector for the final state photon . 

In particular, for the squared amplitude we have 

\Mfi\^ = -L^'^D{qifT^,Tl^ (5.46) 
which is given, more specifically, by 

\M\^ = dz dz' {Ki{z,z')a{z)a*{z')+K2{z,z')P{z)P*{z')) (5.47) 



with Ki and K2 real functions, combinations of the generahzed distributions {H, H, E, E) with 



appropriate kinematical factors. Mixed contributions proportional to a{z)l3*[z') and f){z)a*{z') 
cancel both in their real and imaginary parts and as such do not contribute to the phases. A 
similar result holds also for the pure electromagnetic case. 

After some further manipulations, we finally rewrite the squared amplitude in terms of a P.V. 
contribution plus some additional terms coming from the imaginary parts 



which will be analized numerically in the sections below. In order to proceed with the numerical 
result, it is necessary to review the standard construction of the nonforward parton distribution 
functions in terms of the forward distributions, which is the topic of the next section. 



5.6 Construction of the Input Distributions 

The computation of the cross section proceeds rather straightforwardly, though the construction 
of the initial conditions is more involved compared to the forward (DIS) case. This construction 
has been worked out in several papers [1071 11081 11091 IllOl lllll 11021 1108j in the case of standard 
DVCS, using a diagonal input appropriately extended to the non-diagonal kinematics. Different 
types of nonforward parton distribution, all widely used in the numerical implementations have 
been put forward, beside Ji's original distributions, which we will be using in order to construct 
the initial conditions for our process. 

For our purposes it will be useful to introduce Golec-Biernat and Martin's (GBM) distribu- 
tions |109| at an intermediate step, which are linearly related to Ji's distributions. 

We recall, at this point, that the quark distributions Hq{z^£) have support in z G [—1,1], 
describing both quark and antiquark distributions for z > and z < respectively. In terms of 
GBM distributions, two distinct distributions J^^{X, and J^''{X, () with < X < 1 are needed 
in order to cover the same information contained in Ji's distributions using only a positive scaling 
variable (X). In the region X G {(, 1] the functions T'^ and J^^ are independent, but if X < ( 
they are related to each other, as shown in the (by now standard) plot in Fig. 15.41 

In this new variable (X) the DGLAP region is described hy X > ( {\z\ > (^), and the ERBL 
region hy X < ( {\z\ < (^). In the ERBL region, J^'^ and are not independent. 

The relation between H{z,£,) and jF''{X,C,) can be obtained explicitly [112] as follows: for 

z G {—£,-, 1] we have 




P.V. dz dz' {Ki{z,z')a{z)a*{z')+K2{z,z')/3{z)/3*{z')) 




(5.48) 




(5.49) 



DGLAP 



ERBL 



(1+ k) 



X 



- -1 



(1+^) 



ERBL 



DGLAP 



Figure 5.4: The relationship between JP"^(X, C), T'^{X,C) and Ji's function H''{z,£). 



and for z £ [— 1, 



where i is a flavour index. In our calculations we use a simplified model for the GPD's where 
the dependence can be factorized as follows |113l I107| 

E%z,i,A\Q^) = F^{A')r\z,^,Q^) (5.51) 

where q^{z) and Aq^{z) are obtained from the standard non-polarized and longitudinally polarized 
(forward) quark distributions using a specific diagonal ansatz |114| . The ansatz r*(z, ^) = q^{z, ^) 
is also necessary in order for the quark sum rule to hold |115| . Analogously, in the case of the 
distributions |102l 11161 1117] one can use the special model 

E- = E'' = -e{i - |z|)0^(z/e)5.(A2), 5,(A2) = y </.^(x) = -(1 - x2) 

IE, — A^ 3 

(5.52) 

valid at small A^, where g^J^ = 1.267, M is the nucleon mass and m.„^ is the pion mass, with the 
normalization 

F{(0) =Gi(0) = 1. (5.53) 



Notice that, analogously to the H distributions, q^{z,E,,Q'^) and Ag*(z,^,Q^), which describe 
the A^ = limit of the H-distributions, have support in [—1, 1] and, again, they describe quark 



distributions (for z > 0) and antiquark distributions (for z < 0) 



cfiz,^,Q^) = -q\-z,^,Q^) 

Aq\z,C,Q^) = Aq\-z,tQ^). (5.54) 

Now we're going to estabilish a connection between the q{z,(,,Q'^) and the j^'^{XX) functions, 
which is done using Radyushkin's nonforward "double distributions" |102j . The construction of 
the input distributions, in correspondence of an input scale Qo, is performed following a standard 
strategy. This consists in generating nonforward double distributions f{x,y) from the forward 
ones {f{x)) using a "profile function" Tr{x,y) [101] 

f{y,x) =7r{y,x)f{x), (5.55) 

where we just recall that the 7r{y, x) function can be represented by 

taken to be of an asymptotic shape (see ref. |101l Il08| ) for quarks and gluons. A more general 
profile is given by 

_ r(26 + 2) [{1 - \x\f - y^]" 

na^'y^- 226+lr2(6+l) (l-|x|)2^>+l ^^-^^^ 

and normalized so that 



dy tt{x, y) = 1 . (5.58) 

1+1x1 

b parameterizes the size of the skewing effects starting from the diagonal input. Other choices of 
the profile function are also possible. For instance, the double distributions (DD) defined above 
have to satisfy a symmetry constraint based on hermiticity. This demands that these must be 
symmetric with respect to the exchange y < — > 1 — x — y, and a profile function which respects 
this symmetry constraint is given by [llSj 

6y{l -x-y) 

^(^'?/)= (i-xY • ^^^^ 

This symmetry is crucial for establishing proper analytical properties of meson production am- 
plitudes. We will be using below this profile and compare the cross section obtained with it 



against the one obtained with (I5.57p . 

Now we are able to generate distributions q[z,^,Q'^) in the z variable at = 0, q{z,^,Q'^), 
by integrating over the longitudinal fraction of momentum exchange y characteristic of the double 
distributions 



q{z,tQ^)= dx' dy'5{x' + Cy'-z)f{y',x',Q^). (5.60) 

J-l J-l+\x'\ 




Figure 5.5: GPD's and Hd generated by the diagonal parton distribution with a profile 
function ([5371) at an initial 0.26 GeV^ 



Using (I5.60p and the expression of the profile functions introduced above, the GBM distributions 
are generated by the relation 



X 



/ 2, 



.^'^''^(X, C) = - dx'vr^ i^x', -{X -x')+x' -1] q%x') . 



(5.61) 



with a similar expression for the anti-quark distributions in the DGLAP region X > C (z < —(^) 

(5.62) 



-x+c 

^5,a(x, C) = ^ y dx'n" (x',-^iX + x') +x' + l 



r{\x'\ 



In the ERBL region, X < C, {\z\ < after the integration over y, we are left with the sum of 
two integrals 



c 



X 



J 2, 



dx'-K'i \yx', -{X -x')+x' -I] q^{x') 

dx'TTl fx', j{X -X')+X' -1] (f{\x'\ 

x-( V c 



dxV fx', --(X + x') + x' + 1 ) q'^ix'^ 





J dx'TT" fx', -^{X + x') + x' + ij g"(|x'|) 



(5.63) 



(5.64) 



Solving numerically the integrals we obtain the value of the function J-'^s on a grid, and 
using eqs. (|5.49l) and 15.501 we end up with the numerical form of the H-distributions. We have 
used diagonal parton distribution functions at 0.26 GeV^ [114| and the results of our numerical 



60 



Q =0.26 GeV^ 



40 - 




■20 - 



■40 - 



■60 



-0.04 



-0.02 







0.02 



0.04 



z 



Figure 5.6: GPD's flavour singlet combination at 0.26 GeV^ generated with a profile (|5.57l 




implementation can be visualized in Figs. 15.5 1 and 15. 61 



5.7 The Differential Cross Section 

Our kinematical setup is illustrated in Fig. 15.81 and we choose momenta in the target frame with 
the following parameterizations 

1 = {E,0,0,E) , I' = {E',E' cos(l)^sme^,E' sm(t)^sme^,E' cose^) , 
Pi = (M, 0,0,0) , P2 = (-6^2, 1-P2I cos(/)Arsin6lAr, IP2I sin(/)Arsin6lAr, IP2I cos6lAr) (5.65) 

where the incoming neutrino is taken in the positive z-direction and the nucleon is originally at 
rest. The first plane is identified by the momenta of the final state nucleon and of the incoming 
neutrino, while the second plane is spanned by the final state neutrino and the same l-axis. 
is the angle between the x direction and the second plane, while 4>n '^s taken between the plane 
of the scattered nucleon and the same x axis. We recall that the general form of a differential 
cross section is given by 




da 



1 



\Mfi\\27T)'6^'\l + Pl-P2-l' -q2) 



d^U d^P2 d^q2 



2ro(27r)3 2p0(27r)3 2gO(27r)3 



(5.66) 



and it will be useful to express it in terms of standard quantities appearing in a standard DIS 
process such as Bjorken variable x, inelasticity parameter y, the momentum transfer plus some 
additional kinematical variables typical of DVCS such as the asymmetry parameter ^ and A^. 




Figure 5.7: A pictorial description of the DVNS experimental setup, where the recoiled nucleon 
is detected in coincidence with a final state photon. 




We will be using the relations 



(5.67) 

in the final computation of the cross section. It is also important to note that has to satisfy 
a kinematical constraint 



One possible cross section to study, in analogy to the DVCS case [119], is the following 
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dxdQ^d\A'^\d(^r dxdyd\A^\d(l)r SvrQ^ 



1 + 



AM'^x^ 



(5.69) 



where 4>r is the angle between the lepton and the hadron scattering planes. To proceed we also 
need the relations 
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where x is given by 



X 



2 4^""" 2 2 



(5.70) 



(5.71) 



After some manipulations we obtain a simplified expression for |A^jj| , similarly to eq. ()5.47|1 

fl rl 



= I dz I dz' [Ai{z,z',x,t,Q'^)a{z)a*{z')+A2{z,z',x,t,Q^)p{z)p*{z')] 



X [-2Q2 (4m2 -t){x- 2f {x - 1) x^y + {t - AM^f {x - if xV 
+ Q\x-2)^{2-2y + y^)] 

X [2M^{Mz^ + Q^f{x-2)^[Q^{x-2f-{4M^-t){-l + x)x^fy^]~'^ (5.72) 



where Ai{z, z' ,x,t,Q'^) and A2{z, z' ,x,t,Q'^) are functions of the invariants of the process and 
of the entire set of GPD's. Their explicit form is given in Appendix D. As we have already 
mentioned, the {z,z') integration is be done by using the Feynman prescription to extract the 
phases and then using the distributional identities (|5.43|) and (|5.44l) . For numerical accuracy we 
have discretized the final integrals by finite element methods, as shown in Appendix C. 

We have plotted the differential cross section as a function of Q^, for various values of 
and at fixed x values. We have used both the profile given by (|5.57p (Figs. 15.9115. 12|) and the 
one given in eq. (I5.59|) (Figs. [5TT3] and ISTT^ and their direct comparison in a specific kinematical 
region (Fig. 15.151) . The different profiles generate differences in the cross sections especially for 
larger values. Notice also that the DVNS cross section decreases rather sharply with A^, 
at the same time it increases appreciably with x. The results shown are comparable with other 
cross sections evaluated in the quasi-elastic region (« 10~^ nb) for charged and neutral current 
interactions, and appear to be sizeable. Coherence effects due to neutral current interactions with 
heavy nuclei, in particular with the neutron component may substantially increase the size of 
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Figure 5.9: DVCS cross section at x = 0.1 and center of mass energy ME = 10 GeV^ using the 
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Figure 5.11: DVCS cross section at x = 0.3 and center of mass energy ME = 10 GeV^ using 
the profile ([5371) . 
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Figure 5.13: DVCS cross section at x = 0.3 and center of mass energy ME = 10 GeV^ with 
NPD functions generated by the profile function (|5.59p . 
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Figure 5.14: DVCS cross section at x = 0.3 and center of mass energy ME = 27 GeV^ with 
NPD functions generated by the profile function (|5.59p . 
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Figure 5.15: DVCS cross sections at x = 0.2 and center of mass energy ME = 27 GeV^. a using 
profile (I539D . 

the cross sections, with an enhancement proportional to N'^, where N is the number of neutrons 
[120| . though an accurate quantification of these effects requires a special study }12lj which is 
underway. It is worth to emphasize that in the past this contribution had never been included 
in the study of neutrino-nucleon interactions since very little was known about the intermediate 
energy kinematics in QCD from the point of view of factorization. It seems obvious to us 
that with the new developments now taking place in the study of QCD at intermediate energy, 
especially in the case of the generalized Bjorken region, of which the deeply virtual scattering 
limit is just a special case, it will be of wide interest to quantify with accuracy the role of these 
new contributions for neutrino factories. In general, one expects that electromagnetic effects are 
suppressed compared to the standard (hadronic) deeply inelastic cross section, and this has led 
in the past to a parameterization of the intermediate energy cross section as either dominated 
by the quasi elastic region and/or by the DIS region at higher energies, as we have mentioned 
in our introduction. However, the exclusive cross section has some special positive features, one 
of them being to provide a clean signal for the detection of weakly interacting particles, and we 
expect that this aspect is going to be of relevance at experimental level. 



5.8 Conclusions 



We have presented an extension of the standard DVCS process to the case of one neutral current 
exchange, describing the scattering of a neutrino off a proton in the parton model. We have 
described the leading twist behaviour of the cross section; we have found that this is comparable 
to other typical neutrino cross sections and discussed its forward or DIS limit. We have presented 
a complete formalism for the study of these processes in the parton model. The process is the 
natural generalization of DIS with neutral currents and relies on the notion of Generalized Parton 
Distributions, new constructs in the parton model which have received considerable attention 
in recent years. The possible applications of these new processes are manifold and we hope to 



return in the near future with a discussion of some of the issues not addressed in this thesis. 



5.9 Appendix A 



The colhnear expansion of the internal loop momentum k allows to identify the light cone oper- 
ators appearing in the process at leading twist. We recall at this point that the analysis of the 
hand-bag contribution is carried out exactly as in the electromagnetic case. 

To perform the collinear expansion and isolate the light-cone correlators of DVNS from the 
hand-bag contribution we use the relation 
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inside the expression of T^^'^ in order to obtain 
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Keeping the leading terms for the loop momenta 

fc'^-aA'^ + gf = n^'{k■n + 2a^-2Cj + n^'(^■n-a(M'^ + 
ki" + {I - a)A^ - = n'*(/c • n- 2(1 - a)4 + 20 + (^fe-n + 2(1- a)CM^-^j 
we thus obtain 
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We expand 



k ■ z = (k ■ n) {n ■ z) + (k ■ h) {n ■ z) — k± ■ z± 



and choose a = 1/2. These expansions are introduced in eg. 15.211 and after some manipulations 
the tensor T now becomes 
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We also recall that the expansion of the matrix element M_{k) proceeds also in this case as in 
the electromagnetic case 

M^,{k) = fd^ye'^-y{P'\iif(-ay)4\{l-a)y)\P) = A,^+A2^,^+... (5.74) 



where the ellipses refer to terms which are of higher twist or disappear in the trace of the diagram. 



5.10 Appendix B 

The last ingredients needed in the construction of the input distribution functions are the form 
factors Fl and Fg. From experimental measurements we know, by a dipole parametrization in 
the small region, that 

^^(A^) = (1 + k,)-1G^(A2) = n-^GU^^) = {l-K] , G^(A^) = 0, (5.75) 

where the electric, G^(A^) = F{(A^) + ^^^2(^^)1 magnetic form factors G^j(A^) = 
F{(A^) + F|(A^) are usually parametrized in terms of a cutoff mass my = 0.84 GeV. 
For non-polarized GPD's the valence u and d quark form factors in the proton can be easily 
extracted from F^"'' = 2 Ff + {^^"^^ Ff and given exclicitely by 

2F/(A2) = 2F|'(A2) + FP(A2), F/(A2) = ^^(A^) + 2FP(A2), for / = 1, 2. (5.76) 



This exploits the fact that proton and neutron form an iso-spin doublet. 



At the scale uia = 0.9 GeV one can get 

G\{A')=(l-^\ (5.77) 
for the valence quarks. For the form factors F we obtain 



where 
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' "m2(i-ba^)^(-i + ^) a -bay • 



(1 - BAY (l - iJ.) 



(1 - BAY (l - 



(5.78) 



A = 0.238 5 = 1.417 C = 0.447, 

£) = 0.835 S = 0.477 F = 0.120. (5.79) 



5.11 Appendix C 



In this section we illustrate an analytical computation of the integrals by discretization, using 
finite elements method. We want to approximate with high numerical accuracy integrals of the 
form 
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For this purpose we start by choosing a grid on the interval (—1 = xq, ,Xn+i = define 

J. = f dz^M^ = t r d.'M^ . (5.81) 
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Performing a simple linear interpolation we get 
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(5.83) 



Now, moving to the integral in the interval (^, 1), we introduce a similar grid of equally spaced 
points (^ = yo, ,yn+i = 1) and define the integral 
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As above, after isolating the singularity we obtain 
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Again, performing the integrations we obtain 
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(5.86) 



Collecting our results, at the end we obtain 
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We can use the same strategy for the integrals of "+" type defined as follows 
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This time we call our final integrals Xi and X2. They are given by the expressions 
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with a discretization supported in the (—1 = xq, ,Xn+i = grid, and 
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on the (— ^ = yo) )J/n+i = 1) grid. As a final result for the "+" integral we get 
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5.12 Appendix D 



In this section we will present the full expression of the functions Ai and A2 which appear in 
the squared amplitude 

Aiiz,z'x,t,Q^) = g^Q^ [4gj[E'a{4HdM^ + Edt)x^ 

+ AH'^M\AHd{x - I) + Edx^)] 

+ Agdgu[{AE',HdM'' + AE'^HuM^ + KEdt + E'dEut)x^ 

+ AH'^M\AHd{x - 1) + EdX^) + AH'^M^AH^ix - 1) + Kx^)] 

+ D^U^gdgu[AE'^Edt + AE'^Eut - AE'^Eatx - AE'^E^tx + AE'^EdM^x^ 

+ AE'^EuM^x^ + AE'^HdM^x^ + AE'^HuM^x^ + E'^Edtx^ + E'^E^tx^ 

+ AH'^M\AHd{x - 1) + EdX^) + AH'^M\AHu{x - 1) + E.x^)] 

+ 5^[4£;;£;^t[72 _ 4E'^EutU^x + IGE'^HuM^x^ + AE'^Eutx^ 

+ 4E;£;„M2c/Jx2 + AE'^H^M^U^x'^ + E'^EJU^x"^ + im'^M'^{AH^{x - 1) + ^x^) 
+ AH'^M^U^{AH^{x - 1) + 

+ L>2^^[4iJ^M2(4/fd(x - 1) + ^^dx^) + E'aiAHdM^x'' + ^^^(t (x - 2)^ + 4M2x2))]] 

(5.92) 

and for A2{z, z' ,x,t) we get a similar result 

A2{z, z', X, t, Q2) = ig^Q^ ^gdgu[AE'^Edt + AE'^EJ - 16D,H'^HdM^U, - 16D,H'dHuM^U, - AE'^Edtx 
- AE'^Eutx + IGD^HuHdM^U^x + I&D^H'^HuM'^U^x + AE'^EaM^x^ + AE'^E^M^x^ 
+ AE'^HdM^x^ + AE'dHuM^x^ + E'^Edtx^ + E'dEutx^ + AD^EuHdM^U^x"^ 
+ AD^EdH'^W^U-^x^ + AD^E'^HdM^U^x^ + AD^E'dHuM^U^x^ + D^E'^EdtU^ 
+ D^E'^E^tUy + AH',M\AHdix - 1) + + AH'^M\AH^{x - 1) + 

+ 5d[4i^^M2(4ifd(x - 1) + SdX^) + Dl{E'^{AHdM'' + E^^tjx^ 
+ AH'^M^AHdix - 1) + ^dx^)) + E'^iAHdM^x^ + Ed{t{x - 2f + AM'^x''))] 
+ gl[AH'^M\AH^{x - 1) + ^^^x^) + {E'^{4H^M'' + E^t)x'' 
+ AH'^M'^{AHu{x-1)+Eux'^)) + E'^{AHuM'^x'^ + Eu{t{x-2f + AM'^x'^))\ . 

(5.93) 



Chapter 6 



Leading Twist Amplitudes for Exclusive 
Neutrino Interactions in the Deeply 
Virtual Limit 

6.1 Introduction and Motivations 



In the prevoius chapter we have pointed out |122| that exclusive processes of DVCS-type (Deeply 
Virtual Compton Scattering) could be relevant also in the theoretical study of the exclusive 
neutrino/nucleon interaction. Thanks to the presence of an on-shell photon emitted in the final 
state, this particle could be tagged together with the recoiling nucleon in a large underground 
detector in order to trigger on the process and exclude contamination from other backgrounds. 
With these motivations, a study of the uN uN^ process has been performed in |122| . The 
process is mediated by a neutral current and is particularly clean since there is no Bethe-Heitler 
contribution. It has been termed Deeply Virtual Neutrino Scattering or DVNS and requires 
in its partonic description the electroweak analogue of the "non-forward parton distributions", 
previously introduced in the study of DVCS. 

In this section we extend that analysis and provide, in part, a generalization of those results 
to the charged current case. Our treatment, here, is purposely short. The method that we use for 
the study of the charged processes is based on the formalism of the non-local operator product 
expansion and the technique of the harmonic polynomials, which allows to classify the various 
contributions to the interaction in terms of operators of a definite geometrical twist |123| . We 
present here a classification of the leading twist amplitudes of the charged process while a detailed 
phenomenological analysis useful for future experimental searches will be given elsewhere. 
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6.2 The Generalized Bjorken Region and DVCS 



Fig. 15.11 illustrates the process that we are going to study, where a neutrino of momentum I 
scatters off a nucleon of momentum Pi via a neutral or a charged current interaction; from the 
final state a photon and a nucleon emerge, of momenta q2 and P2 respectively, while the momenta 
of the final lepton is I' . We recall that Compton scattering has been investigated in the near 
past by several groups, since the original works |100l llOll I102| . A previous study of the Virtual 
Compton process in the generalized Bjorken region, of which DVCS is just a particular case, 
can be found in |103| . From the hadronic side, the description of the interaction proceeds via 
new constructs of the parton model termed generalized parton distributions (GPD) or also non- 
forward parton distributions. The kinematics for the study of GPD's is characterized by a deep 
virtuality of the exchanged photon in the initial interaction {u+p — > v + p + 'y) ( ~ 2 GeV^), 
with the final state photon kept on-shell; large energy of the hadronic system [W"^ > 6 GeV^) 
above the resonance domain and small momentum transfers \t\ < 1 GeV^. In the electroweak 
case, photon emission can occur from the final state electron (in the case of charged current 
interactions) and provides an additional contribution to the virtual Compton amplitude. We 
choose symmetric defining momenta and use as independent variables the average of the hadron 
and gauge bosons momenta 

Pi,2 = PT^ qi,2 = q±^, (6.1) 
with A = P2 — Pi being the momentum transfer. Clearly 

A2 

P • A = 0, p2 = - — (6.2) 

and M is the nucleon mass. There are two scaling variables which are identified in the process, 
since 3 scalar products can grow large in the generalized Bjorken limit: , A ■ q, P ■ q. 

The momentum transfer A^ is a small parameter in the process. Momentum asymmetries 
between the initial and the final state nucleon are measured by two scaling parameters, ^ and rj, 
related to ratios of the former invariants 

e = -^ V = ^ (6.3) 
^ 2P-q ' 2P-q ^ ^ 

where ^ is a variable of Bjorken type, expressed in terms of average momenta rather than nucleon 
and gauge bosons momenta. The standard Bjorken variable x = —qf/{2Pi -qi) is trivially related 
to in the t = limit and in the DVCS case rj = 

Notice also that the parameter ^ measures the ratio between the plus component of the 
momentum transfer and the average momentum. 

S,, therefore, parametrizes the large component of the momentum transfer A, which can be 
generically described as 

A = -2^P - A^ (6.4) 



where all the components of Aj^ are O ( yjA^ 



6.3 Bethe-Heitler Contributions 



Prior to embark on the discussion of the virtual Compton contribution, we quote the result for 
the Bethe-Heitler (BH) subprocess, which makes its first appearance in the charged current case, 
since a real photon can be radiated off the leg of the final state lepton. The amplitude of the 
BH contribution for a exchange is as follows 
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where e is the polarization vector of the photon and 
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(6.6) 



for the W case, with D^^{qi)/{A'^ — + ie) being the propagator of the W's and Fi^2 the 
usual nucleon form factors (see also tl22j ). 



6.4 Structure of the Compton amplitude for charged and neutral 
currents 



Moving to the Compton amplitude for charged and neutral currents, this can be expressed in 
terms of the correlator of currents 



V = i| dW^-(P2|T (jJ(x/2)j;^*'^o(-x/2)) , 



(6.7) 



where for the charged and neutral currents we have the following expressions 

J'^^-(-x/2) = -^V^,(-x/2)7^(l - ^')UduM-x/2), 

r^^x/2) = ^a{x/2)Y (-\^ +^^^^^1" {le\ V«(x/2) . 



Here we have chosen a simple representation of the flavour mixing matrix U*^ = Uud = Udu = 
cos9c, where 6c is the Cabibbo angle. 

The coefficients gy and g^ are 

z 1 4 . 2 ^ z 1 

9uV=2^^ ^™ 9uA = "2 

9dV = -\ + \ sin^ Ow 9dA = \, (6-9) 

and 

2 1 , , 

are the absolute values of the charges of the up and down quarks in units of the electron charge. 
A short computation gives 

(P2|r(jJ(x/2)jJo(-x/2))|Pi) = 
(P2K(x/2)5„7,5(x)7^(5fy +5f^75)V.„(-,T/2) - 
V^rf(x/2)<7rf7,5(x)7^(5i. +5|475)^rf(-x/2) + 

^ui^f^Muv + guAl'')S{-x)guivMx/'2) - 
^d{x/2)ji,{gly + g^Al'')S{-x)gdlMx/2)\Pi), 



(6.11) 



(P2|r(jJ(x/2)J^^+(-x/2)) |Pi) = 
(P2|V^„(-x/2)74l - ^^)U^dS{-x)j, i-ga) M^/2) + 
V^Jrr/2)7, {g^) S{x)^^{l - j^)U^dM-^m\Pi), 

(6.12) 

(P2|t(jJ(x/2)J;^-(-x/2)) |Pi) = 

(P2I -?d(x/2)5d7,5(x)7^(l - j')UduM-x/2) + 

V^d(-x/2)7^(l - j')Si-x)UduMxm\Pi) , (6.13) 



where all the factors g/2-\/2 and g/2cos 9w, for semplicity, have been suppressed and we have 



defined 
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Using tiie following identities 
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we rewrite the correlators as 
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(6.18) 



We have suppressed the x-dependence of the operators in the former equations. The relevant 
operators are denoted by 



0^{x/2, -x/2) = Vjx/2)7'^V'a(-.x/2) + t/j,i-x/2)j^M^/2), 
Of (x/2,-x/2) = V^,(x/2)7S'^Va(-x/2) - V^„(-x/2)7S'^^„(x/2), 
Of (x/2, -x/2) =^„(x/2)7'3^„(-x/2) -^,(-x/2)7'^^a(x/2), 

ol''{x/2,-x/2) = Mx/2h'7''M-^/'^) + M-^/'^h'i''M^m, 

df^(x/2, -x/2) = gu^^{x/2)-f''M-x/2) + 5d?n(-x/2)7'3Vd(x/2), 



(6.19) 



0^;^(x/2, -x/2) = 5„V„(x/2)7^^^d(-x/2) - ^^,V„(-a:/2)7^^^d(x/2) 
Of^(x/2, -x/2) = 5„V^Jx/2)7^Vd(-x/2) - <7rfV^J-x/2)7'3v.,(x/2), 
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0;:^(x/2, -x/2) = 5„^„(x/2)757^^rf(-x/2) + J-a^/2)7^^V'd(^/2) 



(6.20) 



and similar ones with u ^ d interchanged. 

We use isospin symmetry to relate flavour nondiagonal operators {Off) to flavour diagonal 



ones (Off) 
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(n|d'^"(x)|p) = (n|d'^°'(x)|n) - (n|d""(x)|n), 



(6.21) 



where 



(6.22) 



are isospin raising/lowering operators expressed in terms of Pauli matrices. 



6.5 Parameterization of nonforward matrix elements 

The extraction of the leading twist contribution to the handbag diagram is performed using the 
geometrical twist expansion, as developed in [iMl [ESI [Ml [l27l |I28l llMl [l30l [H^ adapted 
to our case. We set the twist-2 expansions on the light cone (with x'^ = 0) and we choose the 
light-cone gauge to remove the gauge link 
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with < < 1 a scalar parameter, with 



P, = P,Zl + P2Z2, 



(6.24) 



and 
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(6.25) 



The index {u) in the expressions of the distribution functions F, G has been introduced in order 
to distinguish them from the parameterization given in |124l I103| . which are related to linear 



combinations of electromagnetic correlators. In the expressions above o is a flavour index and we 
have introduced both a vector (Dirac) and a Pauli-type form factor contribution with nucleon 
wave functions (U(P)). The product Pi ■ Pj denotes all the possible products of the two momenta 
Pi and P2, and the measure of integration is defined by |124| 

Dz = ^dzidz2 9{1 - zi) e{l + zi) e{l - Z2) Q{\ + Z2). (6.26) 

In our parameterization of the correlators we are omitting the so called "trace-terms" (see 
ref. [103]), since these terms vanish on shell. In order to arrive at a partonic interpretation 
one introduces variables z+ and z_ conjugated to 2P and A and defined as 

z+ = l/2(zi +Z2), 

Z^ = 1/2(Z2 - Zl), 

Dz = dz+dz^9{l + Z+ + z^)9{l + z+ - z^)e{l - z+ + z^)e{l - z+ - z_). (6.27) 

In terms of these new variables Pz = 2Pz+ + Az_, which will be used below. At this stage, we 
can proceed to calculate the hadronic tensor by performing the x-space integrations. This will 
be illustrated in the case of the current, the others being similar. We define 

/piqx a 
dx \ . (P2|5,..;30°^|Pl) = ~S%, 

zir^ [x^ — le) 

/piqx a 

dx \ . {P2\e,^,^0'^P\Pr) = 

ivr^ (x^ — le) 

dx \ .^2 (-P2lwO'"^l-Pi) = ^t- (6-28) 

27r^ (x^ — le) 

and introduce the variables 

Q2iz) = g° - ^P^, (6.29) 

where (z) is now meant to denote both variables (z+jZ^). The presence of a new variable Q2, 
compared to fl24j . is related to the fact that we are parameterizing each single bilinear covariant 
rather then linear combinations of them, as in the electromagnetic case. 



After some re-arrangements we get 

-Jgi:^ [Pzi^Pzu + Pzt^Qu + Pzyq^, " g^,u{Pz " ?)] | + 



{Qi + ie) 
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p '± _|_ p E. 



M(Q2 + ie) 
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9d / ^^77^;2 



M 



M 



M 



M 



+ 



M{Q'^ + ie) ^ P^^^P^'^ ~'~ -^^m^^^ + PzuQ/i 9tiu{Pz • q)\ \ , 



(6.30) 



with an analogous expressions for S^, that we omit, since it can be recovered by performing the 
substitutions 



paly) Qa{v) ^ pba{v) Qba{v) 



(6.31) 



in (lOnll . 

Similarly, for we get 



e-, = 5n / DzF<^\z) 
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+ ie) 
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+ ie) 
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+ ie) 
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+ ie) 
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M M{Qj + ie) 

M M{Ql + ie) 



. (6.32) 



The expression of can be obtained in a similar way. To compute the tensor we need to 



include the following operators, which are related to the previous ones by , — > 1 



dx 



27r2 (x^ — ie 



-2{P2\Sf_iau(30"''^\Pl) — S'^j, 



27r^ (x^ — ie 
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4 c X /Die _ /o5a/3 1 d \ _ cSa 



dx' 
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dx' 



271^ (x^ — ie) 



27r2 (x^ — ie)^ 
In this case a simple manipulation of (I6.3QP gives 



{P2Ka.(,0''^^\Pi)=ell. 



(6.33) 
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(6.34) 



The expression of is obtained from (|6.34|) by the replacements (|6.31|) 
For the case, a re-arrangement of (|6.32l) gives 

PSq'^h 
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(6.35) 



Also in this case, the expression of the tensor is obtained by the replacements (|6.31|1 . 



6.6 The partonic interpretation 



At a first sight, the functions F^'^\ G^'^\F^^'^\ Q^i'^) do not have a simple partonic interpretation. 
To progress in this direction it is useful to perform the expansions of the propagators 



1 1 1 



Ql + ie 2{P ■q)[z+-C + r]Z- + ie] ' 

Ql + ie ^ ~2{P-q) [z+ + ^ + r]z^ - ie] 

which are valid only in the asymptotic limit. In this limit only the large kinematical invariants 
and their (finite) ratios are kept. In this expansion the physical scaling variable ^ appears quite 
naturally and one is led to introduce a new linear combination 

t = z+ + riz-, (6.37) 

to obtain 

1 1 1 



Q\ + ie 2{P ■q)[t-C + i^V 

111 
QlTi'e'^~2{P-q) [t + C-ie]- ^^'^^^ 

Analogously, we rewrite Pz using the variables {t,Z-} 

Pz = 2Pt + TTz. , (6.39) 

in terms of a new 4- vector, denoted by tt, which is a direct measure of the exchange of transverse 
momentum with respect to P 

7r = A + 2^P. (6.40) 

This quantity is strictly related to A_|_, as given in (16. 4p . The dominant (large) components of 
the process are related to the collinear contributions, and in our calculation the contributions 
proportional to the vector tt will be neglected. This, of course, introduces a violation of the 
transversality of the process of 0(A_|_/2P • q). 

Adopting the new variables {t, z^} and the conjugate ones {P, tt}, the relevant integrals that 



we need to "reduce" to a single (partonic) variable are contained in the expressions 
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(6.41) 



Here H{zj^, 2:_) is a generic symbol for any of the functions. We have similar expressions for the 
integrals depending on the momenta Q2. 

The integration over the Z- variable in the integrals shown above is performed by introducing 
a suitable spectral representation of the function H{t,+^Z-, Z-). As shown in |124]. we can 
classify these representations by the n = 0, 1, powers of the variable Z- , 



With the help of this relation one obtains 

1 



dz^z^^h{ — h ^z^,z- 

T 



(6.42) 



Hn{t,^) = - I dz^z^^H{t + iz^,z^) 
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dT 



signit) 

T 



A-"/in(A,e). 



(6.43) 



The result of this manipulation is to generate single- valued distribution amplitudes from double- 
valued ones. In the single-valued distributions hn{t,^) the new scaling variables t and ^ have 
a partonic interpretation. ^ measures the asymmetry between the momenta of the initial and 
final states, while it can be checked that the support of the variable t is the interval [—1,1]. 
The twist-2 part of the Compton amplitude is related only to the n = moment of Z- . Before 
performing the Z- integration in each integral of Eq. (|6.41|) using Eq. (|6.43l) - a typical example 
is Hq'^(^) - we reduce such integrals to the sum of two terms using the identity 
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^ m— 1 
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(6.44) 



As shown in |103| . after the Z- integration, the integrals in (|6.41|) can be re- written in the form 



ffQ,({) 



2P-qJ_^ {t-C + ie) ' 



"Q^^^^- {2P.qfj_, 
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(6.45) 



where, again, we are neglecting contributions from the terms proportional to ir^, subleading in 
the deeply virtual limit. The quantities that actually have a strict partonic interpretation are the 
^o(^)0 functions, as argued in ref. |133| . The identification of the leading twist contributions is 
performed exactly as in [124j . We use a suitable form of the polarization vectors (for the gauge 
bosons) to generate the helicity components of the amplitudes and perform the asymptotic 
(DVCS) limit in order to identify the leading terms. Terms of 0{l/\/2P ■ q) are suppressed and 
are not kept into account. Below we will show only the tensor structures which survive after this 
limit. 



6.7 Organizing the Compton amplitudes 



In order to give a more compact expression for the amplitudes of our processes we define 

^ ^ +(g-P) + (g-P)' 

/,\ 9u ffd 
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gu , gd 



Pit) 

Calculating all the integrals in the Eqs 
sions of the amplitudes as follows 
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(6.46) 
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. (|6.30|) . (|6.34p and (I6.32p . (|6.35|1 . we rewrite the expres- 
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(6.47) 



where, suppressing all the subleading terms in the tensor structures, we get 
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while for the e"'^'^ expression we obtain 
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(6.49) 



Passing to the S"'^'^ and e"'^^ tensors, which appear in the Zq neutral current exchange, we get 
the following formulas 
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(6.51) 



Obviously the §"■^1^'^^ ga5tj.u^ ^abfiu ^^^^ ^ab^v expressions are obtained by the substitution (I6.3ip . 

At this stage, to square the amplitude, we need to calculate the following quantity, separately 
for the two charged processes 



= l-l-DVNSl + J-DVNSJ^ BH + J^BHJ^ DVNS + KBi^l 



(6.52) 



which simplifies in the neutral case, since it reduces |TDy7V5p |122| . Eqs. (I6.47p - (|6.51l) and their 
axial counterparts are our final result and provide a description of the deeply virtual amplitude 
in the electroweak sector for charged and neutral currents. The result can be expressed in terms 
of a small set of parton distribution functions which can be easily related to generalized parton 



distributions, as in standard DVCS. 



6.8 Conclusions 

We have presented an application/extension of a method, which has been formulated in the 
past for the identification of the leading twist contributions to the parton amplitude in the 
generalized Bjorken region, to the electroweak case. We have considered the special case of a 
deeply virtual kinematics. We have focused our attention on processes initiated by neutrinos. 
From the theoretical and experimental viewpoints the study of these processes is of interest, since 
very little is known of the neutrino interaction at intermediate energy in these more complex 
kinematical domains. 



^ Based on the article il34i 
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